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UNSTEADY  AERODYNAMICS  SUMMARY 


Final  report  for  Work  Unit  2300/FF/07 

This  task  addressed  basic  research  to  better  understand  unsteady,  separated  flows  and  to 
utilize  this  understanding  to  develop  innovative  concepts  to  control  unsteady 
aerodynamic  phenomena  to  enhance  flight  vehicle  maneuverability.  This  effort  sought  to 
expand  the  flight  performance  boundaries  by  developing  a  capability  to  maneuver  in  the 
high-angle-of-attack,  post-stall  flight  regime.  Research  directly  supported  Air  Force 
objectives  for  unsteady  and  separated  flows,  as  outlined  in  the  current  research 
Technology  Area  Plan  for  Project  2307,  Fluid  Mechanics. 

The  research  was  organized  around  two  broad  topical  areas:  (1)  understanding  and 
predicting  the  dynamics  of  unsteady,  separated  flows  and  (2)  controlling  unsteady, 
separated  flows  to  enhance  aerodynamic  performance.  Several  types  of  flow  fields  were 
studied  to  gain  an  understanding  of  unsteady,  separated  flows.  These  included  flow  over 
rectangular  wings,  flow  over  delta  wing,  flow  over  forebodies  and  flow  in  the  vicinity  of 
the  tail  on  a  fighter  type  aircraft.  The  studies  were  both  computational  and  experimental. 
Control  efforts  focused  on  the  pulsed  blowing  to  maintain  attached  flow,  continuous 
blowing  to  control  vortex  breakdown  on  delta  wing,  and  the  use  of  neural  networks  to 
predict  and  control  unsteady  aerodynamic  forces.  Results  for  these  topics  are 
summarized  in  the  next  eight  sections. 

For  a  more  detailed  account  of  the  research,  see  the  attached  Work  Unit  publications  and 
presentations  listing 


Unsteady  Vortex  Dynamics  and  Surface 
Pressure  Topologies  on  a  Finite  Pitching  Wing 

S.  J.  Schreckand  H.  E.  Helin 
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Unsteady  Vortex  Dynamics  and  Surface  Pressure  Topologies 

on  a  Finite  Pitching  Wing 


Scott  J.  Schreck*  and  Hank  E.  Helint 
U.S.  Air  Force  Academy,  Colorado  Springs,  Colorado  80840 


A  straight  wing  having  an  NACA  0015  cross  section  and  rectangular  planform  was  attached  to  a  circular 
splitter  plate.  This  configuration  was  pitched  at  a  constant  rate  to  angles  exceeding  the  static  stall  angle.  The 
unsteady,  vortex-dominated  flow  that  developed  over  the  wing  and  splitter  plate  was  characterized  in  detail 
using  surface  pressure  measurements  and  flow  visualization.  Both  types  of  data  showed  that  the  leading-edge 
vortex  underwent  profound  three-dimensional  alterations  to  cross  section  and  convection  over  the  entire  wing 
span.  These  changes  in  leading-edge  vortex  structure  and  kinematics  were  correlated  with  prominent  spanwise 
variations  in  force  coefficients.  When  appropriately  dissected,  visualization  results  and  pressure  data  suggested 
physical  mechanisms  to  account  for  these  three-dimensional  variations  in  unsteady  forces  and  surface  pressures. 


Nomenclature 

C„  =  normal  force  coefficient 
c  =  wing  chord  length,  cm 

Cp  =  pressure  coefficient 

LE  -  leading  edge 

Re^  =  chord  Reynolds  number 
s  —  wing  span  length,  cm 

t  =  time,  s 

“  nondimensional  time,  tUJc 
=  test  section  velocity,  m/s 
X  =  chordwise  distance  from  leading  edge 
y  =  spanwise  distance  from  splitter  plate 
a  =  instantaneous  angle  of  attack,  deg 
a  —  pitch  rate,  rad/s 

=  nondimensional  pitch  rate,  cafU,, 

Introduction 

HREE-DIMENSIONAL  dynamically  separated  flows 
continue  to  be  intensely  studied.  Energetic  large-scale 
vortical  structures  are  generated  and  transiently  reside  on 
wings  dynamically  pitched  through  the  static  stall  angle  of 
attack.  These  vortical  structures  radically  alter  the  pressure 
distributions  on  a  wing,  thus  producing  greatly  amplified  aero¬ 
dynamic  forces  and  moments.  If  thoroughly  understood  and 
properly  controlled,  three-dimensional  dynamically  separated 
flows  have  the  potential  to  confer  dramatic  performance  en¬ 
hancements  upon  future  aircraft. 

Flow  visualization  studies  have  documented  the  morphol¬ 
ogy  of  portions  of  three-dimensional  unsteady  flows  elicited 
by  pitching  wings. These  studies  concentrated  attention  on 
the  prominent  leading-edge  and  wingtip  vortices,  and  suc- 
cessfijlly  constructed  simplistic  physical  models  based  upon 
vorticity  conservation  to  explain  vortex  anchoring  and  taper¬ 
ing  near  the  wingtip.  Fewer  investigations  have  employed 
surface  pressure  measurements  to  characterize  three-dimen¬ 
sional  unsteady  flowfield  development  on  pitching  wings.'^  ® 
These  studies  found  unsteady  spanwise  pressure  distributions 
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and  is  not  subject  to  copyright  protection  in  the  United  States. 
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Laboratory,  2354  Vandenberg  Dr.,  Ste.  6H79.  Member  AIAA. 

tAssistant  Professor,  Department  of  Aeronautics,  2410  Faculty 
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to  be  augmented  near  the  wingtip  in  comparison  to  the  steady 
spanwise  distribution. 

Unsteady  flows  near  wing-wall  junctures  have  been  shown 
to  be  as  complex  as  those  near  the  wingtip.  Shih'"  examined 
the  dynamically  separated  flowfield  over  a  wing  spanning  a 
test  section  in  which  the  flow  speed  varied  sinusoidally.  Ve¬ 
locity  surveys  revealed  stropg  spanwise  three-dimensionality 
that  exhibited  convective  behavior  suggesting  the  presence  of 
separated  vortical  eddies.  Horner  and  colleagues^*’  employed 
flow  visualization  to  show  that  the  unsteady  flow  near  the 
juncture  between  a  fixed  wall  and  a  pitching  wing  was  two 
dimensional  immediately  following  leading-edge  vortex  ini¬ 
tiation.  However,  as  the  vortex  convected  along  the  wing 
chord,  the  vortex  arched  away  from  the  wing  surface  near 
center  span,  and  symmetrical  counter-rotating  cells  formed 
on  the  wing  surface.  Three-dimensional  disruptions  to  the 
leading-edge  vortex  were  explained  with  a  model  that  relied 
upon  mutual  induction  between  orthogonal  vortex  segments. 

Schreck  and  coworkers and  Klinge  et  al.^-  investigated 
the  dynamically  separated  flow  near  the  juncture  between  a 
wing  and  splitter  plate  undergoing  constant  rate  and  sinusoidal 
pitching,  respectively.  Both  experiments  measured  unsteady 
surface  pressures,  enabling  unambiguous  characterization  of 
leading-edge  vortex  convection  and  quantification  of  time- 
dependent  normal  forces.  In  the  juncture  region,  the  leading- 
edge  vortex  convected  at  a  faster  velocity^and  unsteady  stall 
occurred  earlier.  Preliminary  vorticity  dynamics  models  were 
postulated  to  explain  alterations  to  leading-edge  vortex  con¬ 
vection  near  the  juncture. 

The  current  investigation  combines  quantitative  surface 
pressure  topologies  with  flow  visualization  at  identical  non- 
dimensional  pitch  rates  and  closely  matched  Reynolds  num¬ 
bers.  This  approach  allows  the  presence  and  behavior  of  vis¬ 
ualized  vortex  structures  to  be  unambiguously  confirmed  by 
distinctive  pressure  signatures.  This  methodology  also  enables 
vorticity  generation  sites  and  rates  to  be  inferred  using  surface 
pressure  topologies.  Thus,  vorticity  not  aggregated  into  vor¬ 
tical  structures  can  also  be  accounted  for.  Both  flow  visual¬ 
ization  and  surface  pressure  measurements  show  significant 
spanwise  variability  in  the  unsteady  flowfield  development. 
To  account  for  this  variability,  a  model  is  hypothesized  for 
the  three-dimensional  unsteady  flowfield  on  the  wing-splitter 
plate  configuration. 

Experimental  Methods 
Surface  Pressure  Measurement 

Surface  pressure  measurements  were  performed  in  the  Frank 
J.  Seiler  0.91-  x  0.91-m  low-speed  wind  tunnel  located  at  the 
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U.S.  Air  Force  Academy.  A  rectangular  planform  wing  with 
15.24-cm  chord  length  was  fabricated  from  hollow  aluminum 
NACA  0015  airfoil  stock.  The' basic  wing  was  29.10  cm  long 
and  was  equipped  with  a  fitting  on  the  outboard  end  that 
permitted  arbitrary  length  extensions  to  be  added.  Fifteen 
miniature  pressure  transducers  were  installed  inside  the  hol¬ 
low  basic  wing  model.  These  transducers  were  close  coupled 
to  the  wing  surface  through  pressure  ports  located  along  the 
chord  line,  3.05  cm  inboard  of  the  basic  wing  end.  Pressure 
transducer  signals  were  low-pass  filtered  (300  Hz  cutoff)  and 
amplified  by  a  gain  of  500.  The  resulting  signals  were  then 
sampled  and  digitized  by  the  data  acquisition  system. 

A  circular  aluminum  splitter  plate,  30,48  cm  in  diam  and 
0  64  cm  thick,  was  machined  to  a  sharp  edge  around  the  plate 
perimeter.  The  splitter  plate  had  an  NACA  0015  cutout  cen- 
tered  in  it,  which  allowed  it  to  slide  onto  the  wing  and  be 
positioned  at  arbitrary  span  locations.  To  effectively  move 
the  pressure  ports  along  the  wing  span,  the  splitter  plate  was 
first  positioned  at  the  desired  distance  from  the  pressure  ports. 
Then,  a  tip  extension  of  the  correct  length  was  added  to  the 
basic  wing,  bringing  the  span  length  to  30.48  cm  and  main¬ 
taining  aspect  ratio  constant  at  2.0.  The  chordwise  row  of 
pressure  ports  was  successively  moved  to  11  span  locations, 
effectively  distributing  pressure  ports  over  the  wing  surface 
as  shown  in  Fig.  1.  All  11  tip  extensions  used  in  these  ex¬ 
periments  terminated  in  a  square  tip. 

In  Fig.  1,  spanwise  pressure  port  locations  range  from  0.0 
to  0.80  span  outboard  of  the  splitter  plate.  Chordvyise  pressure 
port  locations  range  from  ”  0.90  to  0.90  chord,  with  0.0  chord 
corresponding  to  the  wing  leading  edge.  Positive  chord  values 
denote  the  wing  upper  surface,  whereas  negative  chord  sig¬ 
nifies  the  lower  surface.  Unsteady  surface  pressures  measured 
at  these  port  locations  were  contour-plotted  using  a  linear 
interpolation  between  adjacent  grid  points  in  both  the  chord- 
wise  and  spanwise  directions. 

Model  pitching  was  driven  by  a  3.5-hp  synchronous  stepper 
motor.  The  wing-splitter  plate  combination  was  mounted  on 
a  steel  shaft  2.86  cm  in  diam  that  was  connected  to  the  stepper 
motor  through  a  gear  linkage  having  a  4:1  reduction  ratio. 
The  independent  variables  explored  in  the  surface  pressure 
experiments  included  nondimensional  pitch  rates  0.05,  0.10, 
and  0.20,  as  well  as  spanwise  pressure  port  locations  0.0,  0.05, 
0.10,  0.15,  0.25,  0.375,  0.50,  0.625,  0.70,  0.75,  and  0.80  span 
outboard  of  the  splitter  plate.  Wing  pitch  axis  was  located  at 
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Fig.  1  Effective  distribution  of  pressure  port  locations  over  the  wing 
planform. 


0.25  chord  for  the  experimental  range  of  nondirnensional 
rate,  and  at  0.33  chord  for  nondimensional  pitch  rate  0.10. 
For  each  parameter  combination,  20  consecutive  wing  pitch 
motions  were  sampled  and  ensemble-averaged.  Test  section 
velocity  was  held  constant  at  9.14  m/s,  corresponding  to  a 
chord  Reynolds  number  of  6.9  x  10-*. 

Flow  Visualization 

Flow  visualization  was  performed  in  the  U.S.  Air  Force 
Academy  0.51-  x  0.38-m  water  tunnel.  Dimensions  of  the 
flow  visualization  wing-splitter  plate  model  were  identical  to 
those  of  the  surface  pressure  measurement  model.  A  0.13- 
mm-wide  slot  was  cut  through  the  leading  edge  of  the  flow 
visualization  wing  along  the  entire  wing  span.  This  slot  al¬ 
lowed  dye  to  flow  from  the  hollow  wing  interior  to  the  wing 
exterior.  Both  wing  and  splitter  plate  were  painted  white  to 
maximize  visibility  of  the  dark  blue  dye.  Dye  was  injected 
directly  into  the  boundary  layer,  visualizing  the  vorticity  in 

the  unsteady  flowfield.  . 

The  wing-splitter  plate  model  used  for  flow  visualization 
was  mounted  on  a  steel  shaft  0.95  cm  in  diam.  The  model 
was  pitched  at  constant  rate  by  a  24-V  dc  gear  motor  that  was 
connected  to  the  model  through  a  gear  linkage  having  a  5:1 
reduction  ratio.  The  visualized  flowfield^was  illuminated  by 
two  100-W  incandescent  lamps.  Flow  visualization  images  were 
recorded  from  the  wing  planform  perspective  at  30  frames/s 
by  a  VHS  video  camera. 

Water  tunnel  flow  visualization  was  performed  at  a  test 
section  velocity  of  0.44  m/s  and  a  water  temperature  of  12.8°C. 
This  yielded  a  chord  Reynolds  number  of  5.6  x  10*  that 
corresponded  closely  to  that  of  6.9  x  10^  for  wind-tunnel 
surface  pressure  measurements.  Flow  was  visualized  for  non- 
dimensional  pitch  rates  of  0.05,  0.10,  apd  0.20,  using  pitch 
axis  locations  of  0.25  and  0.33  chord.  Volume  coefficient  of 
injection*-'  out  of  the  leading-edge  slot  was  modest,  being  0.01. 

Results 

Pitching  the  wing-splitter  plate  configuration  beyond  static 
stall  elicited  dramatic  three-dimensional  modifications  to 
leading-edge  vortex  structure.  Visualized  leading-edge  vortex 
locaticTns  and  structures  correlated  closely  with  minima  in  the 
surface  pressure  topologies.  Three-dimensional  modifications 
to  the  leading-edge  vortex  and  correlation  with  surface  pres¬ 
sure  topologies  persisted  throughout  the  experimental  range. 
Leading-edge  vortex  kinematics  at  all  span  locations  were  also 
well-correlated  with  spanwise  normal  force  loading. 

Angle-of-attack  angle  histories  for  the  three  measured  non- 
dimensional  pitch  rates  of  0.05,  0.10,  and  0.20  are  shown  in 
Fig.  2.  All  three  histories  begin  at  0.0  deg  and  end  at  60.0 
deg.  The  beginning  of  these  histories  coincides  with  the  in¬ 
ception  of  surface  pressure  data  acquisition.  Thus,  the  plots 
in  Fig.  2  can  be  used  in  conjunction  with  subsequent  plots  to 
convert  nondimensional  time  to  instantaneous  wing  angle  of 
attack.  For  surface  pressure  measurements,  the  wing  was 
pitched  at  a  constant  rate  from  0.0  to  60,0  deg  to  collect  upper 
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surface  data,  and  from  0  deg  to  —60.0  deg  to  acquire  lower 
surface  data.  It  should  be  noted  that  none  of  the  angle-of- 
attack  histories  exhibits  prominent  nonlinearity  or  suffers  from 
appreciable  acceleration  or  deceleration  transients. 

Vortex  Structure  and  Surface  Pressure  Topologies 

Figure  3  shows  a  representative  series  of  five  surface  pres¬ 
sure  coefficient  topologies  and  corresponding  dye  flow  visu¬ 
alization  photographs.  Nondimensional  pitch  rate  was  0.10, 
and  the  wing  was  pitched  about  0.33  chord.  Note  that  surface 
pressure  topologies  extend  0.80  span  from  the  splitter  plate, 
whereas  flow  visualization  shows  the  entire  wing  span.  For 
clarity,  the  wing  planform  has  been  outlined  in  the  visuali¬ 
zation  photographs.  Lower  surface  pressure  topologies  were 
characterized  by  mild  two-dimensional  pressure  gradients  that 
underwent  limited  temporal  evolution.  Thus,  lower  surface 
data  will  not  be  presented  here. 

The  frame  1  =  5.28)  surface  pressure  topology  in  Fig. 

3  is  dominated  by  an  arc-shaped  suction  ridge  that  begins  near 
the  intersection  of  the  splitter  plate  and  wing  leading  edge. 
The  suction  ridge  proceeds  outboard  across  the  wing  span, 
curves  back  to  0.33  chord  at  center  span,  and  again  ap¬ 
proaches  the  leading-edge  0.80  span  outboard  of  the  splitter 
plate.  Suction  ridge  position  corresponds  closely  to  the  up¬ 
stream  boundary  of  the  visualized  leading-edge  vortex.  Near 
the  splitter  plate  and  wingtip,  at  0.18  chord,  prominent  lo¬ 
calized  suction  peaks  disrupt  suction  ridge  uniformity.  Both 
suction  peaks  lie  close  to  portions  of  the  vortex  that  exhibit 
pronounced  flexure  in  the  visualization.  Upstream  of  the  suc¬ 
tion  ridge,  two  wedge-shaped  contours,  0.25  and  0.50  span 
from  the  splitter  plate  and  just  behind  the  leading  edge,  de¬ 
note  pressure  maxima  that  lie  on  either  side  of  a  low-pressure 
region.  This  low-pressure  region  corresponds  to  a  chordwise 
dye  streak  that  connects  the  wing  leading  edge  and  the  vortex 
upstream  boundary. 

In  frame  2  =  5.88),  the  suction  ridge  remains  evident 

in  the  surface  pressure  topology.  The  inboard  end  of  the 
suction  ridge  is  still  located  near  the  intersection  of  the  splitter 
plate  and  wing  leading  edge,  and  the  outboard  end  still  ter¬ 
minates  near  the  wing  leading  edge.  Visualization  is  consis¬ 
tent,  showing  the  inboard  and  outboard  ends  of  the  vortex 
still  located  near  the  leading  edge.  Suction  ridge  magnitude 
has  been  substantially  diminished  near  center  span,  where 
visualization  indicates  that  the  leading-edge  vortex  has  arched 
up  and  away  from  the  wing  surface.  Two  local  suction  peaks 
are  still  present  in  the  surface  pressure  topology.  The  inboard 
peak  has  moved  downstream  and  is  located  at  0.36  chord, 
while  the  outboard  suction  peak  remains  at  0.18  chord.  Suc¬ 
tion  peak  locations  lie  close  to  portions  of  the  visualized  vortex 
that  have  undergone  deformation,  but  remain  near  the  wing 
surface. 

The  frame  3  =  6-30)  surface  pressure  topology  indicates 

the  outboard  end  of  the  suction  ridge  still  terminates  near  the 
wing  leading  edge.  However,  the  inboard  end  of  the  ridge 
has  moved  downstream  and  is  now  located  near  midchord. 
Visualization  shows  the  vortex  remains  at  the  leading  edge 
near  the  wingtip,  but  has  begun  to  convect  downstream  near 
the  splitteT  plate.  Two  local  suction  peaks  are  still  present 
near  portions  of  the  vortex  that  appear  to  have  undergone 
deformation.  The  inboard  peak  has  moved  downstream  and 
is  located  at  0.44  chord,  while  the  outboard  suction  peak 
remains  at  0.18  chord.  In  the  ridge  central  portion,  suction 
magnitude  remains  depressed,  and  the  vortex  arch  that  ap¬ 
peared  in  the  previous  flow  visualization  frame  is  still  clearly 
visible.  In  addition,  rotating  cells  have  formed  at  the  juncture 
where  the  legs  of  the  vortex  arch  join  the  remainder  of  the 
vortex.  At  each  location,  cell  rotation  corresponds  to  that  of 
the  vortex,  with  each  cell  rotating  in  a  sense  opposite  the 
other. 

In  frame  4  (r„d  =  6.96),  the  surface  pressure  topology  ex¬ 
hibits  a  suction  ridge  that  begins  near  the  outboard  portion 
of  the  leading  edge  and  has  maximum  magnitude  there.  Vis¬ 


ualization  is  consistent,  showing  the  vortex  still  anchored  to 
the  leading  edge  near  the  wingtip.  The  suction  ridge  then 
extends  aft  and  inboard,  decreasing  in  magnitude,  and  finally 
terminating  near  the  intersection  of  the  splitter  plate  and  trail¬ 
ing  edge.  The  outboard  suction  peak  is  still  visible  at  0.37 
chord,  whereas  the  inboard  suction  peak  no  longer  exists. 
Visualization  is  again  consistent,  showing  the  outboard  ro¬ 
tational  cell  near  midchord  and  the  inboard  one  centered  over 
the  trailing  edge.  These  two  cells  are  connected  by  diffuse 
dye  streaks  resembling  the  previously  well-defined  arch.  Fi¬ 
nally,  a  concentration  of  dye  is  visible  immediately  adjacent 
to  the  splitter  plate.  This  structure  remains  connected  to  the 
inboard  cell  by  a  well-defined  dye  streak. 

The  frame  5  (/nd  =  7.26)  surface  pressure  topology  consists 
principally  of  a  broad,  straight  suction  ridge  extending  from 
the  leading  edge  near  the  wingtip  to  the  trailing  edge  near 
the  splitter  plate.  The  outboard  suction  peak  still  persists  near 
the  wingtip  and  a  new  suction  peak  has  formed  adjacent  to 
the  splitter  plate.  This  new  suction  peak  has  substantially 
lower  magnitude  than  suction  peaks  observed  earlier  in  the 
pitch  motion,  either  inboard  near  the  splitter  plate  or  out¬ 
board  near  the  wingtip.  This  is  once  again  consistent  with  the 
visualization,  which  shows  the  outboard  end  of  the  vortex  lies 
close  to  the  wing  leading  edge.  From  here  the  vortex  proceeds 
inboard  and  aft  until  it  reaches  the  remaining  rotational  cell. 
From  the  rotational  cell,  the  vortex  extends  upward  and  aft 
before  flexing  toward  the  splitter  plate  and  extending  inboard. 
After  contacting  the  splitter  plate,  the  vortex  proceeds  down 
toward  the  wing  surface  and  forward  toward  the  leading  edge. 
It  then  terminates  on  the  wing  surface  adjacent  to  the  splitter 
plate  near  the  new  suction  peak. 

Figure  4  contains  three  flow  visualization  photographs 
showing  the  leading-edge  vortex  at  chose njntermediate  times 
in  the  pitch  motion  for  the  experimental  range  of  nondimen¬ 
sional  pitch  rates.  These  photographs  were  selected  to  capture 
the  initial  disruptions  to  the  leading-edge  vortex  gradual  arc 
shape  that  were  precursors  to  vortex  arching.  In  the  visual¬ 
izations,  disruptions  appeared  as  discontinuities  in  the  ap¬ 
parent  upstream  boundary  of  the  leading-edge  vortex.  Dis¬ 
ruptions  subsequently  developed  into  flexure  points  for  the 
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Fig.  4  Flow  visualizations  of  leading-edge  vortex  at  initial  disruption 
for  experimental  range  of  and  pitch  axis  at  0.33c. 
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vortex  arch,  and  then  rotational  cells.  In  each  photograph, 
spanwise  location  of  both  the  inboard  and  outboard  disrup¬ 
tions  are  marked  by  arrows.  Spatial  and  temporal  relation¬ 
ships  between  these  initial  disruptions  and  the  wing  are  re¬ 
corded  in  Table  1. 

Table  1  indicates  that  higher  nondimensional  pitch  rates 
drive  the  inboard  initial  disruption  closer  to  the  splitter  plate. 
In  addition,  higher  nondimensional  pitch  rates  force  the  out¬ 
board  initial  disruption  farther  from  the  splitter  plate  and 
nearer  the  tip.  When  the  initial  disruptions  occur,  both  in¬ 
board  and  outboard  disruptions  appear  at  the  same  chord 
location-.  However,  increasing  nondimensional  pitch  rate  moves 
initial  occurrence  of  the  disruption  farther  forward  on  the 
wing  chord,  and  delays  it  to  higher  angles  of  attack. 

Figure  5  contains  three  surface  pressure  contour  plots,  iden¬ 
tical  in  format  to  those  presented  in  Fig.  3,  for  nondimensional 
pitch  rates  of  0.05,  0.10,  and  0.20.  These  three  contour  plots 
correspond  in  orientation  and  scale  to  the  three  flow  visual¬ 
ization  photographs  shown  in  Fig.  4.  Thus,  these  topologies 
document  surface  pressure  at  the  time  when  the  leading-edge 
vortex  underwent  initial  disruption.  All  three  contour  plots 
exhibit  a  prominent  arc-shaped  suction  ridge  that  intersects 
the  wing  leading-edge  inboard  near  the  splitter  plate  and  out¬ 
board  near  the  wingtip.  Near  center  span,  the  suction  ridge 
curves  downstream  to  approximately  midchord.  Prominent 
suction  peaks  are  evident  on  the  inboard  portion  of  the  suction 
ridge,  located  at  0.10,  0.10,  and  0.05  span  for  nondimensional 
pitch  rates  of  0.05,  0.10,  and  0.20.  On  the  outboard  part  of 
the  suction  ridge,  a  suction  peak  is  visible  at  0.75  span  for  all 
three  nondimensional  pitch  rates.  Higher  nondimensional  pitch 


Table  1  Spatial  and  temporal  occurrence  of  initial  disruptions 
shown  by  visualization  in  Fig.  4 
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Fig.  5  Surface  pressure  topologies  at  initial  vortex  disruption  for 
experimental  range  of  and  pitch  axis  at  0,25c. 


rates  also  show  denser  contour  distributions,  indicating  steeper 
surface  pressure  gradients  in  both  the  chordwise  and  spanwise 
directions  at  higher  nondimensional  pitch  rates. 

Leading-Edge  Surface  Pressure 

Figure  6  shows  a  representative  contour  plot  documenting 
leading-edge  surface  pressure  coefficient  variation  with  non- 
dimensional  time  and  span  location.  This  plot  corresponds  to 
a  nondimensional  pitch  rate  of  0.10  and  pitch  axis  location  of 
0.33  chord.  Beginning  at  nondimensional  time  0.0,  predom¬ 
inantly  two-dimensional  suction  increase  is  indicated  by  the 
horizontal  contours  in  the  lower  portion  of  the  plot.  However, 
gradual  upward  curvature  at  both  the  inboard  and  outboard 
ends  of  the  contours  indicates  a  lag  in  suction  increase  near 
the  wing  root  and  tip.  After  nondimensional  time  4.00,  lead¬ 
ing-edge  suction  increase  becomes  radically  nonuniform  along 
the  wing  span,  and  three  regions  are  discernible  in  the  plot. 
In  the  central  region,  between  0.15-0.50  span,  leading-edge 
suction  first  attains  a  maximum  at  a  nondimensional  time  of 
4.14.  Suction  subsequently  declines  at  a  nearly  constant  rate 
as  shown  by  the  uniform  contour  spacing  in  the  vertical  di¬ 
rection. 

Inboard,  between  0.0-0.15  span,  the  suction  peak  occurs 
0.89  nondimensional  time  units  later  than  it  did  in  the  central 
region.  Suction  reaches  a  local  maximum  at  0.05  span  and 
nondimensional  time  4.80,  as  indicated  by  the  concentric  con¬ 
tours.  Following  attainment  of  this  maximum,  suction  decline 
is  slow  at  first,  and  then  accelerates,  as  shown  by  the  closer 
contours  after  nondimensional  time  6.0.  Terminal  rate  of  suc¬ 
tion  decrease  is  similar  for  the  central  and  inboard  regions  of 
the  plot,  as  evidenced  by  the  comparable  distances  between 
adjacent  contour  lines.  The  outboard  region  is  dramatically 
different  from  both  the  central  and  inboard  regions.  Here, 
leading-edge  suction  peaks  at  0.80  span  and  nondimensional 
time  6.00,  nearly  1.50  nondimensional  time  units  after  oc¬ 
currence  of  the  maxima  in  either  the  central  or  inboard  region . 
After  reaching  this  outboard  peak,  suction  decreases  at  an 
appreciably  slower  rate  than  it  did  in  either  the  central  or 
inboard  region  of  the  plot. 

Figure  7  is  a  summary  plot  that  records  the  nondimensional 
time  of  leading-edge  suction  maximum  at  11  span  stations, 
for  the  experimental  range  of  nondimensional  pitch  rate  and 
pitch  axis  location.  Data  comprising  this  graph  were  extracted 
from  contour  plots  similar  to  Fig.  6.  Figure  7  shows  that 
leading-edge  suction  reaches  a  maximum  at  earlier  nondi- 


Fig.  6  Representative  contour  plot  of  leading-edge  vs  and  span. 
a*  =  0.10  and  pitch  axis  is  located  at  0.33c. 
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Fig.  7  Summary  plot  documenting  ^nd  of  leading-edge  suction  collapse 
vs  span  for  experimental  range  of 


Fig,  8  Contour  plot  C„  vs  and  span,  a"*"  -  0.10  and  pitch  axis  is 
located  at  0.33c. 

mensional  times  for  higher  nondimensional  pitch  rates,  cor¬ 
responding  to  higher  instantaneous  angles  of  attack.  Moving 
the  pitch  axis  from  0.25  to  0.33  chord  for  nondimensional 
pitch  rate  0.10  also  imposes  a  slight  delay  upon  leading-edge 
suction  peak  occurrence. 

All  four  plots  in  Fig.  7  are  dominated  by  a  broad,  level 
central  region.  This  portion  of  the  curve  shows  that  leading- 
edge  suction  maxima  occurred  first,  and  nearly  simultane¬ 
ously,  over  the  central  portion  of  the  wing  span  for  each 
nondimensional  pitch  rate.  Inboard  and  outboard  of  these 
level  regions,  all  four  curves  bend  sharply  upward.  This  in¬ 
dicates  significant  delays  in  the  attainment  of  leading-edge 
suction  maxima  near  both  splitter  plate  and  wingtip.  All  four 
plots  show  suction  maxima  occurring  first  in  the  broad  central 
region  of  the  wing  span,  followed  by  maxima  inboard  near 
the  splitter  plate  and,  finally,  outboard  near  the  wingtip.  Higher 
nondimensional  pitch  rates  generally  drove  the  spanwise  lo¬ 
cations  of  suction  maximum  delay  closer  to  the  splitter  plate 
and  wingtip. 

Spanwise  Normal  Force  Loading 

Figure  8  is  a  representative  contour  plot  showing  normal 
force  coefficient  variation  with  nondimensional  time  and  span 
location.  Normal  force  initially  increases  uniformly  across  the 
wing  span  with  nondimensional  time  as  shown  by  the  hori¬ 
zontal  contours.  However,  at  nondimensional  time  2.0,  the 
contours  at  the  outboard  span  stations  begin  to  curve  upward 


and  are  spaced  farther  apart,  indicating  a  reduced  rate  of 
normal  force  increase  there.  At  nondimensional  time  5.0, 
similar  upward  curvature  and  reduced  contour  density  occurs 
at  the  inboard  span  stations,  also  indicating  reduced  rates  of 
"normal  force  increase  near  the  splitter  plate.  Prior  to  non- 
dimensional  time  5.0,  normal  force  variation  remained  highly 
uniform  across  the  span,  in  spite  of  the  minor  perturbations 
noted  above. 

After  nondimensional  time  5.0,  spanwise  variation  in  nor¬ 
mal  force  suffered  severe  three-dimensional  disruption.  Fig¬ 
ure  8  exhibits  a  pronounced  inboard  normal  force  maximum 
of  2.95,  as  evidenced  by  the  concentric  contours  centered  at 
0.1  span  and  nondimensional  time  5,94.  A  similar  normal 
force  maximum  of  larger  spatial  and  temporal  extent,  and 
magnitude  3.50,  occurs  outboard  as  indicated  by  the  concen¬ 
tric  contours  centered  at  0.5  span  and  nondimensional  time 
7.62.  Between  these  two  regions  on  the  plot,  from  approxi¬ 
mately  0.2  to  0.4  span,  lies  an  area  characterized  by  sparse 
contours  of  significantly  lower  magnitude.  In  this  central  re¬ 
gion,  normal  force  stalls  at  nondimensional  time  5.58  and 
attains  a  maximum  value  of  2.16.  Normal  force  then  declines 
more  slowly  than  it  does  either  inboard  near  the  splitter  plate 
or  outboard  in  the  tip  region. 

At  nondimensional  time  8.0,  normal  force  coefficient  de¬ 
creased  to  a  minimum  of  approximately  1.25  at  0.0  span,  and 
simultaneously  increased  to  a  maximum  of  3.75  at  0.5  span. 
Between  these  two  span  locations,  at  0.25  span,  a  region  of 
nearly  vertical  contours  indicates  negligible  temporal  change 
in  normal  force.  Subsequently,  normal  force  increased  in¬ 
board,  decreased  outboard,  and  underwent  little  change  in 
the  vicinity  of  0.25  span.  At  nondimensional  time  1 1 .0,  normal 
force  became  uniform  across  the  span  at  a  magnitude  of  ap¬ 
proximately  2.0. 

Figure  9  is  a  summary  plot  containing  information  extracted 
from  contour  plots  similar  to  Fig.  8.  Figure  9  records  the 
nondimensional  time  of  normal  force  maximum,  or  normal 
force  stall,  at  11  span  stations  for  the  experimental  range  of 
nondimensional  pitch  rate  and  pitch  axis  location.  Figure  9 
shows  that  at  any  given  span  station  normal  force  stall  occurs 
at  earlier  nondimensional  times  for  higher  nondimensional 
pitch  rates,  corresponding  to  higher  instantaneous  angles  of 
attack.  For  nondimensional  pitch  rate  0.10,  moving  the  pitch 
axis  from  0.25  to  0.33  chord  slightly  delayed  normal  force 
stall. 

All  four  curves  in  Fig.  9  show  similar  spanwise  stall  pro¬ 
gressions.  Normal  force  stall  occurred  first  at  0.375  span  for 
nondimensional  pitch  rate  0.20,  and  at  0.25  span  for  the  re¬ 
maining  three  conditions.  Stall  occurred  next  inboard,  near 
the  splitter  plate.  Here,  inboard  of  0.25  span,  stall  occurred 
last  at  0.10  span  for  nondimensional  pitch  rate  0.05,  and  at 
0.05  span  for  the  other  three  conditions.  Finally,  normal  force 
stalled  outboard,  near  the  wingtip.  The  spanwise  location  of 
final  stall  moved  inboard  with  higher  nondimensional  pitch 


Fig.  9  Summary  plot  showing  of  C„  stall  vs  span  for  experimental 
range  of 
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Fig.  10  Summary  plot  of  stall  C„  vs  span  for  experimental  range 
of  a’^. 

rate,  progressing  from  0.70  to  0.625  span,  and  finally  to  0.50 
span  for  nondimensional  pitch  rates  of  0.05,  0.10,  and  0.20. 

Figure  10  summarizes  data  taken  from  contour  plots  like 
Fig.  8,  and  documents  stall  normal  force  coefficient  at  11  span 
stations  for  the  experimental  range  of  nondimensional  pitch 
rate  and  pitch  axis  location.  Figure  10  shows  that,  at  any  given 
span  station,  stall  normal  force  coefficient  increases  with  higher 
nondimensional  pitch  rate.  For  nondimensional  pitch  rate  0.10, 
moving  the  pitch  axis  from  0.25  to  0.33  chord  generally  in* 
creases  stall  normal  force  slightly. 

Again,  all  four  curves  in  Fig.  10  display  similar  trends.  Stall 
normal  force  was  lowest  at  0.25  span  for  all  four  conditions. 
For  each  of  the  four  curves,  absolute  maxima  of  stall  normal 
force  generally  occurred  outboard  near  the  wingtip.  Such  ab¬ 
solute  maxima  occurred  at  0.70  and  0.50  span  for  nondimen¬ 
sional  pitch  rates  of  0.20  and  0.10.  However,  for  nondimen¬ 
sional  pitch  rate  0.05,  stall  normal  force  reached  absolute 
maximum  at  0.10  span.  In  addition  to  the  absolute  maxima, 
local  maxima  also  occurred  at  the  opposite  end  of  the  span. 
Local  maxima  occurred  at  0.10  span  for  nondimensional  pitch 
rates  0.05  and  0.10,  and  at  0.05  span  for  nondimensional  pitch 
rate  0.20. 

Discussion 

Flow  visualization,  surface  pressure  topologies,  and  span- 
wise  normal  force  loading  histories  provided  diverse  infor¬ 
mation  regarding  the  unsteady  flow  over  the  pitching  wing- 
splitter  plate  configuration.  Correlating  these  data  enabled 
formulation  of  a  cohesive  account  of  vortex  kinematics, 
underlying  vorticity  dynamics  and  spanwise  normal  force 
loading. 

Temporal  and  spatial  fluctuations  of  considerable  magni¬ 
tude  were  observed  in  flow  visualization,  surface  pressure 
topologies,  and  spanwise  normal  force  loading  histories.  In 
contrast.  Fig.  2  showed  that  the  rigid  wing-splitter  plate  com¬ 
bination  pitched  at  constant  rate  in  uninterrupted  fashion. 
Thus,  kinematic  or  geometric  anomalies  did  not  contribute 
to  these  prominent  temporal  and  spatial  variations. 

Dye  flow  visualization  revealed  a  large,  energetic  leading- 
edge  vortex  that  persisted  during  a  significant  interval  of  the 
pitch  motion  and  extended  over  a  large  area  of  the  wing 
planform.  Following  initiation  near  the  wing  leading  edge, 
the  vortex  convected  downstream  toward  the  wing  trailing 
edge  while  growing  in  cross  section.  Vortex  growth  and  con¬ 
vection  were  accompanied  by  radical  three-dimensional  de¬ 
formation.  Visualized  vortex  structures  and  locations  were 
corroborated  by  prominent  pressure  minima  that  dominated 
surface  pressure  topologies,  consistent  with  observations  made 
by  Walker  and  colleagues'^  for  surface  pressure  distributions 
on  two-dimensional  airfoils.  Modifications  to  vortex  structure 
or  location  were  accompanied  by  corresponding  alterations 
to  surface  pressure  topologies. 

Soon  after  a  predominantly  two-dimensional  initiation.  Fig. 
3  (frame  1,  =  5.28)  showed  that  the  leading-edge  vortex 


had  experienced  limited  convection  near  the  splitter  plate  and 
wingtip.  Here,  pressure  minima  occurred  in  the  surface  pres¬ 
sure  contours  consistent  with  enhanced  vortex-surface  inter¬ 
action  due  to  vorticity  straining.  The  central  portion  of  the 
vortex  had  convected  downstream,  forming  a  moderate  vortex 
arc  lying  near  the  wing  surface.  This  limited  three  dimen¬ 
sionality  was  reflected  in  the  Fig.  8  normal  force  contour  plot 
as  a  mild  disruption  to  the  previously  two-dimensional  contour 
lines. 

The  gradual  vortex  arc  was  initially  disrupted  as  visualized 
in  Fig.  4,  and  confirmed  by  surface  pressure  contours  in  Fig. 

5.  Initial  vortex  disruptions  moved  farther  apart,  and  closer 
to  both  splitter  plate  and  wingtip,  with  higher  nondimensional 
pitch  rate.  The  two  delayed  leading-edge  suction  maxima  of 
Fig.  7  also  moved  apart,  shifting  the  inboard  one  closer  to 
the  splitter  plate,  and  the  outboard  one  closer  to  the  wingtip 
as  nondimensional  pitch  rate  increased.  Comparing  Table  .1 
to  Fig.  7  showed  that  leading-edge  suction  collapse  preceded 
initial  vortex  disruption  by  approximately  1. 0-2.0  nondimen¬ 
sional  time  units,  and  that  this  intervening  period  shortened 
with  higher  nondimensional  pitch  rate.  Table  1  also  showed 
that  the  vortex  was  located  closer  to  the  wing  leading  edge 
at  initial  disruption  with  increasing  nondjmensional  pitch  rate. 

These  correlations  indicate  that  leading-edge  suction  col¬ 
lapse  indicated  the  demise  of  pressure  gradients  responsible 
for  vorticity  production.  Vorticity  production  curtailment,  in 
turn,  disrupted  the  vorticity  sheet  feeding  the  downstream 
leading-edge  vortex.  Acharya  and  Metwally'^  have  shown  that 
disruption  of  vorticity  production  in  the  leading-edge  region 
imposes  pronounced  modifications  upon  leading-edge  vortex 
kinematics  downstream. 

That  portion  of  the  vortex  between  .the  two  initial  vortex 
disruptions  subsequently  arched  over  the  wing  surface  as  vis¬ 
ualized  in  Fig.  3  (frame  2,  =  5.88).  Vortex  arching  kin¬ 

ematics  in  the  current  investigation  were  consistent  with  those 
visualized  by  Schreck  et  al.,“  Freymuth,''  and  Horner  et  al.“’ 
Vortex  arching  was  accompanied  by  suction  collapse  on  the 
wing  area  vacated  by  the  vortex.  Here,  vortex-surface  sepa¬ 
ration  distance  increased,  reducing  vortex-surface  interaction 
and  attenuating  suction  as  demonstrated  by  Panaras.*^ 
Following  vortex  arching,  the  apex  of  the  arch  convected 
downstream  at  a  significantly  higher  speed  than  the  remainder 
of  the  vortex.  Vortex  convection  velocity  was  locally  accel¬ 
erated  when  the  arched  portion  of  the  vortex  encountered 
stronger  freestream  influence. However,  vortex  arch  height 
above  the  wing  prevented  vortex  convection  from  strongly 
impacting  either  surface  pressures  or  normal  forces.'*'  This 
was  indicated  by  the  sparse  and  then  vertical  contours  im¬ 
mediately  inboard  of  center  span  in  Fig.  8.  Vorticity  straining 
in  the  remainder  of  the  vortex  was  augmented  by  the  vortex 
arch  and  continued  to  amplify  vortex-surface  interaction,  both 
near  the  splitter  plate  and  the  wingtip.  This,  in  turn,  prompted 
corresponding  maxima  in  surface  pressure  topologies  and 
spanwise  normal  force  loading  topologies. 

At  the  two  sites  where  the  vortex  flexed  to  accommodate 
arching,  counter-rotating  cells  appeared  at  the  wing  surface, 
as  visualized  in  Fig.  3  (frame  3,  —  6.30).  Initially,  these 

cells  were  symmetric  in  the  spanwise  direction  about  a  line 
just  inboard  of  the  wing  center  span.  Similar  cell  pairs  have 
been  observed  under  dynamic  conditions  by  Shih,*^  Schreck,'"^ 
and  Horner  et  al.,'“  and  in  the  static  regime  by  Winkelmann 
and  Barlow.  Rotational  sense  of  the  cells  was  consistent  with 
that  of  the  contiguous  leading-edge  vortex  that  had  arched 
up  over  the  wing  surface.  At  this  time,  the  leading-edge  vortex 
near  the  splitter  plate  began  to  convect  away  from  the  leading 
edge,  locally  deforming  the  vortex.  Figure  8  showed  approx¬ 
imately  equal  values  of  normal  force  inboard  and  outboard. 
However,  normal  force  had  peaked  and  was  decreasing  in¬ 
board  near  the  splitter  plate  while  continuing  to  increase  out¬ 
board  near  the  wingtip. 

After  the  symmetric  pair  of  counter-rotating  cells  appeared, 
these  cells  and  the  associated  vortex  convected  asymmetri- 
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cally.  Figure  3  (frame  4,  =  6.96)  showed  that  the  inboard 

cell  reached  the  trailing  edge  well  in  advance  of  the  outboard 
one,  even  though  both  began  at  the  same  chord  location.  This 
observation  is  consistent  with  observations  made  by  Schreck 
et  al."  who  observed  accelerated  vortex  convection  near  the 
wing'SpIitter  plate  juncture.  Inboard,  vortex  convection  began 
shortly  after  leading-edge  suction  finally  collapsed  near  the 
splitter  plate,  as  recorded  in  Fig,  7.  The  influence  of  asym¬ 
metric  convection  was  consistent  with  normal  force  stall  times 
shown  in  Fig.  9,  which  showed  that  normal  force  stall  occurred 
inboard  significantly  earlier  than  it  did  at  outboard  span  sta¬ 
tions.  Asymmetric  vortex  convection  was  prompted  by  a  com¬ 
bination  of  influences.  Outboard,  near  the  wingtip,  vortex 
convection  was  retarded  by  pinning  at  the  tip  region  as  vis¬ 
ualized  by  Freymuth,**  and  quantitatively  corroborated  by 
Robinson  and  coworkers.^  Inboard,  vortex  convection  was 
accelerated  by  mutual  induction  acting  between  the  inboard 
leg  of  the  vortex  arch  and  the  vortex  image  presented  by  the 
splitter  plate. 

Finally,  Fig.  3  (frame  5,  =  7.26)  showed  the  outboard 

portion  of  the  leading-edge  vortex  began  to  convect  down¬ 
stream  and  away  from  the  leading  edge  near  the  wingtip.  Soon 
after  this  time  Fig.  8  showed  normal  force  reaching  a  maxi¬ 
mum  near  the  wingtip  and  beginning  to  decline.  Experiments 
by  Lorber  et  al.’*^  have  recorded  similar  lift  maxima  restricted 
to  the  wingtip  region.  Inboard,  visualization  showed  the  vor¬ 
tex  above  and  behind  the  wing  surface,  but  still  apparently 
connected  to  the  wing  surface  by  two  vortex  segments  ema¬ 
nating  from  the  wing  surface.  By  this  time,  the  inboard  portion 
of  the  wing  had  stalled  and  nearly  reached  minimum  normal 
force. 

Figure  7  showed  that  leading-edge  suction  collapse,  near 
both  the  splitter  plate  and  wingtip,  was  significantly  delayed 
relative  to  central  span  locations.  Outboard,  down  wash  pro¬ 
duced  by  the  prominent  wingtip  vortex  delayed  leading-edge 
suction  collapse  compared  to  center  span.  Similar  delays  in 
leading-edge  suction  collapse  inboard  suggest  a  concentration 
of  streamwise  vorticity,  of  opposite  sense  to  that  at  the  wing¬ 
tip,  is  also  present  near  the  wing-splitter  plate  juncture. 

Surface  pressure  gradients  on  the  wing-splitter  plate  con¬ 
figuration  gave  rise  to  localized  concentrations  of  streamwise 
vorticity.  Low  pressures  produced  on  the  wing  during  pitching 
would  induce  secondary  flows  over  the  splitter  plate  toward 
the  wing.  Thus,  streamwise  vorticity  would  be  generated  and 
convected  over  the  splitter  plate  toward  the  wing,  collecting 
in  the  wing-splitter  plate  juncture.  Although  no  streamwise 
vortex  was  observed  in  the  juncture  region,  a  concentration 
of  streamwise  vorticity  not  coalesced  into  a  vortex  would  also 
induce  down  wash  and  preserve  leading-edge  suction.  Juncture 
vortices  have  been  visualized  by  Klinge  et  al.*-  near  the  junc¬ 
ture  of  a  pitching  wing-splitter  plate  configuration. 

The  amount  of  streamwise  vorticity  present  at  either  the 
wingtip  or  the  juncture  dictated  the  influence  that  it  exercised. 
In  Fig.  7,  leading-edge  suction  collapse  delays  were  longer 
and  extended  farther  inboard  from  the  wingtip  than  delays 
near  the  splitter  plate.  This  asymmetry  suggests  that  larger 
amouTits  of  streamwise  vorticity  were  present  near  the  wingtip 
than  at  the  juncture  at  leading-edge  suction  collapse.  Figure 
7  also  showed  that  leading-edge  suction  collapse  delay  was 
displaced  closer  to  the  wingtip  and  juncture  with  increasing 
nondimensional  pitch  rate.  This  is  consistent  with  the  model 
proposed  by  Luttges  and  Kennedy,'*’  wherein  higher  nondi¬ 
mensional  pitch  rates  reduce  total  vorticity  production  and 
prompt  more  condensed  vorticity  distributions. 

Conclusions 

A  generic  wing-splitter  plate  configuration  was  pitched  at 
constant  rate  to  angles  exceeding  the  static  stall  angle.  The 
resulting  three-dimensional  unsteady  flowfield  was  charac¬ 
terized  using  dye  flow  visualization  and  unsteady  surface  pres¬ 
sure  measurements.  Visualized  three-dimensional  deforma¬ 
tions  of  the  leading-edge  vortex  were  corroborated  by  complex 


surface  pressure  topologies  that  were  similarly  three-dimen¬ 
sional  in  nature.  Spanwise  variations  in  normal  force  loading 
were  consistent  with  three-dimensional  vortex  deformations. 
Mechanisms  based  upon  vorticity  dynamics  were  postulated 
to  account  for  these  observed  kinematics. 

The  initially  two-dimensional  unsteady  vortex  rapidly 
underwent  three-dimensional  deformation  along  the  entire 
wing  span.  Three-dimensional  deformation  began  when  vor¬ 
tex  convection  was  arrested  near  the  splitter  plate  and  wingtip, 
but  continued  near  center  span.  Soon,  disruptions  appeared 
in  the  vortex  upstream  boundary,  and  were  quickly  followed 
by  vortex  arching  above  the  wing  surface  near  center  span. 
Subsequently,  vortex  convection  near  the  splitter  plate  out¬ 
paced  that  near  the  tip,  with  the  inboard  portion  of  the  vortex 
being  shed  prior  to  that  outboard. 

Vortex  kinematics  were  driven  by  an  identifiable  set  of 
vorticity  dynamics.  Streamwise  vorticity  accumulated  near  the 
splitter  plate  and  wingtip.  These  vorticity  accumulations  tem¬ 
porarily  arrested  vortex  convection  and  delayed  vortex  arch¬ 
ing.  Mutual  induction  acting  between  the  vortex  arch  and  the 
image  of  it  presented  by  the  splitter  plate  accelerated  vortex 
convection  near  the  splitter  plate. 

Spanwise  normal  force  loading  waj  consistent  with  vortex 
kinematics.  Two-dimensional  and  moderately  three-dimen¬ 
sional  vortex  structure  yielded  normal  force  loading  that  was 
uniform  along  the  span.  Prominently  three-dimensional  vor¬ 
tex  structure  prompted  correspondingly  nonuniform  spanwise 
normal  force  loading.  Normal  force  stall  times  varied  sub¬ 
stantially  across  the  span,  and  were  closely  correlated  to  vor¬ 
tex  arching  or  shedding. 

This  experiment  has  examined  three-dimensional  dynami¬ 
cally  separated  flow  over  a  generic  wjng-splitter  plate  config¬ 
uration,  Many  aspects  of  this  flow  have  been  successfully 
characterized  and  explained,  though  only  for  a  limited  pa¬ 
rameter  range.  As  such,  the  current  investigation  provides  a 
basis  for  further  investigations  employing  configurations  of 
more  practical  interest  and  encompassing  broader  parameter 
ranges. 
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Abstract 

The  three-dimensional  Navier-Stokes  equations  were 
solved  for  the  unsteady  flow  about  a  generic  wing-wall 
configuration  for  constant  rate  pitching  from  0°  to  60® 
nngU  of  attack  at  non-dimensional  pitdi  rates  of  0.2  and 
0.1.  The  configuration  consisted  of  a  straight  wing  with  a 
NACA  0015  aoss  section  bounded  on  one  end  by  a  flat 
flitter  plate  orthc^onally  attadied  to  the  wing  root  and  by 
a  free  tip  on  the  other.  The  numerical  results  compare  fa¬ 
vorably  with  experimental  flow  visualization  and  surface 
pressure  measurements,  particularly  at  the  higher  pitch 
rate.  The  strong  interaction  of  the  wing-tip  vortex  with 
tiie  dynamic  stall  vortex  and  the  resulting  anqilification  of 
the  normal  force  near  the  wing  tip  are  predicted  in  good 
agreemmit  with  the  eiqieriment.  fu  contrast,  the  viscous 
intArarrifin  at  the  wing-wall  juncture  is  much  less  signifi¬ 
cant  in  altering  the  overall  force  distribution  on  flie  wing. 
Spanwise  and  chordwise  vorticity  contours  and  instantane¬ 
ous  limiting  streamlines  are  combined  with  surface  pres¬ 
sure  and  normal  force  plots  to  provide  a  detailed  under¬ 
standing  of  the  temporal  evolution  of  the  massively  sqia- 
rated,  three-dimensional,  unsteady  flow  structure. 

Introduction 

When  a  wing  or  other  lifting  surface  is  rapidly  pitched 
beyond  its  static  stall  angle  of  attadc,  the  motion  is  accom¬ 
panied  by  intense  vorticity  generation  at  the  wing  surface. 
Due  to  a  disparity  between  the  rate  of  production  and  the 
rate  of  transport,  the  vorticity  accumulates  on  the  leeward 
side  of  die  wing  near  the  leading-edge  and  the  boundary 
layer  remains  attached  well  beyond  the  point  at  which 
separation  occurs  under  static  conditions.  Eventually,  the 
boundary  layer  separates  in  response  to  an  increasingly 
adverse  pressure  gradient  or  the  upstream  propagation  of  a 
trailing-edge  separation  and  the  concentrated  vorndty 
rolls  up  to  fonn  a  discrete  vortex,  the  dynamic  stall  vortex, 
whidi  detaches  and  convects  downstream  over  the 
upper  surface  of  the  wing.  The  accumulated  vorticity  in 
the  attached  boundary  layer  and,  later,  the  dynamic  stall 
vOTtex  alters  the  circulation  and  the  resulting  pressure 
distphiirirm  on  the  wing,  producing  transient  forces  and 
mftmpntR  that  are  much  larger  than  those  realized  under 
static  condirinns.  The  entire  process,  known  as  dynamic 
tfali  is  important  in  many  aerodynamic  applications.  The 
dynamir,  M  ovetshoot  that  accompanies  the  formation  of 
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the  dynamic  stall  vortex  potentially  can  be  exploited  to 
greatly  pnbnnri*  the  maneuvering  performance  of  combat 
aircraft.  Dynamic  stall  is  also  a  limiting  factor  in  the 
performance  of  helicopter  rotors  and  wind  turbines  due  to 
vibration  and  fatigue  problems  induced  by  the  large  tran¬ 
sient  aerodynamic  loads.  The  development  of  control 
methodologies  that  effectively  exploit  dynamic  stall  to  en¬ 
hance  aerodynamic  pafonnance  as  well  as  the  modifica¬ 
tion  or  avoidance  of  adverse  unsteady  effects  require  a 
thorough  understanding  of  the  physics  of  unsteady  sepa¬ 
rated  flows. 

Reffiynrft  1  provides  a  recent  overview  and  assess¬ 
ment  of  the  issues  related  to  dynamic  stall  and  unsteady 
separated  flows.  As  a  result  of  extorsive  experimental  and 
computational  study,  an  understanding  of  many  aspects  of 
two-dimensional  dynamic  stall  is  relatively  mature.  A 
CTini'iar  iwiHarctanHing  of  tbree-dimensloual  dynamic  stall 
is  at  a  much  earlier  stage  of  development.  A  fimte  wing, 
for  example,  admits  a  variety  of  unsteady  flow  interactions 
not  present  in  tire  two-dimensitmal  case  including  wing-tip 
effects,  wall  effects,  and  spanwise  flow.  These  fliree-di- 
niangifmal  interactions  significantly  complicate  the  two- 
rfiniOTicinnal  dynamic  stall  model  described  earlier. 

Robinson  and  Wissler^  experimentally  studied  the 
large-ampUtude  pitching  motion  for  a  straight  rectangular 
wing  and  found  that  the  interaction  of  the  wing-tip  and 
dynamic  stall  vortices  significantly  amplified  the  sectional 
lift  mrffiriftnts  near  the  wing  tip.  Schreck  et  al.^,  Klinge 
et  al.^,  and  Schreck  and  Helin®  investigated  the  unsteady 
separated  flow  for  a  finite  wing  undergdng  ^usoidal  and 
ffinctant  rate  pitching-  Both  surface  pressure  measure¬ 
ments  and  flow  visualization  showed  significant  ^anwise 
variation  with  respect  to  leading-edge  vortex  aoss  section, 
convection  rate,  and  normal  force  coefficient,  correspond¬ 
ing  to  prominent  wing-tip  and  wing-juncture  interactions. 
Lc^t^  examined  the  wing-tip  interaction  fa  a  pitching 
finite  wing  at  higher  Mach  and  Rejmolds  numbers  and 
noted  similar  trends  in  wing-tip  force  amplification  and 
dynamic  stall  propagation  patterns. 

Numaotts  two-dimensional  Navia-Stokes  solutions 
have  been  presented  fa  an  airfoil  undergoing  either  oscil¬ 
latory  a  constant  rate  pitdiing  motions.  Representative 
woks  include  those  of  Visbal'^’*®,  Rumsey  and 
Anderson^^  Patterson  and  Laber^^,  Ghia  et  al.^^,  and 
TCnight  and  Qioudhuri^^.  Three-dimensional,  Navier- 
Stokes  solutions  fa  the  unsteady  flow  goierated  by  the 
forced  motion  of  a  finite  wing  include  the  woik  of 
Guruswamy^^  fa  a  swept  wing  in  pitching  motion. 


Chadeijian  and  Giiruswamy^^  and  Yang  and 
Przekwas^^*^^  for  the  sinusoidal  oscillation  of  rectangular 
and  forward-swept  wings,  Visbal^^  and  Visbal  and 
Gordnier^^  foe  a  pitching  delta  wing,  and  Gordnier^^  for  a 
rolling  delta  wing. 

In  Ref.  22,  Newsome  numerically  studied  the  three- 
dimensional,  unsteady,  separated  flow  in  the  juncture  of  a 
wing-wall  configuration  corresponding  to  the  experimental 
model  of  Schreck  et  al.^*^.  The  numerical  results  provided 
a  detailed  view  of  the  influence  of  the  wall  in  inducing  a 
complex  spanwise  secondary  flow  on  the  wing  in  the 
juncture  region.  Because  interest  was  focused  on  the  wing 
juncture,  the  wing-tip  interaction  was  not  modeled,  thus 
limiting  a  comparison  with  the  experiment  to  the  juncture 
region.  However,  the  experimental  results  of  Schreck^’^ 
and  others^*^  clearly  show  that  the  wing-tip  flow  is  the 
predominant  three-dimensional  interaction.  The  present 
calculations  model  both  the  wing  tip  and  wing  juncture  but 
focus  on  wing  tip  and  the  interaction  of  the  wing-tip  vor¬ 
tex  and  the  dynamic  stall  vortex  leading  to  locally  ampli¬ 
fied  sectional  force  coefficients  in  the  tip  region.  Two 
different  pitch  rates  are  considered  foe  large  amplitude 
pitching  to  60^  angle  of  attack,  corresponding  to  maneu¬ 
vers  representative  of  those  likely  to  be  employed  by  future 
combat  aircraft 


Governing  Equations 


The  governing  equations  are  the  three-dimensional, 
compressible,  time-dependent,  Navier-Stokes  equations, 
written  in  generalized  coordinates  and  conservation  form 
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and  p  is  doisity,  (u,v,w)  are  Cartesian  velocities  in  the 
inertial  coordinate  frame,  F  is  the  total  energy  per  unit 
volume,  p  is  pressure 


P  =  (Y-1)[f-^p(u2+v2+w2)]  (3) 

and  contravarient  velocities  (Uy,W)  are 

V  =  +  TlyV  +  +  Tl,  (4) 


The  generalized  coordinates  (^,11,0  correspond  to  the  di¬ 
rections  along  the  wing  ^an  (^),  the  wing  chord  circum¬ 
ferential  direction  (ti),  and  the  direction  normal  to  the 
wing  surface  (0-  In  the  thin-layer  assumption,  only  de¬ 
rivatives  normal  to  non-slip  wall  surfaces  are  retained, 
greatly  simplifying  the  viscous  terms  in  Eqn.  1. 
Consistent  with  the  results  of  Ref.  22  in  which  the  effects 
of  the  viscous  ctoss  derivative  terms  were  found  to  be 
small  and  then  confined  to  the  juncture  boundary  layer 
region,  only  the  thin  layCT  terms  in  the  two  viscous  direc¬ 
tions  (§,  0  are  retained.  Thus  3  G,/a  Ti  ■  0  and  the  vis¬ 
cous  term,  ,  may  be  written  as 
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A  corresponding  e}q}ressiQn  for  is  obtained  by  substi¬ 
tuting  ^  and  W  in  place  of  ^  and  U  in  Eqn.  5.  In  the 
above,  the  equations  are  non-dimensionalized  in  terms  of 
fise  stream  density,  p„.  and  sound  speed,  c„.  M„  and 
are  the  fieestream  Mach  number  and  Reynolds  munber 
and  F^  is  the  Prandtl  number.  Equation  1  is  closed  by 
Stoke's  hypothesis  for  bulk  viscosity  and  Sutherland's  law 
for  molecular  viscosity  (p).  All  calculations  are  for  lami¬ 
nar  flow. 


Computational  Method 

NASA.  Langley’s  computer  code.  CFL3D^’^,  was 
utilized  as  the  flow  solver.  The  code  is  an  upwind-biased, 
finite-volume,  approximate  factorization  algorithm  em¬ 
ploying  either  Van  Leer's  flux  vector  splitting  or  Roe's  flux 
difference  splitting.  Details  of  the  algorithm  are  given  in 
Ref.  23.  hi  Ref.  22,  several  modifications  woe  made  to 
improve  the  accuracy  and  efSdency  of  the  code  for  un¬ 
steady  applications.  The  diagonalized  algorithm  was 
generalized  for  unsteady  flows  with  moving  grids  by  in- 
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eluding  the  time-dependeat  transfonnatiem  derivatives. 
^  j.T]  ,  in  the  defi^tion  of  the  eigenvalues  correspond¬ 
ing  to  the  diagonalized  operators.  To  ensure  temporal, 
conservadon  and  achieve  second  order  tune  accuracy, 
(<I’-K)>  three  point  backward  tiine  differencing  and  a 
subiteradon  capability  were  added  to  die  baseline  algo¬ 
rithm.  The  resulting  algorithm,  used  in  the  present  calcu- 
ladons,  solves  direedy  for  the  primidve  variables 
q  -  (p.  u,  V,  w,  pf  as: 
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posed.  The  far  field  boundary  condition  specification  was 
based  upon  Riemann  invariants  for  one-dimensional  flow 
nonnd  to  the  far  field  boundary.  For  the  O-grid  wing-tip 
boundary  and  the  C-grid  wake  cut,  continuation  boundary 
conditions  were  applied  implicitly  by  reordering  the  ap¬ 
proximately  factored  matrix  equations  (Eqn.  6a  and  Eqn. 
6c)  so  that  these  boundaries  are  treated  as  regular  interior 
points.  All  other  boundary  conditions  were  applied  ex¬ 
plicitly.  Unsteady  pitching  flow  solutions  were  started 
from  a  nominally  steady  solution  at  zero  degrees  angle  of 
attack.  Calculations  were  performed  on  the  Cray  C90 
computer  at  the  U.S.  Army  Corps  of  Engineers  Waterways 
Experiment  Station  at  Vicksburg,  Mississippi. 
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In  the  above,  M  =  dQjdq .  7  is  the  Jacobian  of  the  coordi¬ 
nate  transfarmation,  R  is  the  steady  state  residual,  A.  and 
T  are  the  eigmivalues  and  eigenvectors  of  the  inviscid  fhix 
Jacobians,  and  5>v  is  a  temporary  variable  defined  above 
by  context  To  obtain  second  order  temporal  accuracy,  an 
inner  iteration,  m,  at  die  n+1  time  level  is  employed  such 
that  Ecpis.  6a,  6b,  6c  are  solved  iteratively  until  the  right 
hand  side  of  Eqit  6a.  (conesponding  to  die  residual  of  the 
unsteady  equation)  is  less  than  a  prescribed  tolerance, 
usually  two  to  three  orders  of  magnitude  below  the  steady 
state  residual,  R. 

Roe's  ffaix  difference  splitting  was  used  to  evaluate  the 
inviscid  terms.  A  third-order  upwind-biased  interpolation 
of  the  state  variables,  q,  was  ukd  to  obtain  the  required 
flux  values  at  cell  intmfaces.  Second-order  central  differ¬ 
encing  was  used  to  evaluate  die  viscous  terms. 


Experiment 

Schreck^'^  experimentally  smdied  the  flow  for  a 
straight  wing  having  a  NACA  0015  cross-section  and  a 
rectangular  planfoim.  The  wing  chord  and  ^an  were  6 
and  12  inches  respectively.  A  circular  splitter  plate  12 
inches  in  diameter  was  ordioganaHy  attached  to  die  root  of 
the  wing.  The  wing  tip  was  a  square  end  with  no  round¬ 
ing.  The  model  was  pitched  from  0  to  60*^  about  the 
and  chord  locations  for  five  non-dimensional  pitch 
rates.  =  QC/(Af„  c«)  =  0.01. 0.02, 0.05, 0.10, 0.20.  In 

the  above,  Q  is  the  pitch  rate  in  rad/sec.  C  is  the  wing 
chord  and  and  are  the  freestteam  Mach  number 
and  speed  of  sound.  The  configuration  was  studied  in  the 
wind  tunnel  and  water  tunnel  faciUties  of  the  F.  J.  Seiler 
Research  Laboratory  and  the  USAF  Academy,  combining 
surface  pressure  measurements  and  flow  visualization  to 
obtain  a  more  complete  understanding  of  the  three-di¬ 
mensional,  unsteady  vortex  dynamics.  Uie  wind  tunnel 
model  was  instrumented  to  obtain  surface  pressure  at  fif¬ 
teen  chordwise  locations  on  both  the  iqiper  and  lower  wing 
surfaces  at  11  spanwise  locations  (0.0,  0.05,  0.1,  0.15, 
0.25,  0.375,  0.5,  0.625,  0.7,  0.75,  0.8).  All  pressure 
leadings  were  ensemble  averaged  over  20  different  pitch 
cycles.  The  test  section  velocity,  17*  =  30  ft/sec,  produced 
a  Reyndds  number  based  upon  chord  of  about  69.000  and 
a  fireestream  Mach  number,  » 0.026.  In  addition  to 
surface  pressure  measurmnents,  water  tunnel  flow  visuali¬ 
zation  was  accomplished  by  injecting  dye  into  the  bound¬ 
ary  layer  at  the  leading  edge  dong  the  span  of  the  wing. 

corresponding  chord  Reynolds  numbor  was  56,000. 


The  grid  was  rigidly  attached  to  the  wing-wall  con¬ 
figuration  and  rotated  with  the  configuration  without  de¬ 
formation.  The  Navier-Stdees  equations  were  solved  with 
respect  to  inertial  coordinates  at  the  coordinate  locations 
determined  by  the  moving  grid.  On  the  wing  and  wall 
surfaces,  no  slip  boundary  conditions  were  applied  for  the 
moving  surfaces.  An  adiabatic  wall  was  assumed  and  an 
expressirm  for  the  nonnal  pressure  gradient  in  terms  of  the 
w^  surface  acceleration  provided  surface  pressure.  On 
the  wing  root  boundary  beyond  the  wall  surface,  symmetry 
conditions  conesponding  to  an  invisdd  slip  wall  were  im- 


During  most  of  die  pitching  motion,  flow  visualization 
revealed  a  large,  energetic,  leading-edge  vortex  that  ex¬ 
tended  over  a  broad  area  of  the  wing  planfoim.  Both  sur¬ 
face  pressure  and  flow  visualization  data  show  that  the 
leadi^-edge  vortex  underwent  significant  three-dimen¬ 
sional  alteration  in  convection  rate  and  cross-section. 
These  modifications  were  accompanied  by  prominent 
spanwise  variations  in  stall  angle  and  normal  force  coeffi¬ 
cient.  Three  distinct  spanwise  regions  were  identified 
corresponding  to  the  wing-tip,  wing-juncture,  and  center 
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span  regions.  The  wing-tip  vortex  successfully  maintained 
leading-edge  suction  well  beyond  regions  further  inboard. 
The  leading-edge  vortex  remained  anchored  to  the  wing 
tip  until  late  in  the  pitching  cycle  when  it  detached  and 
convected  downstream  to  produce  a  broad  pressure  mini¬ 
mum  on  the  outboard  wing  planform.  Overall  normal 
force  stall  was  associated  with  the  departure  of  the  lead¬ 
ing-edge  vortex  from  the  wing  tip.  The  flow  near  the 
wing-wall  juncture  was  initially  more  diffuse  due  to  the 
presoice  of  the  wall.  Later,  however,  the  pitching  motion 
induced  steep  transverse  pressure  gradients  in  die  juncture 
regimi.  The  most  prominent  feature  associated  with  the 
center  span  region  was  the  tendoicy  of  the  leading-edge 
vortex,  for  certain  test  conditions,  to  arch  up  and  form  a 
horseshoe-shaped  vortex  stretching  from  the  "pinned'' 
wing-tip  region  to  the  splitter-plate  wall.  Two  counter¬ 
rotating  cells  were  observed  on  the  wing  surface  near  the 
base  of  the  legs  of  the  arched  vortex.  All  of  the  effects  de¬ 
scribed  above  were  more  pronounced  at  higher  non-di¬ 
mensional  pitch  rates. 

Numerical  Results 

Three-dimensional  Navier-Stokes  solutions  were  com¬ 
puted  for  the  wing-wall  configuration  pitching  about  the 
quarter  chord  for  non-dimensional  pitch  rate  constants  of 
Q*  =  0.2  and  0.1.  An  analytical  expression  fw  the  non- 
dimensional  pitch  rate,  similar  to  that  of  Visbal’,  was  de¬ 
rived  as  a  best  fit  to  the  experimental  pitch  rate  schedule 


and  computational  results.  Etqterimentally,  no  informa¬ 
tion  was  available  concerning  the  turbulent  state  of  the 
flow.  As  a  practical  matter,  appropriate  turbulence  models 
for  massively  separated,  three-dimensional,  unsteady  flows 
do  not  currently  exist.  Accordingly,  the  flow  was  assumed 
to  remain  laminar  and  calculations  w^  made  at  the  lower 
experimental  Reynolds  number  of  56,000  used  in  the  water 
tunnel  flow  visualization  studies. 

Fig.  2  shows  the  baseline  grid  for  the  wing-wall 
configuration.  The  grid  trilogy  is  a  C  grid  along  the 
chord  and  a  1/2  O  grid  along  the  wing  span.  Two  modifi¬ 
cations  were  made  to  the  configuration  to  simplify  the  grid 
generation  process.  First,  the  wing  tip  was  rounded  with  a 
circular  tip  cap  corresponding  to  the  local  chord  thickness. 
Second,  the  aft  end  of  the  circular  splitter  plate  was 
changed  to  a  rectangular  shape  to  conform  to  the  C-grid 
topology.  The  grid  attended  15  chords  outward  and  con¬ 
sisted  of  213  circumferential  points,  75  radial  points,  and 
61  spanwise  points.  It  was  constructed  to  be  nearly  or¬ 
thogonal  to  the  wing  and  splitter  plate  surfaces  and  to 
conform  to  the  semi-circular  outer  boundary  of  the  splitter 
plate.  The  minimum  cell-centered  spacing  adjacent  to  the 
wing  and  wall  surfaces  was  set  at  As/C  «=  0.00005,  corre¬ 
sponding  to  roughly  17-18  points  in  the  steady  state 
boundary  layer.  "The  wing  surface  was  defined  by  152  cir¬ 
cumferential  grid  points. 

Numerical  Accuracy 


fl-e 


V 


(7) 


where  the  nondimensional  time,  f’  =  tc^lC,  is  valid  over 
the  range  0^  C  In  the  above,  is  the  non-di¬ 
mensional  end  time,  and  t*  and  ^  are  constants  chosen  to 
match  the  experimental  pitch  rate  time  history. 


33.0 

12.6133 

18.7648 

53.75 

2.94423 

3.45023 

Table  1  Constants  for  Pitch  Rate  Schedule,  Eqn.  7 

Figure  1  provides  a  comparison  of  die  resulting  pitch 
angle  schedule  (obtained  by  integrating  EqrL  7)  with  the 
experimmital  schedule. 

Experimentally,  botii  the  wind  and  water  turmel  flow 
regimes  correspond  to  incompressible  flow. 
Computationally,  a  higher  Mach  number  (M„  =  0.2)  was 
spedfied  in  order  not  to  degrade  the  convergence  of  the 
compressible  flow  solver.  Even  at  this  low  Mach  number, 
comiH^ssibility  effects  are  significant  foe  the  pitch  rates 
considered  and  do  impact  the  comparison  of  experimental 


Because  of  the  large  resource  requirements  for  three- 
dimensional  unsteady  solutions,  preliminary  temporal  and 
spatial  accuracy  assessments  were  accomplished  in  two 
dimensienrs  for  the  equivalent  pitching  akfoil.  Due  to  the 
numerous  approximations  made  in  the  implicit  t^iators, 
the  algorithm  represented  by  Eqn.  6  is  not  unconditionally 
stable.  For  the  three-dimensional  giid  shown  in  Fig.  2,  the 
global  unsteady  time  step  size  was  limited  by  the  local 
stability  limit  of  the  grid  volumes  at  the  trailing  edge  of 
the  wing  tip.  However,  tins  stability-limited  time  step  was 
still  less  restrictive  than  that  required  fa:  temporal  accu¬ 
racy,  even  for  two-dimensional  flows.  Thus  a  mote  accu¬ 
rate  sdieme  employing  larger  but  more  expensive  time 
st^s  may  be  more  efficient  than  a  less  accurate  scheme 
requiring  more  but  less  expensive  time  steps.  Figures  3a 
and  3b  demonstrate  this  point  for  a  two-dimensional  airfoil 
pitching  about  the  quarter  chord  at  a  non-dimensional 
pitch  rate  of  =  0.2  according  to  tiie  pitch  rate  schedule 
of  Eqn  7.  The  grid  was  identical  to  the  three-dimensional 
grid  at  the  wing  root.  Figure  3a  compares  solutions  em¬ 
ploying  three  point  backward  time  differencing  and  subit¬ 
eration  to  achieve  second  order  temporal  accuracy.  The 
least  expensive  option.  At'*' =  0.002,  employing  two 
subiteratiois.  corre^y  resolved  the  temporal  dynamics  of 
the  flow.  Figure  3b  shows  the  results  df  first-order  time- 
accurate  solutions  without  subiteration  for  the  same  time 
step  sizes  as  Fig.  3a  and  for  time  step  sizes  four  times 
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smaller,  coiresponding  to  an  equivalent  computational  cost 
of  a  single  larger  Hme  stqp  with  three  subiterations.  There 
is  a  si^ificant  amount  rf  scatter  in  the  results  with  the 
smaller  time  steps  consistently  approaching  but  not 
equaling  the  more  accurate  solutions  with  subiteration.  In 
fhig  particular  case,  the  second-order  solution  with  two 
subiterations  was  more  accurate  than  the  first-order  solu¬ 
tion  without  subiteradon  but  using  a  time  step  eight  times 
smaller.  The  first-order  solution  was  obtained  at  a  compu¬ 
tational  cost  2.66  (8/3)  times  higher.  Furthermore,  since 
second-order  accuracy  was  obtained  with  the  diagonal  al¬ 
gorithm  which  is  2.6  times  less  expensive  than  the  block 
inversion  algorithm,  the  total  improvement  in  code  effi¬ 
ciency,  compared  to  the  first-order,  time-accurate,  block 
inversion  algorithm,  was  nearly  a  factor  of  seven.  In  all 
the  results  presented  in  Fig.  3,  the  C-grid  wake  boundary 
condition  was  updated  implicitly.  As  noted  in  Ref.  22,  the 
rapid  convergence  of  the  suibiteration  force  coefficient 
histories  was  not  achieved  unless  this  boundary  condition 
was  updated  implidtly. 

The  effect  of  grid  resolution  on  the  two-dimensional 
solution  is  shown  in  Fig.  4  in  which  the  baseline  213x75 
grid  was  uniformly  refined  to  produce  a  425x149  grid.  A 
comparison  of  normal  and  tangential  force  coefficients 
indicates  that  the  baseline  two-dimensional  grid  accurately 
follows  the  overall  fine  grid  prediction  but  tends  to  smooth 
the  smaller  fluctuations  present  in  the  fine  grid  calcula¬ 
tion.  Although  finer  grid  resolution  is  desirable,  when 
extrapolated  to  three-dimensional  unsteady  flow,  this  level 
of  grid  refinement  greatly  increases  computer  memory  and 
processing  requirements  for  relatively  modest  accuracy 
improvements.  The  three-dimensional  grid  was  con¬ 
structed  to  provide  equal  viscous  resolution  on  the  wing 
and  wall  surfaces  and  to  provide  fine  spanwise  resolution 
of  the  wing-tip  region. 

Unsteady  Flow  *  Initial  Condition 

The  pitching  wing  computations  were  initiated  from  a 
nominally  converged  steady  state  solution  at  0^  angle  of 
attack.  In  fact,  the  flow  over  a  NACA  0015  airfoil  at  the 
Reynolds  numbers  considered  (56,000-69,(XX))  does  not 
exhibit  a  steady  state  solution,  but  rather  is  characterized 
by  periodic  vortex  shedding  and  a  lift  coefficient  that  osdl- 
l^s  around  a  zero  mean  value.  The  computation  devel¬ 
oped  as  an  attached  flow  solution,  followed  by  trailing 
edge  boundary  layer  separaticm,  prr^agation  of  the  sepa¬ 
rate  regirm  upstream,  and  eventually,  asymmetric  separa¬ 
tion  and  periodic  vortex  shedding.  In  the  experiment, 
differences  in  starting  conditions  were  minimized  by  en¬ 
semble  averaging  data  over  20  pitch  cycles.  In  the  compu¬ 
tation,  a  phase-neutral  starting  solution  was  obtained  by 
terminating  the  calculation  after  the  forward  propagation 
of  the  separated  region  but  prior  to  the  develqpment  of 
significant  flow  asymmetry.  Figure  5  shows  die  limiting 
streamlines  on  the  upper  surface  of  the  wing  and  the 


splitter  plate  at  QP  angle  of  attack.  The  primary  separation 
Hnft  is  evident  at  a  chord  location  c£  xlC  ^  0.55.  Near  the 
juncture,  the  separation  line  was  swept  forward  to  a  chord 
location  xlC  •  0.29.  The  distinctive  recirculation  pattern 
observed  downstream  of  the  separation  line  near  the  junc¬ 
ture  was  produced  by  a  forward  unstable  node,  a  trailing 
stable  focus  and  two  nearby  saddle  points.  Additionally, 
tiie  flow  separated  on  the  splitter  plate  at  a  position  0.1 
chord  upstream  of  the  wing  leading  edge,  producing  a 
"hOTseshoe  vortex*'  that  wrapped  around  the  wing  and  was 
then  convected  downstream,  introducing  longitudinal  vor- 
tidty  into  the  juncture  region.  At  the  wing  tip,  separation 
was  suppressed  by  the  spanwise  flow  around  the  tip  and 
the  separation  line  was  displaced  almost  to  the  wing  trail¬ 
ing  e(tee. 

Unsteady  Flow  -  Pitch  Rate  of  =  0.2 

For  the  pitching  motion  of  a  three-dimensional  wing, 
the  influence  of  the  wing-tip  vortex  on  the  dynanuc  stall 
vortex  is  the  most  interesting  and  most  significant  interac¬ 
tion  in  terms  of  its  effect  on  the  forces  and  moments  im¬ 
parted  on  the  wing.  As  noted  by  Lorber^,  before  stall,  the 
wing-tip  vortex  lowers  the  effective  angle  of  attack,  re¬ 
duces  lift,  and  delays  the  onset  of  stall.  During  stall,  the 
interaction  of  tiie  wing-tip  vortex  prolongs  the  residence  of 
the  dynamic  stall  vortex  on  the  wing,  and  increases  the 
unsteady  increments  to  lift,  drag,  and  pitching  moment. 

The  wing-splitter  plate  configuration  was  pitched  from 
QO  to  60^  angle  of  attack  according  to  the  pitch  rate 
schedule  of  Eqn.  7  for  =0.2,  starting  from  the  initial 
solution  shown  in  Fig.  5.  A  constant  non-dimensional 
time  step  of  =0.002  with  two  subiterations  per  time 
step  was  maintained  throughout  the  pitch  cycle.  Figures  6, 
8,  and  9  graphically  depict  the  evolution  of  the  unsteady 
flow  field  on  the  lee  side  of  the  wing.  In  Fig  6,  contours  of 
constant  vorticity  magnitude  at  nine  spanwise  stations  (ylC 
=  0.05,  0.25,  0.50,  0.75,  1.0,  1.25,  1.50,  1.75,  2.0)  are 
plotted  by  slicing  the  three-dimensional  vorticity  field. 
Hg.  8  focuses  on  the  wing-tip  region,  displaying  vorticity 
magnitude  contours  at  five  chordwise  planes  (x/C  =  0.05, 
0.25,  0.5,  0.75,  1.0).  Hg.  9  displays  the  mireh  larger 
surface  vorticity  magnitudes  on  the  wing  and  splitter  plate. 
Fig.  6  provides  a  good  view  of  the  global  flow  features.  At 
1(P  angle  of  attack,  the  primary  separation,  originally  near 
mid  chord,  has  propagated  further  forward  and,  at  20®, 
reached  the  wing  leadiog  edge.  At  30®,  the  dynamic  stall 
has  formed  and  extends  spanwise  from  the  wing  root  to  a 
position,  ylC  1,5.  At  40®,  the  dynamic  stall  vertex  ap¬ 
pears  bowed,  with  maximum  displacement  above  the  wing 
midspan,  curving  downward  toward  the  wing  root  and 
wing  tip.  Near  the  center  span,  secondary  and  tertiary 
vortices,  characteristic  of  the  two-dimensional  pitching 
airfoil  are  evident.  At  the  tip,  the  dynamic  stall  vortex 
formation  was  suppressed  by  ^  wing-tip  vortex.  At  50® 
and  60®,  the  dynanaic  stall  vortex  near  the  wing  root  that 


previously  lagged  center  span  propagation  has  accelerated 
to  a  downstream  position  beyond  the  center  span.  As  the 
wing  was  pitched  to  60®  angle  of  attack,  the  position  of  the 
dynamic  stall  vortex  varied  uniformly  with  span  from  a 
"piimed"  position  at  the  wing  tip  to  a  position  weU 
downstream  at  the  wing  root.  Figure  7  presents  water 
tunnel  flow  visualization  results  at  40®,  51®,  and  59®  an¬ 
gles  of  attack  for  the  same  flow  conditions  (pitch  rate, 
pitch  axis,  and  Reynolds  number).  The  photographs  show 
that  the  simulation  accurately  predicted  the  principal  flow 
features.  In  particular,  the  photographs  confirm  the  bowed 
appearance  rf  the  dynamic  stall  vortex  at  40®,  the  accel¬ 
erated  convection  near  the  wing  root  at  50®  and  60®,  and 
the  vortex  ’’pinning"  at  the  wing  tip. 

As  mentioned  previously,  a  rounded  tip  was  used  for 
the  computations  while  a  square  tip  was  used  in  the  ex¬ 
periment.  The  primary  motivation  for  this  change  was  a 
desire  to  simplify  the  grid  generation  process  and  to  model 
the  computational  domain  with  a  single  grid  in  which  all 
interior  boundary  conditions  were  updated  implicitly.  In 
contrast  to  the  square  wing  tip,  the  separation  line  for  the 
wing-tip  vortex  was  determined  by  the  boundary  layer  dy¬ 
namics  and  not  fixed  at  the  sharp  wing-tip  edge.  Lorber^ 
experimentally  compared  rounded  and  square  wing-tip 
caps  and  found  that,  in  the  immediate  vicinity  of  the  wing 
tip,  the  square  tip  cap  produced  higher  unsteady  loads  pre¬ 
sumably  because  of  a  more  concentrated  wing-tip  vortex. 
Globally,  the  difference  in  the  tip  geometry  is  expected  to 
be  small.  La  any  event,  the  experimental  data  of 
Schreck^*^  did  not  include  surface  pressure  measurements 
beyond  the  yiC  =1.6  span  location  so  there  were  no 
measurements  to  compare  with  in  the  immediate  tip  re¬ 
gion.  Figure  8  shows  the  development  of  the  wing-tip 
vortex  and  its  effect  upon  the  dynaimc  stall  vortex.  At 
10®,  the  wing-tip  vortex  first  appeared  at  the  wing  trailing 
edge.  With  increasing  angle  of  attack,  it  propagated  for¬ 
ward  and  grew  in  size.  The  effect  on  the  dynamic  stall 
vortex  is  quite  evident.  At  30®  and  40®  angle  of  attack, 
the  flow  at  the  leading  edge  near  the  wing  tip  remained 
attached  with  the  wing-tip  vortex  effectively  terminating 
and  "pinning”  the  dynamic  stall  vortex  to  the  wing  tip. 
The  flow  did  not  fiilly  separate  at  the  first  chordwise  sta¬ 
tion  until  the  wing  reached  60®  angle  of  attack.  It  is  inter¬ 
esting  to  note  that  the  flow  visualization  results  of  Fig.  7 
also  show  the  dynamic  stall  vortex  terminated  inboard  of 
the  wing  tip  at  40®  as  predicted  by  the  computation. 

Figures  9  and  10  examine  the  surface  vorticity  and 
instantaneous  limiting  streamlines  on  the  wing  surface. 
Throughout  the  pitching  cycle,  the  highest  levels  of  vor¬ 
ticity  generation  occurred  at  the  wing  leading  edge.  For 
the  lower  angles  of  attack,  this  vorticity  was  relatively 
uniform  along  the  wing  span.  However,  at  30®  and  be¬ 
yond,  the  wing-tip  vortex  forced  the  leading-edge  flow  to 
remain  attached  near  the  tip,  increasing  the  intensity  of 
vorticity  generation  in  this  region.  At  20®,  numerous 


spanwise  vorticity  streaks  are  evident.  Comparison  with 
Fig.  10  shows  that  these  streaks  were  the  result  of  numer¬ 
ous  snrall-scale  separations  and  xeattachments  in  the  wing 
midspan  region.  At  30®,  a  strong  ridge  of  concentrated 
vorticity  can  be  seen  along  the  entire  span  of  the  wing 
from  the  wing-tip  leading  edge  to  the  wing  root.  This 
ridge  is  associated  with  the  separation  line  of  the  dynamic 
stall  vortex.  At  40®  and  50®,  this  ridge  weakened  near  the 
wing  root  as  the  dynamic  stall  vortex  propagated  away 
from  the  wing.  At  20®,  30®,  and  40®,  there  was  another 
strong  vorticity  region  associated  with  the  secondary 
separation  of  t^  wing-tip  vortex.  At  50®,  these  two  vor¬ 
ticity  concentrations  apparently  merged  and,  at  60®,  dis¬ 
appeared  as  the  flow  was  completely  separated  over  the 
entire  wing  surface. 

Instantaneous  limiting  streamlines  provide  a  final 
representation  of  the  three-dimensional  unsteady  flow 
structure  on  the  wing  surface  and,  through  inference,  the 
flow  field  above  the  surface.  Instantaneous  limiting 
streamlines  are  obtained  by  "freezing”  the  unsteady  flow 
field  at  a  particular  time  instant  and  integrating  the  sur¬ 
face  trajectories  of  massless  particles  introduced  at  a  limit¬ 
ing  distance  above  the  wing  surface.  The  frame  of  refer¬ 
ence  adopted  here  is  the  body  fixed  frame  that  rotates  with 
the  wing  about  the  1/4  chord  position.  Considerable  re¬ 
search^ has  been  devoted  to  the  identification  of  ciitical 
points  and  the  establishment  of  topological  rules  which 
govern  to  type  and  connectivity  of  critical  points  as  an  aid 
to  undei’standing  three-dimensional  flow  sepaiation. 
Interest  here  is  restricted  to  identifying  lines  of  separation 
(S)  upon  which  other  limiting  streamlines  coavejrge  md 
lines  of  reattachment  (R)  from  which  other  limiting 
streamlines  diverge.  In  Fig.  10,  the  limiting  streamlines 
are  plotted  on  a  flooded  surface  map  of  total  pressure  loss 
AFj  /Fj. ,  a  prcperty  directly  related  to  entropy  and,  conse¬ 
quently,  to  vOTticity.  The  limiting  streamlines  remforce 
many  of  the  flow  features  already  identified.  In  particular, 
the  limiting  streamlines  clearly  show  the  movement  of  the 
primary  separation  toward  the  leading  edge  and  the  con¬ 
current  formation  and  forward  pr<q)agation  of  the  wing-tip 
vortex.  At  30®  angle  of  attack,  the  curved  lines  of  separa¬ 
tion  and  reattachment  associated  with  the  "piniung"  of  the 
dynamic  stall  vortex  are  obvious.  The  wing-tip  vortex  and 
the  diordwise  vortical  structures  were  separated  by  a  line 
of  reattachment  that  extended  from  the  wing  tip  inward 
toward  the  trailing  edge.  Figure  10  clearly  shows  that  the 
regions  of  high  vorticity  production  and  total  pressure  loss 
evident  in  both  Hgs.  9  and  10  correspond  to  separation 
lines.  In  addition  to  these  primary  features,  both  the  wing- 
tip  vortex  and  chordwise  vortical  structures  were  charac¬ 
terized  by  very  complex  patterns  of  flow  separation  and 
reattachment. 

A  comparison  of  the  computed  and  experimental  pres¬ 
sure  coefficients  is  provided  in  Fig.  11,  The  comparison  is 
limited  to  spanwise  locations  less  than  ylC  =  1.6,  the  last 
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st^doii  at  which  experimental  measurements  were  taken. 
The  computation  conectiy  predicted  all  primary  flow  fea¬ 
tures  including  the  suction  pressure  ridge  seen  at  30^  inci¬ 
dence  along  the  span  of  the  wing  resulting  from  the  for¬ 
mation  of  the  dyiiamic  stall  vortex.  At  40^,  both  the  com¬ 
putation  and  experiment  show  the  downstream  propaga¬ 
tion  of  this  ridge  except  at  the  wing  tip  where  propagation 
was  delayed  and  the  suction  pressure  minimum  intensified 
due  to  the  interaction  of  the  dynamic-stall  and  wing-tip 
vortices.  At  50^  and  beyond,  the  computation  and  ex¬ 
periment  show  the  weakening  and  disappearance  of  the 
low  suction  pressure  region  from  the  wing  midspan  to 
wing  root  as  the  dynamic  stall  vortex  propagated  further 
away  from  the  wing.  Figure  12  presents  a  similar  com¬ 
parison  of  the  computational  and  experimental  normal 
force  coefficients.  The  two  results  are  in  good  agreement 
regarding  the  location  and  shape  of  the  broad  region  in 
which  normal  force  was  amplifled  due  to  vortex  interac¬ 
tion  as  well  as  a  much  smaller  normal  force  amplification 
near  the  wing  root  from  35®  to  45®  angle  of  attack. 
However,  the  computation  appears  to  undeipredict  the 
strength  of  the  suction  pressure  minimum  and  the  corre¬ 
sponding  normal  force  ampliflcation.  This  undeipredic- 
tion  is  most  likely  caused  by  inadequate  grid  resolution  or 
compressibility  e^cts.  A  two-dimensional  grid  resolution 
study  (Fig.  4)  showed  that  further  grid  enhancement  did 
not  significantly  change  the  intensity  of  the  pressure 
TniniTniim  or  increase  the  predicted  normal  force  coeffi¬ 
cient.  However,  it  is  well  known^^’^^  that  compressibility 
decreases  the  peak  suction  pressure  and  unsteady  force 
amplification  when  compared  to  values  for  very  low  Mach 
number  or  incompressible  flows.  For  the  computed 
freestream  Mach  number  of  0.2,  local  Mach  numbers  in 
the  wing-tip  region  near  the  leading  edge  wctc  well  into 
the  compressible  range  (M  ^  0.65),  thus  supporting  the 
hypothesis  that  compressibility  was  the  factor  most  re¬ 
sponsible  for  the  underprediction. 

Although  the  normal  force  amplification  near  the 
wing  tip  is  easily  traced  to  the  interaction  of  the  wing-tip 
and  dynamic-st^  vmtices,  the  source  of  the  smaller  nor¬ 
mal  force  local  maximum  near  the  wing  root  is  less  obvi¬ 
ous.  To  assess  viscous  wall  efrects  at  the  wing  juncture,  a 
computation  was  made  in  which  the  viscous  wall  boundary 
was  modeled  as  an  invisdd  surface  on  a  grid  that  did  not 
provide  spanwise  viscous  resolution  of  the  splitter  plate 
wall.  The  grid  was  otherwise  identical  to  the  baseline 
grid.  In  this  case,  the  normal  force  coefficient  showed  no 
evidence  of  amplification  near  the  wall.  Additionally,  the 
dynamic  stall  vortex  was  not  bowed  inward  toward  the 
wing  surface  in  the  juncture  region  as  in  the  viscous  wall 
case.  Both  features  appear  to  be  consequences  of  three- 
dimensional  viscous  interaction.  Further  evidence  for  this 
hypothesis  is  provide  in  Fig.  13  which  examines  the  ve¬ 
locity  field  in  the  wing  juncture  under  the  dynamic-stall 
vortex  at  a  time  corresponding  to  a  35®  pitch  angle.  This 
particular  angle  of  attadc  was  near  the  point  at  which  the 


wing  root  suction  pressure  was  a  minimum.  Both  three- 
dimensional  and  surface  tangent  velocities  are  shown  fts*  a 
chordwise  slicing  plane.  The  surface  tangent  velocities,  in 
particular,  show  a  well  developed  region  of  longitudinal 
vorticity  in  the  juncture.  Apparently,  this  additional  vor- 
tidty  in  the  juncture  region  increased  the  vortex  strength 
and  the  suction  pressure  under  the  vortex  in  comparison 
with  spanwise  stations  further  outboard.  This  longitudinal 
vorticity  also  pulled  the  dynamic  stall  vortex  inward  to¬ 
ward  the  wing  surface  near  the  juncture.  Both  effects 
would  act  to  intensify  the  local  suction  pressure  minimum 
and  increase  the  sectional  normal  force  coefficient  near  the 
wing  root.  The  same  phenomenon  was  also  observed  in 
earlier  calculations  that  modeled  the  wing  juncture  but  not 
the  wing  tip.^^ 

Figure  14  compares  the  overall  normal,  tangential, 
and  lift  force  coef^enls  generated  by  the  two-dimen¬ 
sional  airfoil  and  the  three-dimensional  wing.  As  ex¬ 
pected,  over  most  of  the  pitching  motion,  the  fimte  wing 
force  coefficients  were  less  than  those  of  the  airfoil  because 
the  tip  vortex  lowered  the  effective  angle  of  attack  near  the 
wing  tip.  Above  50®,  however,  because  the  wing-tip  vor¬ 
tex  "pinned"  the  dynamic  stall  vortex  to  the  tip,  both  the 
lift  and  normal  force  coefficients  of  the  fimte  wing  were 
higher  than  the  airfoil. 

Unsteady  Flow  •  Pitch  Rate  of  Q*  =  0.1 

A  second  computation  was  made  for  the  lower  pitch 
rate  cf  O'"  =0.1.  In  general,  the  lower  pitch  rates  are 
more  difficult  to  compute.^^  At  higher  pitch  rates,  the 
vigorous  forcing  motion  produces  an  inertially  dominated 
flow  that  is  less  sensitive  to  secondary  effects.  In  contrast, 
at  lower  pitch  rates,  relatively  small  inaccuracies  in  the 
prediction  of  critical  events,  when  magnified  over  the 
larger  time  scale,  can  result  in  significant  differences  in 
the  global  flow  features  as  well  as  flow  details. 

The  Q*  =  0.1  pitch  rate  case  used  the  same  grid  and 
initial  condition  as  the  £2**^  =  0.2  case.  The  pitch  rate 
schedule  was  again  governed  by  Eqn.  7  with  the  ^pro- 
priate  constants  for  the  lower  pitch  rate.  A  larger  time 
step,  At'^  =  0.(K)3  with  two  subiterations,  was  used  after  a 
comparison  with  the  smaller  time  step,  A/''’  =  0.(X)2,  re¬ 
vealed  identical  solutions  for  the  equivalent  two-dimen¬ 
sional  airfoil. 

Figure  15  provides  vorticity  contours  for  the  lower 
pitch  rate  at  selected  aoss  sectional  planes,  in  a  format 
identical  to  Fig.  6.  The  same  general  flow  features  are  ob¬ 
served.  The  dynamic  stall  vortex  was  again  pinned  to  the 
wing  leading  edge  until  very  late  in  the  pitch  cycle.  Near 
the  root,  the  vortex  propagation  was  first  retarded  and  then 
accelerated  as  the  pitch  cycle  progressed.  In  contrast  to 
the  higher  pitch  rate,  the  dynamic  stall  vortex  at  30®  was 
larger  and  had  propagated  further  downstream  to  the 
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midchord  posidon.  At  50°  incidence  and  above,  the 
original  dynamic  stall  vortex  has  propagated  beyond  the 
field  of  view  in  Fig.  15.  However,  a  new  vortex  is  evident 
at  50°  and  by  the  time  60°  has  been  reached,  this  vortex 
was  a  dominant  flow  feature  at  midspan.  A  comparison  of 
die  predicted  and  experimental  normal  force  coefficients  is 
provided  in  Hg  16.  Overall,  the  agreement  is  not  as  good 
as  in  the  higher  pitch  rate  case.  As  before,  the  magnitude 
of  the  normal  force  coefficient  was  underpredicted.  In 
addition,  the  normal  force  distribution  appeared  to  lead  the 
experiment  by  a  time  interval  corresponding  to  about  5°  in 
an^e  of  attack.  A  two-dimensional  grid  resolution  study. 
Fig.  17.  indicated  that  additional  grid  resolution  did  not 
significantly  change  either  the  magnitude  or  the  timing  of 
die  predict^  normal  force  time  history.  Compressibility  is 
again  the  most  likely  source  of  the  discrepancy.  The  most 
direct  support  for  this  view  was  provided  by  Patterson  and 
Lorber^^  who  compared  a  compressible  flow  solver  at 
M„  =  0.2  with  an  incompressible  flow  solver  for  the  same 
two-dimensional  airfoil  pitching  at  non-dimoisional  pitch 
rate  of  0.02.  They  found  that  the  compressible  flow  solver 
predicted  the  formation  of  die  dynamic  stall  vortex  earlier 
than  did  the  incompressible  code  and  tiiat  it  underpre- 
dicted  the  magnitude  of  the  resulting  force  coefficients. 
Both  observations  are  in  agreement  with  the  present  re¬ 
sults. 

Conclusions 

A  numerical  study  of  wing-tip  and  wing-juncture  in¬ 
teractions  for  a  finite  wing-splitter  plate  configuration 
pitching  through  extreme  angles  of  attack  was  conducted 
and  the  solutions  compared  with  corresponding  ex¬ 
perimental  wind  tunnel  and  water  tunnel  data.  At  the 
higher  of  two  pitch  rates  examined,  the  computations 
agree  well  with  available  experimental  data  and  provide 
new  insight  into  the  massively  stalled,  three-dimensional, 
unsteady  flow.  At  the  lower  pitch  rate,  the  computation 
again  correctly  predicted  the  overall  flow  features  but  pre¬ 
dicted  the  formation  and  propagation  of  the  dynamic  stall 
vortmt  prematurely.  Evidence  was  presented  to  suggest 
tiiat  this  discrepancy  resulted  from  tiie  difference  in  the 
computational  and  experimental  Mach  numbers.  The 
computational  fieestream  Mach  number  selected. 

=  0.2,  corresponds  to  a  generally  accepted  lower  limit 
for  reasonable  convergence  behavior  of  a  compressible 
flow  solver.  Since  compressibility  effects  were  significant 
at  this  Mach  number,  it  appears  that  an  incompressible 
flow  solver  would  be  required  to  rigorously  match  the  ex¬ 
perimental  test  conditions. 

Conversely,  most  unsteady  flows  involving  real  flight 
vehicles  are  compressible.  Ultimately,  the  unsteady  flow 
regime  cf  interest  to  maneuvering  aircraft  is  further  char¬ 
acterized  by  high  Reynolds  numbers,  turbulent  flow,  and 
three  dimensionality.  All  pose  significant  challenges  to 
experimental  and  numerical  investigations.  The  present 


study  has  focused  on  the  latter  characteristic-three-dimen- 
sional,  unsteady,  separated  flows.  The  computational  re¬ 
sults  provided  a  detailed  description  of  the  three-dimen¬ 
sional  interactions  at  the  wing  tip  and  wing-juncture  of  a 
generic  pitching  wing.  At  the  wing  tip,  the  numerical  re¬ 
sults  delineated  the  interaction  of  the  wing-tip  vortex  and 
the  dynamic  stall  vortex,  the  influence  of  the  wing-tip 
vortex  in  "pinning"  the  dynamic  stall  vortex  to  the  leading 
edge,  and  the  resulting  amplification  of  the  sectional  force 
r/^ffidentx  in  the  leading-edge  region.  At  the  wing  junc¬ 
ture,  the  numerical  results  provided  insight  into  the  three- 
dimensional  viscous  interaction  leading  to  a  smaller  local 
normal  force  maximum  near  the  wing  root. 
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Fig.  1  Angle  of  attack  histories  for  non-dimmsional 
pitch  rates  =  0.1,  =  0.2 
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Fig.  2  Grid  for  wing-splitter  plate  configuration  -  213x75x61 


Fig.  3  Temporal  Accuracy  for  equivalent  2D  pitching  airfoil,  Q'*'  =  0.2,  Re«  ^  56000 
a.  With  subiteration  b.  Without  subiteration 
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Fig.  7  Water  tunnel  flow  visualization,  Q*  =  0.2,  Re,  =  56000 
Dye  injection  at  leading  edge 
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Fig.  9  Surface  vorticity  flooded  contours  on  wing  and  splitter  plate 
Q*=  0.2,  Re,  =56000 
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Fig.  10  Instantaneous  limiting  streamlines,  =0.2,  Re.  =  56000 
Plotted  over  normalized  total  pressure  loss  flooded  contours 
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Fig.  15  Vorticity  contours  at  9  spanwise  stations,  Q*  =  0.1,  Re_  =  56000 
y  /  c  =  {0.05, 0.25, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0} 
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Fig.  16  Comparison  of  computational  and  experimental  normal  force  coefficients 
Q*=0.1,  Re„- 56000 


Fig.  17  Spatial  accuracy  for  equivalent  2D  airfoil,  Q.*  =  0.1,  Re.  =  56000 
Grid  Sizes:  213x75  425x149 
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Abstract 

Flow  visualization  data  over  the  suction 
surface  of  a  65°  swept  delta  wing  were  recorded  at  a 
30°  incidence  angle  for  a  range  of  static  roll  angles  (- 
42°  to  42°)  and  Reynolds  number  of  9.0  x  10^.  The 
data  were  reduced  to  provide  information  on  the 
vortex  relative  position  and  the  burst  point  location. 
The  results  were  compared  to  previously  recorded 
rolling  moment  data,  to  provide  an  understanding  of 
the  relationship  between  flow  behavior  and 
aerodynamic  loading.  The  combination  of  roll  about 
the  body  axis  and  an  axis  incidence  angle  with  the 
flow  induced  a  sideslip  angle  relative  to  the  wing. 
Combining  the  sideslip  angle  with  the  nominal 
leading  edge  sweep  of  the  wing  provided  an  effective 
leading  edge  sweep  angle.  Using  the  effective  sweep 
angle,  along  with  the  previous  data  for  iinyawed  delta 
wings,  the  flow  behavior  over  the  either  side  of  the 
wing  was  predicted.  Extending  this  analysis 
suggested  the  presence  of  non-zero  trim  positions  in 
roll,  based  on  wing  movement  in  and  out  a  stalled 
regime.  These  trim  positions  have  been  recorded  in 
free-to-roll  experiments. 

Symbols 

s  =  local  semispan 

X  =  chordwise  position 

y  =  spanwise  position 

z  =  position  normal  to  the 

wing  surface 

A  =  sweep  angle 

a  =  angle  of  attack  between 

the  wing  and  the  freestream 
A  =  increment 

a  =  incidence  angle  between 

the  root  chord  and  the  wind  tunnel  floor 
(|)  =  roll  angle  about  the  body 

axis  subscripts 
Subscripts 

bp  =  burst  point 

eff  =  effective 

0  =  nominal 


Introduction 

The  flow  over  the  delta  wing  planform 
has  been  of  interest  since  the  1950’s  due  to  the  delta 
wing  design  used  on  high  speed  aircraft.  The  delta 
wing  provides  a  unique  type  of  flow  field  with  the 
formation  of  large  helical  vortices  over  the  upper 
wing  surface.  This  flow  structure  allows  the  delta 
wing  to  maintain  lift  to  high  angles  of  attack  when 
compared  to  the  conventional  straight  wing. 

With  the  current  emphasis  on  improving 
aircraft  maneuverability,  the  flight  envelope  is  being 
pushed  to  higher  limits,  even  into  the  post  stall 
regime.  In  these  flight  regimes  the  stall 
characteristics  of  the  wing  become  important.  For  j 
the  delta  wing,  a  precursor  to  stall  on  the  delta  wing 
is  the  movement  of  the  burst  point  to  a  position  over 
the  wing.  Therefore  the  position  of  the  vortex  burst 
point  can  be  an  indicator  for  wing  performance. 

To  evaluate  the  fluid  dynamics  and 
aerodynamic  loading  of  the  delta  w^ing  undergoing 
high  angle  of  attack  maneuvers,  the  Canadian 
National  Research  Council  (NRC)  developed  a 
comprehensive  test  progran^for  testing  a  65°  swept 
delta  wing  in  a  wind  tunnel  .  Force  balance,  surface 
pressure,  and  flow  visualization  data  were  recorded 
for  the  flow  over  a  delta  wing,  which  was  held  at  an 
incidence  angle  of  30°and  was  forced  statically  and 
dynamically  through  a  range  of  roll  angles  about  the 
body  axis. 

The  forcj  balance  results  have  been 
discussed  previously  .  In  this  paper  the  flow 
visualization  will  be  discussed  with  emphasis  on  the 
behavior  of  the  burst  point  and  the  changing  structure 
of  the  vortex  flow  for  different  geometric 
orientations.  This  presentation  is  limited  to  the  static 
results.  The  results  were  analyzed  according  to  the 
effective  sweep  angles  on  leading  edges  of  the  wing. 
The  effective  sweep  angle  was  determined  by  the 
relative  orientation  of  the  leading  edge  to  the 
fi-eestream.  This  analysis  provides  some  physical 
interpretation  of  the  rolling  moment  results  that  have 
been  previously  published. 


Methods 


Results 


The  delta  wing  model  the  tunnel  set 
up  have  been  documented  by  Hanff  (Fig.  1).  The 
model  was  65°  swept  back  flat  plate  model  with  a  25 
°  knife  edge  bevel  and  a  root  chord  length  was  24.48 
inches.  The  root  chord  axis  of  the  model  was  held  at 
a  30°  incidence  zmgle  relative  to  the  floor  of  the  wind 
tunnel.  Roll  angles  were  measured  about  this  axis. 
Kerosene  smoke  was  injected  forward  of  the  apex 
and  was  entrained  into  the  two  vortices  of  the  suction 
surface  of  the  wing.  A  movable  set  of  optics  was 
designed  to  produce  a  laser  sheet  normal  to  the  root 
chord  axis.  The  laser  sheet  was  moved  along  the 
model  from  the  trailing  edge  to  the  apex  illuminating 
cross  section  cuts  of  the  conical  vortices. 

The  flow  visualization  was  reduced 
using  an  im^e  analysis  program  documented  by 
McLaughlin  .  The  image  analysis  program,  first, 
identified  the  wing  and  its  roll  orientation.  Then  the 
position  and  size  of  both  leading  edge  vortices,  as 
well  as  the  burst  condition,  were  determined  for  each 
analyzed  frame  of  the  video  data.  Vortex  size  was 
determined  by  the  area  of  the  smoke  illuminated  over 
either  side  of  the  wing.  The  vortex  position  was 
initially  estimated  as  the  centroid  of  that  illuminated 
smoke  region. 

The  burst  condition  of  the  vortex  was 
determined  by  the  presence  or  absence  of  coherent 
core  flow,  denoted  by  a  small  darkened  circle  near 
the  center  of  the  illuminated  smoke  region  (Fig.  2, 
Frame  a.).  Reasons  for  the  dark  appearance  of  the 
core  flow  are  discussed  in  Nelson  .  If  the  core  was 
detected,  the  core  center  was  used  to  estimate  the 
position  of  the  vortex. 

Other  res^^;pjiers  have  discussed  the 
nature  of  vortex  burst  ,  either  as  the  sudden 
expansion  of  the  core  (bubble  breakdown)  or  as  a 
sudden  core  trajectory  change  from  a  line  to  a  helix 
(helical  breakdown).  In  either  case  the  breakdown 
was  manifested  in  the  flow  visualization  as  the 
expansion  of  the  darkened  core,  followed  by  its 
disappearance  as  the  core  flow  diffuses.  Once  the 
general  spatial  region  of  core  diffusion  was 
determined  from  the  image  analysis  program,  the 
region  was  reviewed  manually  to  determine  where 
core  expansion  initiated.  The  burst  location  was 
recorded  as  the  chordwise  position  where  the 
darkened  core  began  to  expand.  A  similar  burst 
identificatjpn  criteria  was  used  previously  by 
Kegelman  .  Since  the  burst  location  is  not  fixed,  this 
technique  determined  the  most  upstream  burst 
location. 


Examples  of  the  flow  visualization  data 
(viewed  from  the  front^)are  shown  in  Fig.  2  for  roll 
angles  from  -14°  to  42°.  These  results  were  recorded 
for  the  60%  chord  position.  Since  the  vortex  remains 
conical  over  the  surface,  the  relative  spanwise 
position  of  the  vortex  (y/s)  was  nearly  constant,  so  a 
single  position  was  representative  of  the  vortex  flow 
position  over  the  whole  wing.  The  frames  show  the 
smoke,  entrained  in  the  vortices  on  either  side  of  the 
rolled  wing  and  illuminated  by  a  laser  sheet 
perpendicular  to  the  roll  axis.  The  wing  can  be  seen 
highly  illuminated  at  the  bottom  of  each  frame.  The 
hump  in  the  center  of  the  wing  is  the  housing  for  the 
force  balance  and  mounting  apparatus. 

The  images  show  that  flow  over  the 
leeward  (upper)  side  produced  a  more  compact 
vortex  with  a  more  aft  burst  location.  The  aft  burst 
location  was  indicated  by  the  presence  of  a  core  in 
the  leeward  vortex  only.  At  the  higher  roll  angles  the 
core  on  the  leeward  side  was  present,  but  was  washed 
out  by  the  high  intensity  of  the  illuminated  smoke. 
The  leeward  vortices  also  exhibited  vortical 
ornamentation,  evidenced  by  the  wavy  outline  of  the 
illuminated  smoke.  The  ornamentation  is  the 
discretization  of  the  outer  shear  layer  into  a  s^jies  of 
vortices  similar  to  that  shown  by  Gad-el-Hak  (Fig. 

2, 2nd  column). 

For  roll  angles  greater  than  7°  the  flow 
on  the  windward  wing  no  longer  exhibited  the 
organized  outboard  flow  along  the  surface,  indicative 
of  a  vortex.  This  flow  was  appropriately  referred  to 
as  separated  rather  than  as  vortical  flow.  As  the  roll 
angle  increased  past  14°,  the  smoke  on  the  windward 
side  did  not  fill  the  separated  region  completely,  but 
remained  in  the  inboard  part  of  the  region.  The 
height  of  the  separated  region  could  still  be 
estimated,  however.  From  the  estimation  the 
separation  height  reached  a  maximum  at  21°  roll, 
declining  with  higher  angles. 

The  trend  of  having  a  small,  coherent 
vortex  on  the  leeward  side  and  a  burst  vortex  or 
separated  flow  on  the  windward  side  was  consistent 
in  all  the  frames.  The  negative  roll  angles  were 
included  in  Fig.  2  to  show  the  symmetry  of  the  tests. 
The  lack  of  antisymmetry  between  the  flow  fields  for 
14  and  -14°  suggests  that  the  roll  angles  may  have 
had  a  slight  offset  in  the  zero  roll  position,  possibly 
because  of  support  interference^  \ 

From  the  flow  visualization  data  the 
vortex  positions,  as  well  as  the  vortex  burst  point 
locations,  were  tracked.  The  relative  spanwise 
positions  of  the  vortex,  y/s,  at  the  60%  chord 
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position  are  shown  in  Fig.  3  for  the  range  of  roll 
angles  from  -42°  to  42°.  As  discussed  above,  these 
data  represent  the  spanwise  core  positions  when  a 
vortex  core  was  apparent.  When  a  vortex  core  was 
not  apparent,  i.  e.  the  vortex  was  burst  or  the  flow 
was  separated,  the  data  were  the  spanwise  centroid 
positions  of  the  illuminated  smoke.  The  spanwise 
position  is  shown  because  the  far  inboard  spanwise 
position  corresponded  with  separation  of  the  flow. 
Therefore  these  data  provide  a  indication  of  the 
changing  structure  of  the  flow. 

The  data  on  the  upper  half  of  the  graph 
correspond  to  the  right  wing  shown  in  Fig.  2,  while 
the  lower  half  data  correspond  to  the  left  wing.  The 
general  trend  of  the  vortex  position  was  to  move 
inboard  as  the  wing  moved  into  a  more  windward 
position.  The  vortices  reached  a  most  inboard 
position  near  35°  and  -35°  roll  for  the  left  and  right 
wings,  respectively.  A  lack  of  antisymmetry  was 
-apparent  at  (|)  =  0°  with  the  vortex  on  the  left  moving 
slightly  further  inboard.  This  may  also  be  a  result  of 
support  interference^*  as  suggested  earlier. 

Several  regions  have  been  labeled  on  the 
graph  to  show  the  relation  between  vortex  position 
and  burst  condition.  When  the  burst  point  was  aft  of 
the  trailing  edge  or  over  the  surface,  there  is  a  near 
linear  relationship  between  the  roll  angle  and  vortex 
position.  When  the  burst  point  reached  the  apex,  the 
character  of  the  flow  changed  from  vortex  to 
separated  flow.  Also  the  smoke  did  not  reach  the 
outboard  separated  region,  as  was  described  above. 
Both  of  these  factors  contributed  to  the  inboard 
movement  of  the  position  data. 

The  burst  point  data  for  either  wing  are 
shown  in  Fig.  4  with  the  relative  chordwise  burst 
location  plotted  versus  roll  angle.  Similar  to  the 
position  data,  these  data  show  the  inverse  behavior  of 
the  burst  point  data  from  one  wing  to  the  other.  As 
the  roll  angle  increased  the  burst  point  on  the  right 
wing  moved  aft  while  the  burst  point  moved  forward 
on  the  left  wing.  The  intersection  of  the  two  curves 
occurred  at  a  slightly  negative  roll  angle. 

The  graph  shows  that  the  burst  point  was 
only  over  the  surface  for  small  portion  of  the  test 
range.  When  the  burst  point  was  at  or  aft  of  the 
trailing  edge,  the  data  show  the  burst  point  at  the 
trailing  edge.  Since  the  laser  sheet  did  not  move  past 
the  end  of  the  model,  it  was  not  possible  to  track  the 
burst  point  aft  of  the  wing.  The  data  show  the  burst 
point  at  the  apex  for  all  cases  where  no  vortex  core 
was  apparent  over  the  wing. 

Discussion 


In  order  to  discuss  these  results,  a 
conceptual  model  for  delta  wing  flow  is  suggested. 
For  the  model,  it  is  postulated  that  the  organization  of 
vortices  over  the  delta  wing  is  a  function  of  the 
sweep  back  angle,  the  angle  of  attack  and  the 
influence  of  the  vortex  on  the  opposite  side  of  the 
wing.  The  sweep  angle  and  the  angle  of  attack 
determine  the  orientation  of  the  leading  edge  relative 
to  the  flow.  The  orientation  should  determine  the 
vorticity  flux  fi*om  the  leading  edge  into  the  delta 
wing  vortex.  The  sweep  angle  and  the  angle  of 
attack  also  determine  the  region  shadowed  from  the 
freestream  in  which  vorticity  can  accumulate.  The 
sweep  angle  and  angle  of  attack  are  important 
parameters  in  the  determination  of  vortex 
organization. 

The  influence  of  the  opposite  vortex,  for 
delta  wings  with  sweep  angle  less  than  70°,  would  be 
to  impose  a  downward  velocity  component  on  the 
vortex.  At  higher  sweep  angles  the  vortex-vortex 
interactions  become  more  complicated.  The 
magnitude  of  the  influence  is  determined  by  the 
strength  and  thg  proximity  of  the  two  vortices.  Roos 
and  Kegelman  have  shown  that  the  spanwise 
position  of  the  vortex  core  is  roughly  60%  of  the 
semispan  from  the  root  chord,  a  function  of  the 
sweep  angle.  For  a  Iqwer  swept  wing  the  vortices  are 
farther  apart.  Nelson  has  shown  that  the  circulation 
contained  in  the  vortex  is  also  proportional  to  the 
sweep  angle.  Tlie  lower  swept  wing  produces  higher 
circulation  values.  Since  these  two  effects  act 
counter  to  one  another,  as  a  first  approximation  for 
analysis,  the  influence  of  the  opposite  vortex  is 
assumed  to  be  constant.  Changes  in  vortex 
development  are  therefore  assumed  to  be  only  a 
function  of  the  sweep  angle  and  angle  of  attack.  This 
appears  to  be  true^^om  the  half  models  results 
reported  by  Roos  which  showed  excellent  flow 
similarity  to  full  delta  wing  models  with  sweep 
angles  of  60°  and  70°. 

The  implication  for  a  yawed  delta  wing 
is  that  the  flow  will  develop  according  to  the 
effective  sweep  of  the  wing  leading  edge.  The 
effective  sweep  would  be  a  combination  of  the 
nominal  sweep  of  the  wing  and  the  sideslip  angle. 
The  flow  over  the  windward  wing  would  develop  as 
if  the  sweep  angle  were  reduced  and  vice  versa  for 
the  leeward  wing.  The  trend  for  this  behav^^r  was 
seen  in  previously  reported  burst  point  data  over 
yawed  delta  wings.  The  burst  point  on  the  windward 
wing  moved  forward  with  an  increase  in  yaw, 
consistent  with  a  decrease  in  the  leading  edge  sweep 
angle. 
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In  the  current  tests,  the  incorporation  of 
roll  onto  the  delta  wing  at  an  angle  of  attack  induced 
a  sideslip  angle  while  reducing  the  angle  of  attack. 
Equations  for  calculating  the  reduced  angle  of  attack 
and  the  sideslip  aij^le  as  functions  of  (j)  and  a  are 
found  in  Ericsson  . 

a((t))  =  tan“^ftan  a  cos 

AA  (([))  =  tan"Htan  a  sin 

The  effective  sweep  angle  of  the  leading  edge  would 
be  determined  by  a  combination  increment  of  the 
sweep  angle  due  to  the  induced  side  slip  angle  and 
the  nominal  sweep  angle,  A. 

windward  wing:  A^ff  =  Ao  -  AA 

leeward  wing  :  A^ff=  Aq  +  AA 

Using  the  calculated  effective  sweep  angle  and  angle 
of  attack,  combined  with  previously  reported  burst 
point  data  the  burst  point  locations  for  the  rolling 
delta  wing  could  be  predicted. 

In  order  to  address  the  accuracy  of  this 
,simple  estimation  tec^^jcjji^  burst  data  were  taken 
from  several  sources  for  a  range  of  sweep 

angles  and  angles  of  attack.  These  data  are  plotted  in 
Fig.  5.  These  data  were  collected  over  as  wide  of  a 
range  of  sweep  angles  as  possible.  On  the  low  end, 
researchers  have  reported  no  evidence  of  a  coherent 
vortex  over  the  surface  for  a  sweep  angle  of  45°  or 
less,  regardless  of  angle  of  attack.  The  upper  end  is 
obviously  limited  to  less  than  90°.  The  angle  of 
attack  range  was  determined  by  the  presence  of  the 
burst  point  over  the  surface.  For  any  model  the  burst 
point  was  over  the  surface  for  10°  to  15°  angle  of 
attack  range.  For  the  85°  delta  wing  the  burst  point 
would  reach  the  apex  at  about  70°  angle  of  attack. 

The  different  symbols  on  the  graph 
represent  the  different  sweep  angles.  The  separation 
of  the  symbols  shows  a  shift  to  higher  angles  of 
attack  with  increased  sweep  angle.  The  scatter  in  the 
data  was  ^e  to  differences  in  the  model  shape. 
Kegelman  has  shown  that  this  scatter  can  be  reduced 
by  comparing  wings  with  the  same  leading  edge 
shape,  j^he  flagged  data  points  were  taken  from 
Wentz  .  Their  test  model  had  a  flat  surface  and  a 
convex  surface.  Data  taken  using  the  convex  surface 
as  the  upper  surface  produced  a  definite  change  in  the 
burst  behavior.  Therefore  the  data  from  the  convex- 
up  tests  are  flagged.  All  other  data  were  taken  with 
flat  plate  models. 

Visual  inspection  of  the  data  indicated 
that  the  burst  point  location  had  a  similar  function  of 
the  angle  o^^ttack,  regardless  of  the  sweep  angle. 
Lamboume  had  previously  shown  that  the  burst 
point  location  data,  under  different  sweep  conditions, 
collapsed  on  a  single  curve  when  plotted  against  cos" 
l(cos  a  sin  A).  Plotting  the  collection  data  from  Fig. 


5  against  this  function  did  not  produce  the  same  level 
of  collapse  that  Lamboume  had  reported  so  another 
function  was  developed  using  curve  fitting 
techniques. 


0.00741A- 0.287  ^ 

a  -(0.00866A-0.3573)J 
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It  is  stressed  that  the  model  is  not  based  on  a  physics 
of  the  flow  field.  Rather,  a  function  for  the  burst 
point  location  was  sought  to  model  sweep  and  attack 
angle  effects  in  order  to  compare  the  nominal 
behavior  of  the  flow  over  delta  wings  to  that  under 
yawed  conditions.  The  results  of  applying  this 
function  to  the  A  and  a  conditions  for  Fig.  5  are 
shown  in  Fig.  6. 

The  burst  point  fixation  was  applied  to 
the  data  taken  from  McKeman  in  Fig.  7.  The  burst 
point  location  data  were  recorded  for  a  70°  swept 
wing  over  a  range  of  angles  of  attack  and  sideslip 
angles.  The  data  were  presented  for  the  windward 
side  of  the  wing  only.  The  burst  point  location  was 
estimated  by  substituting  A^ff  in  the  burst  point 
function.  The  results  show  excellent  agreement 
between  the  predicted  vortex  burst  position  and  the 
measured  values  suggesting  that,  at  least  in  this 
limited  range,  the  flow  field  on  either  side  of  the 
delta  wing  can  be  determined  by  the  effective  sweep 
angles  and  angle  of  attack. 

The  same  function  was  applied  to  the 
current  data  for  a  delta  wing  with  a  combination  of 
pitch  and  roll.  The  comparison  between  the 
measured  data  and  that  predicted  by  the  function  are 
shown  in  Fig.  8.  For  this  comparison  the  data  from 
Fig.  4  were  replotted  versus  the  effective  sweep  angle 
of  each  side  of  the  wing.  For  these  cases,  the  burst 
point  function  predicted  the  experimental  data  within 
20%  chord  for  effective  sweep  angles  below  70°.  A 
20%  chord  change  in  the  burst  location  with  different 
leading  edge  shapes  was  not  unreasonable,  since 
most  of  the  reference  data  were  collected  on  models 
with  flat  upper  surfaces,  and  the  upper  surface  of  this 
delta  wing  model  was  beveled  at  the  edges  Also, 
it  was  difficult  to  determine  burst  locations  forward 
of  20%  chord  because  the  bright  illumination  of  the 
dense  smoke  in  that  region  washed  out  the  small, 
dark  core  flow.  It  is  noted  that  the  excellent 
comparison  with  McKeman’s  data,  for  the  windward 
wing  only,  was  also  for  effective  sweep  angles  less 
the  70°. 


These  results  show  that  a  burst  point 
equation  based  on  the  nominal  wing  configurations 
provided  accurate  burst  point  estimates  for  effective 
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sweep  angles  less  than  70°.  Above  70°  the  trend  was 
predicted  but  the  error  increases.  This  may  be 
indicative  of  the  scatter  in  the  data  on  which  the 
predictive  equation  was  based  or  it  may  indicate  a 
breakdown  in  the  premise  of  opposite  vortex 
influence  at  higher  sweep  angles. 

The  first  possibility  could  be  resolved 
with  a  parametric  set  of  flow  visualization  tests  over 
delta  wings  with  common  leading  edge  shapes  and  a 
wide  range  of  sweep  back  angles.  As  for  the  second 
possibility,  it  has  been  shown  that  the  vortex  sheets 
separating  from  both  sides  of  the  wing  come  in 
contact  prior  to  vortex  breakdown'^  for  sweep  angles 
greater  than  75°.  For  angles  of  attack  above  initial 
vortex  sheet  contact,  the  vortex  cores  can  be  forcibly 
displaced  off  the  surface.  This  is  presumably  an 
effect  of  the  interaction  between  the  vortices  due  to 
their  proximity.  This  core  displacement  may  explain 
stall  for  higher  sweep  angles  prior  to  the  burst  point 
reaching  the  apex^^.  For  the  case  of  the  yawed  wing, 
with  lower  nominal  sweep  than  the  effective  value, 
opposite  vortices  would  not  have  the  same  proximity 
and  would  not  interact  in  the  same  manner.  This 
may  be  why  this  predictive  technique  becomes  less 
accurate  for  effective  sweep  angles  greater  than  70°. 

This  flow  field  dependence  on  the 
effective  sweep  and  attack  angles  was  used  to 
estimate  the  burst  point  behavior.  The  roll  angle 
range,  within  which  the  burst  point  was  over  the 
surface,  was  only  a  small  portion  of  the  tested  roll 
angle  range.  In  the  remainder  of  the  range  this 
effective  sweep  angle  concept  allows  for  speculation 
about  the  static  loading  behavior.  In  particular, 
trends  in  the  rolling  moment  data  could  be  estimated. 

Using  the  effective  angle  estimates 
discussed  above,  the  effective  leading  edge  sweep 
and  attack  angles  were  calculated  for  a  65°  swept 
wing  with  a  30°  nominal  angle  of  attack  and  a  range 
of  roll  angles.  The  results  are  plotted  as  the  solid 
curve  in  Fig.  9.  The  upper  half  of  the  curve 
represents  the  leeward  wing  and  the  lower  half,  the 
windward  wing. 

Also  included  on  the  graph  (open 
circles)  are  th|2j^1j|ljQangles  for  a  range  of  delta  wing 
sweep  angles  ’  .  These  data  were  taken  from  the 

flow  over  complete  delta  wing  models.  The  flat  plate 
rectangular  wing  was  assumed  to  stall  at  about  7° 
angle  of  attack.  The  dashed  line  on  the  graph 
represents  an  average  stall  angle  for  the  different 
sweep  angles  and  divides  the  graph  into  stalled  and 
unstalled  regions.  The  effective  sweep  angle  analysis 
would  predict  that  movement  of  the  windward  wing 
into  the  stalled  region  would  cause  the  flow  over  that 
wing  to  stall,  along  with  a  loss  of  lift. 


Finally,  a  line  was  included  at  the  90° 
effect! ve^weep  position  and  is  denoted  as  wingtip. 
Erickson  has  shown,  on  highly  swept  wings, 
imparting  a  sideslip  angle  to  the  wing  can  cause  the 
vortex  to  move  off  of  the  surface.  The  90°  position 
provided  a  limit  for  this  type  of  behavior. 

The  results  on  Fig.  8  can  be  related  to 
the  rollin^poments  in  the  following  manner. 
Eamshaw  showed  that  the  lift  curve  slope  for  a 
lower  swept  wing  is  higher  than  for  a  higher  swept 
wing  if  neither  has  stalled.  On  the  figure  the 
windward  wing,  with  the  lower  effective  sweep 
angle,  would  produce  more  lift  then  the  leeward  wing 
at  the  same  angle  of  attack.  If  the  moment  arms 
about  the  root  chord  are  assumed  to  be  nearly  equal, 
this  loading  behavior  would  provide  a  stabilizing 
rolling  moment  and  a  trim  position  at  0°  of  roll. 

Movement  of  the  windward  wing  into 
the  stalled  region  would  cause  the  lift  produced  by 
the  windward  wing  to  decrease.  The  leeward  lift 
could  then  dominate  and  produce  an  unstable  rolling 
moment.  A  limit  factor  to  instability  could  be  the 
movement  of  the  leeward  wing  past  the  90°  effective 
sweep  angle,  which  would  cause  the  vortex  to  move 
off  the  surface,  reducing  the  lift  due  to  the  presence 
of  the  vortex.  The  result  would  be  an  increase  in  the 
correcting  moment  back  to  0°  roll.  Similarly 
movement  of  the  windward  wing  back  into  the 
unstalled  region  at  high  roll  angles  should  again 
produce  a  correcting  moment  back  toward  0°. 

For  the  test  data  (A  ==  65°,  o  .30°)  the 
effective  angle  curve  is  in  the  unstalled  region  at  0° 
roll.  As  described  above  a  statically  stable  roll  point 
would  be  expected  at  0°.  This  point  was  2 

demonstrated  in  free-to-roll  tests  reported  by  Hanff  . 

The  effective  angle  curve  intersects  the 
estimated  stall  curve  at  approximately  10°  roll.  If  the 
effective  angle  analysis  holds,  a  destabilizing  rolling 
moment  would  be  produced  for  higher  roll  angles. 
When  the  roll  angle  was  increased  to  50°,  however, 
the  windward  wing  returns  to  the  unstalled  region 
and  a  stabilizing  rolling  moment  would  again  be 
produced.  A  zero  rolling  moment  and  non-zero 
stable  roll  position  must  then  occur  somewhere 
between  10°  and  50°  of  roll.  Rolling  moment  and 
free-to-roll  data  for  a  65°  delta  wing,  reported  by 
Hanff  showed  this  general  trend  in  the  moment 
data  with  non-zero,  static  stability  (critical )  points  at 
±21°. 

Further  resets,  also  for  the  65°  model, 
were  reported  by  Hanff  for  attack  angles  of  20°,  25° 
,  35°  and  40°.  The  corresponding  stability  points 
were  at  0°,  1 1.5°  | ,  1 1 1°  | ,  and  0°,  respectively. 
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Using  the  current  analysis,  effective  angle  diagrams 
were  plotted  for  these  cases  in  Figs.  1 1-14. 

In  Fig.  10  the  entire  effective  angle 
curve  was  well  out  of  the  stalled  region  for  all  cases. 
Only  a  stability  point  at  0°  would  be  expected  for  this 
case,  consistent  with  the  Hanff  s  experimental  results. 
The  same  results  would  be  expected  from  Fig.  11. 
Hanff  reported  a  slight  deviation  from  that  point. 

For  an  attack  angle  of  35°  (Fig.  12),  the 
0°  roll  position  on  the  curve  has  moved  into  the 
stalled  region.  No  rolling  moment  would  be  expected 
under  these  symmetric  conditions,  however,  a 
increase  in  the  roll  angle  moves  the  leeward  wing 
into  the  unstalled  region  producing  an  destabilizing 
moment.  No  zero  stability  point  was  expected  for 
this  case.  This  analysis  would  again  suggest  the 
presence  of  a  non-zero  stability  point.  It  does  not 
provicj^  the  1 1 1°  |  trim  point^yalue  reported  by 
Hanff  and  again  by  Jenkins  .  This  analysis  would, 
at  first,  might  suggest  that  the  stability  point  would 
be  at  a  position  greater  than  the  21°  value  reported 
for  an  attack  angle  of  30°  since  the  return  to  the 
unstalled  region  occurs  at  a  higher  roll  angle.  Other 
considerations,  however,  would  have  to  be  the  roll 
angle  where  the  leeward  wing  passes  90°  effective 
sweep  and  that  the  wing  is  stalled  for  0°  roll.  Also 
the  flow  behavior  on  either  side  of  the  wing ,  once 
one  side  has  stalled,  can  no  longer  be  considered  as 
independent.  These  factors  make  it  more  impossible 
predict  where  the  trim  point  would  be  located  using 
this  simple  analysis. 

In  Fig.  13  (a  =  40°)  the  wing  is  well  into 
the  stalled  region  for  ^  =  0°.  In  this  case  the  0°  roll 
position  may  become  neutrally  stable.  Again  it 
would  be  difficult  to  make  accurate  estimates  of  any 
non-zero  statically  stable  roll  positions. 

Increasing  the  wing  sweep  angle  could 
move  the  effective  angle  curve  up  and  out  of  the 
stalled  region.  Roljjng  moment  data  from  free-to-roll 
tests  by  Thompson  for  a  80°  sweep  wing  shows  no 
signs  of  the  non-zero  trim  points.  Increasing  the 
sweep  angle,  however,  moves  the  models  into  a 
regime  where  wing  rock  can  occur,  a  dynamic 
problem  which  will  not  be  addressed  here. 

For  sweep  angles  less  than  70°,  the 
model  predicts  three  stable  point  cases;  one,  two  or 
three  stable  roll  positions.  One  point  can  occur  at  0° 
roll  when  the  effective  angle  curve  is  outside  the  stall 
region.  Two  statically  stable  points  can  occur  when 
the  0°  point  is  slightly  in  the  stalled  region.  Three 
points  can  occur  when  the  0°  roll  position  is  slight 
outside  the  stalled  region  and  the  effective  angle 
curve  passes  through  the  stalled  region. 


Conclusions 

The  results  show  that  the  vortex 
behavior  on  either  side  of  a  delta  wing  can  be  related 
to  the  effective  sweep  and  the  attack  angles  of  the 
leading  edge.  Applications  of  an  empirical  equation 
for  the  burst  point  based  on  imyawed  data  showed 
good  comparisons  to  yawed  data  for  effective  sweep 
angles  below  70°.  The  differences  may  be  due  to  the 
influence  of  the  vortex  on  the  opposite  wing.  This 
issue  could  be  addressed  in  a  parametric  study 
comparing  the  flow  over  yawed  and  unyawed  wings. 
Comparison  of  existing  data  from  several  authors  is 
difficult  due  to  the  sensitivity  of  the  flow  to  model 
leading  edge  shape. 

The  empirical  equation  was  not  based  on 
physical  phenomenon  but  rather  was  a  curve  fit  for 
comparison  between  yawed  and  unyawed  data.  A 
similar  formula  based  on  the  fluid  dynamic  should  be 
developed.  Such  an  equation  could  then  be  extended 
to  apply  to  the  dynamic  data  for  which  the  reported 
data  was  a  baseline.  The  importance  of  the  effective 
sweep  and  attack  angles,  shown  here,  may  help  to 
develop  this  equation. 

Finally  the  non-zero  stability  points, 
shown  by  Hanff,  have  been  discussed  from  a 
different  point  of  view.  The  simple  analysis  using 
effective  sweep  and  attack  angles  provides  and 
explanation  for  the  general  trends.  An  extended 
study  to  try  to  define  the  sweep  and  attack  angle 
limitations  of  this  behavior  and  to  predict  these  points 
would  be  valuable. 
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Figure  2.  Caricature  of  llie  flow  visualization  results. 


Figure  3.  How  visualization  images;  recorded  at  the  60%  chord  position,  o  -  30  , 
(j)  =  -14°  to  42°. 
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Figure  4.  Vortex  relative  spanwise  positions;  recorded  at  the  60%  chord  position, 
a  =  30°,  (|)  = -42°  to  42°. 
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Figure  5.  Vortex  burst  point  relative  chordwise  positions,  a  =  30°,  (|>  =  -42°  to  42°. 
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Figure  6.  Relative  vortex  burst  positions  versus  angle  of  attack  for  A  -  50®  to  85® 
(Ref ercnces  1 0, 1 3, 1 4, 1 5). 
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Figure  8.  Burst  point  position  data  taken  on  a  yawed  70®  swept  delta  wing  compared 
lo  Uie  curve  fit  approximation  (Reference  McKeman'^). 
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Figure  9.  Burst  point  position  from  Fig.  5  compared  to  burst  point  function  approximation. 
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Figure  1 1 .  An  effective  sweep  angle  and  effective  angle  of  attack  curve  plotted  over^the 
effective  angle  map  showing  the  stalled  regime  for  the  delta  wing;  a  =  20  , 

A  =  65°,(l)  =  0°to80°. 
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Figure  12.  An  effective  sweep  angle  and  effective  angle  of  attack  curve  plotted  over  the 
effective  angle  map  showing  the  stalled  regime  for  the  delta  wing;  a  =  25°, 

A  =  65°,  (!)  =  0°  to  80“. 


Figure  13.  An  effective  sweep  angle  and  effective  angle  of  attack  curve  plotted  over  the 
effective  angle  map  showing  the  stalled  regime  for  the  delta  wing;  a  =  35°, 

A  =  65°,  <|)  =  0°  to  80°. 
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Figure  14.  An  effective  sweep  angle  and  effective  angle  of  attack  curve  plotted  over  the 
effective  angle  map  showing  the  stalled  regime  for  the  delta  wing;  a  =  40°, 

A  =  65°,  d)  =  0°  to  80°. 
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Abstract 


Introduction 


The  wake  development  behind  an  ogive  cylinder 
was  investigated  using  smoke-wire  visualization 
methods.  Both  steady  and  unsteady  flow  structures 
were  examined  over  a  wide  range  of  incidence  angles 
and  reduced  frequencies.  Results  show  a  qualitative 
similarity  between  the  static  and  dynamic  wake 
development.  At  lower  reduced  frequencies,  the 
dynamic  wake  development  appears  to  be  a  phase- 
shifted  variation  of  static  wake  development.  Finally, 
an  analytic  model  is  proposed  to  explain  the  wake 
development  behind  a  slender  body. 

Nomenclature 


Modern  fighter  designs  arc  characterized  by  long, 
slender  noses  well  forward  of  the  aircraft’s  center  of 
gravity.  This  configuration  complicates  aircraft 
maneuverability  through  the  post-stall  fiighl  regime. 
Vortical  flows  generated  the  extended  forcbody 
produce  non-linear  aerodynamic  loading  which  make 
predictable  control  extremely  difficult.  Missiles 
launched  in  a  highly  dynamic  environment  or  at  large 
incidence  angles  experience  similar  problems. 
Ultimately,  in  order  to  utilize  both  dynamic  and  liigli-a 
environments,  comprehensive  control  paradigms  need 
to  be  developed  based  upon  an  understanding  of  the 
vortical  flow  phenomenology. 
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w 

oe(jj 
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Subscripts 

A 

N 


Maximum  body  diameter 
Oscillation  frequency  (Hz) 

Distance  from  nose  tip  to  vortex 
separation  location 
Reduced  frequency  (tt  f  Woo) 
Distance  from  tip  to  pivot  point 
Overall  model  length 
Reynolds  number,  Woo  Up 
Wake  width 
Oscillation  amplitude 
Static  or  instantaneous  angle  of 
attack 

Mean  pilch  angle 


Aftbody 

Nose 
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The  aerodynamics  associated  with  slender  bodies 
of  revolution  at  high  angles  of  attack  have  been  the 
subject  of  numerous  investigations^  At  sufficiently 
higli  angle  of  attack,  an  asymmetric  vortex  wake  devel¬ 
ops  which  produces  significant  side  force  loading.  At¬ 
tempts  to  quantify  and  control  the  magnitude  and 
direction  of  the  wake  asymmetr)'  have  met  with  varying 
degrees  of  success.  The  reader  is  referred  to  reviews  by 
Hunt^  and  Ericsson  &  Reding^  for  thorough 
discussions  of  high  angle-of-attack  asymmetric  flow 
separation  studies. 

Since  most  maneuvers  involve  rapid  attitude 
changes,  it  is  necessary  to  both  predict  and  control  the 
resulting  wake  struciures.  Smith  and  Nunn"*  found  a 
sliglU  increase  in  normal  force  associated  with  a 
slender  body  pitched  from  0°  to  90°.  They  attributed 
the  effect  to  a  motion-induced  delay  in  wake  transition. 
Gad-el-Hak  and  Ho^  investigated  the  wake  structure  of 
an  ogive-cylinder  oscillating  in  pilch.  The  model  was 
pitched  about  the  mid-chord  while  traversed  through  a 
towing  lank.  Flow  visualization  results  were  reported 
for  a  variety  of  reduced  frequencies  at  angles  of  attack 
ranging  from  0°  to  30°.  They  noted  a  strong  hysteresis 
in  vortex  position  not  unlike  that  seen  for  pitching 
della  wings. 


Monlividas  Reiscnlhal  and  Nagib^  conducted 
smoke  wire  visualizations  of  llie  wake  structure 
generated  behind  a  pitcliing  ogive  cylinder.  A  variety 
of  cylinder  lengths  and  pitch  rates  were  investigated. 
The  dynamic  nose  vortex  separation  locations  data 
collapsed  when  plotted  against  the  pitch  rate  non- 
dimensionlized  by  the  distance  from  the  tip  to  the  pitch 
point  rather  than  body  diameter,  Tliis  definition  of 
non-dimensional  pitch  rate  is  adopted  in  this 
investigation. 

Stanek  and  Visbal^  developed  a  Navier-Stokes 
based  model  of  the  vorticity  field  development  over  a 
3.5  caliber  ogive  forcbody  undergoing  a  0®  to  20®  pitch 
up  maneuver.  Strong  shear  layers  and  their  roll-up 
into  dynamic  stall-like  vortical  structures  was  observed 
as  the  body  decelerated  Subsequent  vorlex/surface 
interactions  produced  secondary  separation  regions, 
ejection  of  vorticity  from  the  surface  and  localized 
regions  of  high  suction. 

The  most  recent  investigations  of  dynamic  slender 
bodies  were  completed  by  Panzer,  Rediniotis  and 
Tclionis^  and  Hoang,  Wetzel  and  Simpson^.  Panzer  et 
aL  mapped  the  velocity  flow  field  behind  a 
hemisphcrc-c>1indcr  using  lascr-Doppler  vclocimctry 
and  seven-hole  probes.  Hoang  ei  ai  measured  the 
surface  pressures  of  a  spheroid  undergoing  a  pitch-up 
maneuver.  The  delay  in  wake  development  was  noted 
in  both  investigations. 

The  wake  stnicture  behind  a  slender  body  at 
various  angles  of  attack  is  quite  complicated.  A 
number  of  very  different  flow  characteristics  evolve  as 
the  pitch  angle  varies  from  0®  to  90®.  Each  change  in 
wake  structure  reflects  a  different  methodology  to  shed 
vorticity  generated  in  the  boundary  layer.  A  more 
complete  understanding  of  these  characteristics  is 
necessary  in  order  to  develop  effective  control  and 
prediction  strategics  related  to  the  underlying  physics 
governing  slender  body  wake  characteristics. 

This  investigation  documented  the  large-scale  flow 
structures  generated  by  an  ogive  cylinder  at  angles  of 
attack  ranging  from  0®  to  90®.  First,  static  flow 
visualization  data  were  collected  to  quantify  the  various 
flow  structures.  These  data  provided  the  metrics  for 
subsequent  dynamic  tests.  Quantitative  comparisons 
were  then  made  between  the  steady  and  unsteady  test 
conditions  in  order  to  evaluate  any  change  in  tlie 
characteristics  of  the  wake  structures  produced. 


Experimental  Facilities 

The  experiments  were  conducted  in  the  University 
of  Colorado's  large  subsonic  (2’x2’)  wind  tunnel.  A 
variable-speed  DC  motor  coupled  to  a  fixed-pitch  fan 
allows  test  section  velocities  from  0  to  65  fps.  One 
entire  side  wall  is  constructed  of  clear  cast  acry'lic  to 
allow  unobstructed  visualization.  A  circular 
polycarbonate  window  in  the  test  section  ceiling  is  used 
for  top-view  visualization.  The  freeslrcam  turbulence 
level  was  previously  measured  at  less  than  0.05%. 

The  model  used  in  all  tests  consisted  of  an  ogive 
nose  attached  to  a  cylindrical  aftbody  (Figure  1)  The 
overall  model  dimensions  were  13.9"  long  and  1.65"  in 
diameter.  The  nose  was  8.4"  long  and  the  aftbody  was 
5.5"  long  providing  a  length  to  diameter  ratios  of  L^/D 
of  5.0  and  La/D  =  3.4,  respectively.  The  nose  tip 
bluntness  ratio,  defined  as  the  local  nose  radius  divided 
by  the  body  diameter  was  7.5®/o.  Tlie  models  were 
tested  at  frccstream  velocities  of  6.6  fps  and  15  fps. 
The  Reynolds  numbers  based  on  body  diameter  were 
5,000  and  12,000.  Pitch  frequencies  were  varied 
between  2.2  Hz  and  5  Hz.  Pitch  amplitude  and  pivot 
point  w'crc  held  constant  at  10°  and  10.5"  from  the  tip, 
respectively. 

A  double  smokc-wirc  located  upstream  of  the 
model  provided  flow  visualization.  The  smoke  wires 
were  stretched  vertically  and  offset  approximately  1 
mm  to  either  side  of  the  model  centerline.  A  mixture 
of  theatrical  fog  fluid  and  2()W50  weight  motor  oil 
coated  on  the  wires  was  ohmically  healed  to  produce 
uniformly  dense  streaklincs  with  a  duration  of 
approximately  2-3  seconds.  The  smoke  wires  could  be 
independently  healed  to  permit  visualizations  of  Just 
one  side  of  the  model. 

The  resulting  vortical  flow  structures  were 
photographed  with  a  35mm  camera  located  outside  of 
the  tunnel  test  section.  Camera  ix)silions  above  and  to 
the  sides  of  the  test  model  produced  orthogonal 
perspectives.  High-inlcnsily  strobes  mounted 
downstream  of  the  test  section  illuminated  the 
strcaklines  with  7  /isce  flash  pulses.  For  the  dynamic 
motion  cases,  the  strobe  flash  was  phase-locked  to 
specific  angles.  Some  photos  were  also  taken  with  the 
camera  positioned  inside  the  test  section  at  a  point 
directly  downstream  of  the  model.  In  this 
configuration,  strobes  were  mounted  on  the  lop  of  the 
test  section  directly  above  tlic  base  of  the  model.  The 
flash  was  constricted  to  a  narrow  plane  to  produce  a 
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light  sheet  which  allowed  examination  of  the  wake 
cross-sectional  properties. 

A  Kodak  Ektapro  1000  higli-specd  video  camera 
was  used  to  capture  the  time  dependent  nature  of 
certain  flow  events.  Flow  events  were  videotaped  at 
500  full-frames  per  second.  The  Kodak  videotaped 
images  were  downloaded  to  standard  VHS  format  tapes 
for  later  analysis. 

Results 

Steady  State 

Photographs  were  taken  from  the  side  and  top  of 
the  test  section  at  model  pitch  angles  ranging  from  0“ 
to  90°.  The  Reynolds  numbers  based  on  maximum 
body  diameter  were  5,000  and  12,000.  The  smoke 
wires  were  heated  independently  during  the  side  view 
photographs  so  that  the  near  and  far  sides  of  the  wake 
(looking  from  the  side)  could  be  differentiated  and 
analyzed  separately.  The  right  and  left  sides  of  the 
wake  structure  are  more  easily  distinguishable  from 
above  the  model,  thus,  both  smoke  wires  were  heated 
simultaneously  for  top  view  photos.  This  latter 
approach  also  produced  symmetric  thermal  heating  in 
the  flow  about  the  model  thereby  reducing  asymmetric 
flow  perturbations  from  being  introduced  upstream. 

Tlie  wake  asymmetry  has  been  shown  to  be  highly 
sensitive  to  model  nose  roll  angle'**.  Fortunately,  the 
wake  tends  to  favor  a  'regular'  state  over  wide  ranges  of 
nose  roll  angle".  The  roll  angle  sensitivity  of  the 
model  used  in  this  study  was  determined  and  the 
regular  states  were  identified.  The  nose  roll  angle  was 
set  at  approximately  the  midpoint  of  one  of  these 
regular  states  and  care  was  taken  throughout  testing  to 
ensure  that  the  nose  roll  angle  remained  at  that  angle. 

Wake  sensitivity  to  minor  variations  in  frecstream 
direction  was  also  investigated.  The  pitch  angle  was 
varied  from  0°  to  90°  at  yaw  angles  of  j8  =  ±2°.  No 
noticeable  differences  in  large-scale  wake  structures 
compared  to  the  zero  sideslip  case  were  evident,  thus 
eliminating  minor  freestream  variations  as  possible 
causes  of  flow  asymmetry  at  the  specific  roll  angle 
investigated. 

Figure  2  is  representative  of  the  wake  structure 
development  as  the  model  angle  of  attack  is  increased 
from  0°  to  90°.  The  first  column  contains  actual 
visualizations  while  the  second  contains  caricatures  to 
clarify  and  highliglit  the  main  flow  features.  Althougli 


the  wake  structure  developed  from  a  static  slender  body 
have  been  extensively  studied,  static  results  are 
presented  here  for  later  comparison  with  dynamic 
results. 

The  photo/caricature  pair  in  the  top  row  is 
representative  of  the  wake  structure  at  very  small 
incidences.  From  a=0°  to  10°,  the  flow  over  the 
model  remains  attached  to  the  surface.  The  only  large- 
scale  flow  disruption  is  the  bluff  body  wake. 

At  approximately  a  =  10°,  a  pair  of  counter¬ 
rotating  vortices  develop  over  the  leeward  side  of  the 
model  (Figure  2,  second  row).  Tliese  nose  vortices 
grow  larger  as  the  incidence  angle  is  increased.  Base 
shedding  is  also  evident. 

Between  a  =  25°  and  a  =  30°,  the  near  side  vortex 
separates  from  the  model  surface  (Figure  2,  third  row). 
Separation  of  the  far  side  vortex  occurs  at 
approximately  35°.  Asymmetric  nose  vortex 
separation  locations  have  been  linked  to^the  onset  of 
significant  asymmetric  side  force  loading  .  The  onset 
of  asymmetric  separation  locations  in  this  investigation 
is  consistent  with  the  previously  documented  data  for 
this  model  geometry. 

Figure  3  is  a  enlargement  of  the  shear  layer 
between  the  separated  near  side  vortex  and  the  model 
surface.  The  shear  layer  clearly  exhibits  a  pronounced 
and  well-behaved  instability.  The  nature  of  this 
instability  is  unknown,  however,  the  frequency 
corresponds  to  periodic  variations  in  the  density  of  the 
freestream  smokelines.  To  the  authors'  knowledge,  litis 
instability  has  not  been  oreviously  visualized. 
However,  Dcgani  and  Ziliac'^  noted  fluctuations  in 
surface  pressure  measurements  in  the  region  of  the 
nose  vortex  separation  point  which  they  attributed  to  a 
shear  layer  instability. 

The  last  photograph  and  caricature  of  Figure  2 
contain  all  the  events  present  in  the  wake  from  a  =  45° 
to  70°.  Downstream  of  the  nose  vortex  separation 
locations,  aflbody  vortices  develop  which  grow  larger 
with  angle  of  attack.  These  aflbody  vortices  exist  over 
about  10°  angle  of  attack  range  before  transforming 
into  periodic  aflbody  shedding.  The  nose  vortex 
streaklines  remain  stable  and  coherent  up  to  a  =  55°. 
From  a  =  55°  to  a  =70°,  the  nose  vortices  are  still 
distinctly  separate  from  the  aflbody  shedding,  however, 
they  also  appear  to  shed  structures  in  a  discrete, 
periodic  fasliion. 
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Beginning  at  a  =70®,  the  nose  vortices  alternate 
bclwccn  two  shedding  methods:  at  some  limes  the 
nose  vortices  are  distinctly  separate  from  the  lower 
wake,  but  at  other  limes  they  merge  with  the  lower 
wake.  This  transition  from  one  type  of  shedding  to  the 
other  occurs  randomly.  Tlie  vortex  transition  region 
extends  up  to  a  =  80®,  at  which  point  the  wake  sheds 
periodically  from  one  side  to  the  other,  much  like 
Karman  vortex  shedding. 

Dynamic 


The  ogive-cylinder  was  pitched  about  the  center  of 
gravity  of  the  model  and  support  system  to  minimize 
inertia  effects.  The  distance  from  the  lip  to  the  center 
of  gravity  was  11”.  All  photographs  were  taken  with 
the  shutter  open  while  the  five  strobe  flash  events, 
phase  locked  to  pitch  motion,  were  collected.  The 
mean  angle  of  attack  was  varied  from  =  10®  to 
=  70®,  liowcver,  the  pitch  amplitude  cases  was  held 
constant  at  10®.  Testing  was  conducted  at  frecslream 
velocities  of  6.5  ft/sec  and  15  fi/sec  and  at  pilch 
frequencies  of  2  Hz  and  5  Hz  to  produce  reduced 
frequencies  of  k  =  0.4,  0.8  and  2.0.  The  non- 
dimcnsionalizing  length  was  the  distance  from  the  tip 
to  the  pivot  point. 

Figure  4  depicts  the  typical  sequence  of  events  in 
the  near  .side  fiow  structure  development  behind  the 
ogive-cylinder  for  the  three  non-zero  reduced 
frequencies  considered.  Static  results  are  presented  for 
comparison.  The  unsteady  results  were  photographed  a 
mean  pitch  angle,  of  30®. 

At  the  minimum  angle  of  attack,  a  =  20®  (top 
row),  the  dominant  feature  at  both  k  =  0.4  and  k  =  0.8 
is  the  nose  vortex  which  separates  from  the  surface 
approximately  halfway  between  the  lip  and  the  base. 
The  wake  appears  similar  at  k  =  0.4  and  0.8,  although 
at  the  higlier  pitch  rale,  the  nose  vortex  appears  more 
turbulent  The  flow  development  at  very  high  pilch 
rales  is  shown  in  the  fourth  column.  Clearly,  the  wake 
development  methodology  is  different  than  at  the  lower 
pitch  rales.  The  shape  of  the  wake  above  the  model 
appears  similar  to  that  of  the  nose  vortices  at  the  lower 
reduced  frequencies,  however,  the  smokelines  have 
been  highly  diffused,  indicating  a  greater  degree  of 
turbulence  relative  to  the  other  cases.  Also,  this  vortex 
appears  to  originate  at  a  point  downstream  from  the 
tip,  unlike  the  other  cases,  where  the  nose  vortex 
appears  to  originate  at  or  very  near  the  tip.  Upstream 
of  the  separation  location,  the  flow  appears  to  be 
completely  attached  to  the  model  surface  as  if  the  nose 
vortex  originates  at  the  separation  location.  For 


convenience,  this  structure  will  be  referred  to  as  the 
separated  nose  vortex.  In  the  steady  state  case,  the 
nose  vortex  is  still  completely  attached  at  this  angle. 

By  a  =  30®  during  the  upstroke  (Figure  4,  second 
row),  the  nose  vortices  have  reattached  along  the  entire 
length  of  the  model  and  remain  so  during  the 
remainder  of  the  upstroke  at  k  =  0.4  and  0.8.  Note  that 
at  in  the  steady  state  case,  the  nose  vortex  has 
separated  at  a  =  30°.  At  k  =  2.0,  the  separated  nose 
vortex  convects  downstream  while  the  attached  flow 
region  continues  to  grow. 

Tlic  third  row  of  photographs  in  Figure  4  depicts 
the  wake  structure  at  the  maximum  angle  of  attack,  a  = 
40®.  In  the  steady  stale  case,  the  nose  vortex 
separation  location  has  moved  closer  to  the  nose  while 
an  aftbody  vortex  is  evident  downstream  of  the 
separation  location.  At  k  =  0.4,  the  wake  appears  very 
similar  to  the  steady  stale  case.  However,  the  nose 
vortex  remains  completely  attached  in  the  k  =  0.8  case. 
The  nose  vortex  does  not  separate  at  this  higher  pilch 
rate  until  a  =  35®  during  the  downslroke  (not  shown). 
At  k  =  2.0,  the  separated  nose  vortex  has  completely 
detached  itself  from  llte  model  surface.  The  flow  over 
the  model  is  now  attached  along  the  entire  surface  with 
the  exception  of  a  small  disturbance  near  the  model 
base. 

At  a  =  30®  during  the  downslroke  (Figure  4),  the 
nose  vortex  separation  location  at  k  =  0.4  and  0.8  is 
closer  to  the  nose  than  in  the  steady  state  case.  At  k  = 
0.8,  the  nose  vortex  appears  to  have  separated  from  the 
middle  portion  of  the  body  first,  ratlicr  from  the  base  of 
the  model  as  was  the  situation  at  k  =  0  and  0.4.  The 
distance  from  the  model  surface  to  the  vortex  core  at  k 
=  0.8  indicates  a  greater  lag  in  the  development  of 
vortex  separation,  relative  to  the  sialic  case,  than  at  k  = 
0.4.  At  k  =  2.0,  a  vortex  develops  along  the  entire 
length  of  the  model  which  will  be  labeled  a  cylindrical 
vortex  as  it  is  roughly  cylindrical  in  shape.  As  the 
model  approaches  the  minimum  angle  of  attack,  the 
cylindrical  vortex  begins  to  convecl  downstream,  thus 
giving  rise  to  the  separated  nose  vortex  present  during 
the  upstroke. 

The  effect  of  mean  pitch  angle  can  be  seen  in 
Figure  5.  The  reduced  frequencies  in  these  photos  are 
again  k  =  0,  0.4,  0.8  and  2.0,  however,  the  mean  angle 
of  attack  is  now  50®.  During  the  upstroke,  the  nose 
vortex  separation  location  moves  farther  from  the  tip  as 
the  reduced  frequency  increases.  Afibody  vortices  and 
aftbody  shedding  are  evident  downstream  of  the 
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separation  location  in  all  cases.  However,  the  aftbody 
wake  at  k  =  0.8  appears  to  be  as  large  as  that  at  k  =  0.4. 
This  is  probably  due  to  tlic  fact  that  the  pitch  point  is 
located  forward  of  the  model  base.  Therefore,  a 
portion  of  the  altbody  actually  encounters  an  increased 
effective  velocity  during  the  upstroke. 

During  the  downstrokc,  the  nose  vortex  separation 
location  at  k  =  0.4  and  0.8  is  close  to  the  nose  than  the 
static  case.  The  aflbody  wake  in  both  cases  appears 
larger  than  in  the  static  case.  At  k  =  2.0,  the  wake 
appears  to  have  evolved  into  a  cylindrical  vortex  just  as 
it  did  at  a,„  =  30®.  Near  the  tip,  a  small  disruption  in 
the  smokelines  is  evident.  This  disruption  later 
develops  into  the  nose  vortex  during  the  upstroke.  The 
detached  nose  vortex  present  during  the  preceding 
upstroke  can  be  seen  just  downstream  of  the  model  mid 
chord. 

Figure  6  shows  the  near  and  far  side  wake 
visualization  for  the  model  at  an  instantaneous  angle  of 
attack  of  50®  and  at  k  =  0.0,  0.4  and  2.0.  Clearly,  the 
separation  locations  of  the  nose  vortices  are  different  in 
all  three  cases,  indicating  the  presence  of  a  non-zero 
side  force.  Further  evidence  of  wake  asymmetry  is 
presented  in  rear  view,  cross-sectional  photographs  of 
the  cross-sectional  wake  geometry  at  a  =  25°  shown  in 
Figure  7.  Tlie  top  photo  depicts  the  static  wake,  while 
the  bottom  photo  is  at  k  =  0.8,  =  30®  during  the 

downstroke.  Tliis  data  show  that,  at  least  for  strictly 
planar  motions,  motion  induced  effects  do  not 
overpower  the  origins  of  asymmetry  and  that 
asymmetric  side  force  alleviation  will  require  attention 
in  the  unsteady  environment  as  well  as  the  steady 
environment. 

Steady  Vs.  Dynamic 

The  detailed  flow  visualization  conducted  for  this 
investigation  allows  for  the  identification  of  all  large- 
scale  structures  shed  into  tlie  wake  of  a  static  and 
pitching  slender  body.  These  structures  can  now  be 
categorized  into  different  regimes  and  plotted  as  a 
function  of  angle  of  attack  as  shown  in  Figure  8.  The 
static  regimes,  indicated  the  shaded  boxes,  are  an 
elaboration  of  previously  defined  wake  regimes  behind 
a  static  slender  body*^. 

At  k  =  0.4  (Figure  8a),  the  regimes  occur  at 
roughly  the  same  angles  of  attack  as  the  static  regimes, 
however,  the  motion  induced  effects  are  manifested  as 
a  shifting  of  the  regime  boundaries.  This  phase- 
shifting  is  to  higlier  angles  of  attack  during  the 
upstroke  and  to  lower  angles  during  the  downstroke. 


Note  that  no  new  regimes  have  developed  at  this 
reduced  frequency  and  that,  with  the  exception  of  the 
shear  layer  instability,  all  regimes  developed  ty  a  static 
slender  body  are  also  generated  by  the  pitching  slender 
body.  Tlie  shear  layer  instability  also  occurs  in  the 
dynamic  instance,  howe\'cr,  its  presence  was  only 
captured  very  briefly  and  so  was  not  noted  on  the 
regime  map.  The  regime  map  for  k  =  0.8  is  similar  in 
nature,  although  the  phase  shifts  are  vary  slightly. 

However,  at  k  =  2.0  (Figure  8b),  two  dynamic 
regimes  emerge:  the  detached  vortex  phenomena  and 
the  development  of  a  c>’lindrical  vortex  during  the 
doivnstroke.  Neither  of  these  phenomena  arc  present 
in  the  static  case,  nor  are  they  evident  at  the  lower 
pitch  rates.  The  existence  of  these  two  regimes  at  very 
high  pitch  rates  demonstrates  the  truly  dynamic  effects 
of  model  motion  on  the  wake  development. 

Nose  vortex  width  obtained  from  the  light  sheet 
photos  is  plotted  in  Figure  9  for  k  =  0.  The  total  width, 
defined  as  the  near  side  vortex  width  plus  the  far  side 
vortex  ividth,  increases  in  a  fairly  linear  fashion  up  to 
a  =  25°.  As  the  angle  of  attack  increases,  the  far  side 
vortex  width  increases  dramatically  to  a  value  equal  to 
about  one  and  a  half  times  the  model  diameter.  The 
near  side  vortex  width  increases  at  approximately  the 
same  rate,  reaching  a  maximum  of  three-quarters  of 
the  body  diameter  at  a  =  35°.  An  interesting  feature  of 
this  data  is  that  the  near  side  vortex  separates  from  the 
surface  when  the  total  vortex  width  approaches  the 
body  diameter.  Unfortunately,  the  synchronization  of 
the  strobe  flash  to  model  motion  prevented  the  use  of 
higli  flash  rates  necessary  to  provide  sufficient 
illumination  for  the  dynamic  cases. 

Nose  vortex  separation  location  provides  an  easily 
measurable  index  of  vortex  strength*'*.  Static  and 
dynamic  separation  locations  of  the  near  side  nose 
vortex  arc  plotted  in  Figure  10.  The  data  are  plotted  as 
a  function  of  the  non-dimensionalized  pitch  cycle  witli 
the  minimum  angle  of  attack  considered  to  be  the 
beginning  of  the  cycle  while  the  maximum  angle  of 
attack  occurs  at  the  mid-point  of  the  cycle.  During  the 
upstroke,  the  separation  location  is  located  further  from 
the  tip  than  the  corresponding  static  location  for  all 
pitch  rates  considered.  At  the  maximum  angle  of 
attack,  the  separation  location  quickly  moves  upstream 
towards  the  model  nose  at  k  =  0.4  and  0.8.  As 
mentioned  previously,  at  k  =  2.0,  the  nose  vortex 
detaches  from  the  model  at  the  maximum  angle  of 
attack.  During  the  downstroke,  the  dynamic  nose 
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vortex  separation  location  is  located  nearer  to  the  nose 
than  the  corresponding  static  location. 

The  nose  vortex  is  presumably  the  principal 
method  of  convection  for  vorticily  generated  upstream 
of  the  separation  point.  Downstream  of  the  separation 
location,  the  vorticity  generated  must  be  removed  from 
the  model  surface  by  some  other  means.  For  this 
particular  geometry,  the  allbody  vortices  emerge  over  a 
relatively  small  angle  of  attack  range  before  developing 
into  allbody  shedding.  For  convenience,  the  aflbody 
vortices  and  allbody  shedding  will  be  grouped  together 
and  referred  to  as  the  allbody  wake.  The  length  of  the 
aflbody  wake  measured  parallel  to  the  model  axis  along 
a  line  one  body  diameter  downstream  of  the  model 
surface  is  plotted  in  Figure  1 1  as  a  function  of  the  pilch 
cycle.  As  expected,  the  aflbody  length  is  inversely 
proportional  to  the  nose  vortex  separation  location;  as 
the  separation  location  moves  away  from  the  tip,  the 
aflbocfy  wake  length  decreases  and  vice  versa. 

Discussion 

The  flow  visualization  and  kinematic  analysis 
serve  to  characterize  the  wake  development  behind  a 
slender  ogive-cylinder.  In  general,  an  increase  in  angle 
of  attack,  either  statically  or  dynamically,  allows  for 
the  formation  of  nose  vortices  which  separate  from  the 
surface  at  a  well-defined  and  highly  repeatable  angle  of 
attack.  Aflbody  vortices  and  aflbody  shedding  become 
significant  structures  at  higher  angles.  Dynamically 
forced  model  motion  phase-shifted  the  boundaries 
between  the  various  regimes,  with  the  amount  of  phase 
shifting  dependent  upon  the  pilch  rate.  At  very  higli 
pitch  rales,  the  truly  dynamic  effects  of  motion  history 
are  sufficiently  exaggerated  to  allow  for  direct 
visualization.  Based  on  the  data  presented,  general 
observations  can  be  made  concerning  the  mechanics  of 
wake  development  behind  a  pitching  slender  body. 

First,  except  for  very  high  reduced  frequencies,  all 
the  data  indicate  that  the  wake  development  behind  a 
(lynamically  forced  pitching  slender  body  is  very 
similar  in  nature  to  the  wake  development  behind  a 
static  slender  body.  The  regime  maps  demonstrate  the 
similarities  as  well  as  the  phase  shifting  which  arises 
from  model  motion.  The  fact  that  dynamic  motion 
elicits  no  new  shed  structures  suggests  that  the  static 
vorticity  convection  mechanisms  are  able  to  modify 
themselves  to  accommodate  the  motion  induced  change 
in  vorticity  generation. 

The  kinematic  behavior  of  the  nose  vortex 
separation  location  hysteresis  loops  is  further  evidence 


of  the  slalic/dynamic  wake  similarities.  The  separation 
location  plotted  against  the  cycle  show  that  the  nose 
vortices  experience  a  “catch-up/overshool”  cycle.  At 
either  extremum  of  the  pitch  c>clc,  the  nose  vortex 
separation  locations  lag  the  steady  state  values.  As  the 
model  begins  to  pilch  up  or  down,  the  separation 
location  catches  up  to  and  then  surpasses  the  sialic 
location.  The  process  repeals  as  the  model  reaches  the 
other  extremum  of  the  cycle.  The  static  separation 
location  seems  to  act  as  an  anchor  for  the  dynamically 
forced  separation  locations.  The  length  of  the  aflbody 
wake  follows  an  opposite,  but  otherwise  similar  trend. 

The  hysteresis  loops  formed  by  the  separation 
location  as  a  function  of  instantaneous  angle  of  attack 
clearly  depict  the  dependence  of  the  dynamic 
separation  location  on  the  static  value.  At  all  mean 
angles  considered,  the  hysteresis  loops  are  centered 
around  the  static  curve.  Similar  trends  were  reported 
by  Gad-el-Hak  &  Ho^  for  pitching  ogive  cylinder  and 
by  Montividas  et  for  an  ogive  cylinder  undergoing  a 
single  pilch  up  or  pitch  down  motion. 

A  second  observation  concerns  the  nature  of  the 
wake  development  at  the  highest  reduced  frequency 
considered.  Any  rapid  changes  in  model  incidence  will 
most  likely  elicit  rapid  changes  in  wake  circulation 
strength.  In  accordance  with  the  laws  of  Helmholtz, 
these  circulation  changes  nuisl  accompanied  by  a 
new  vortex  system.  Figure  12  indicates  this  type  of 
contiguous  vortex  system.  Although  smoke  was 
introduced  on  just  the  near  side,  smoke  entrained  to  the 
far  side  illuminated  the  entire  vortex  ring.  This  vortex 
system  was  formed  during  the  latter  portion  of  model 
pilch  down  with  the  vorticity  contained  within  the 
structure  generated  during  the  down  stroke.  This  is  the 
same  vortex  system  whicli  detached  from  the  model 
near  the  maximum  angle  of  attack. 

Since  the  vortex  ring  develops  during  the  pitch 
down  and  detaches  at  the  end  of  the  upstroke,  it 
appears  that  the  pilch  down  lends  to  increase  wake 
circulation.  In  contrast,  pilch  up  leads  to  a  decrease  in 
wake  circulation. 

Finally,  the  discrete,  or  ‘quantized*  behavior  of  the 
wake  structure  (the  nose  vortex  in  particular)  suggests 
that  the  discrete,  attached  vorlex/wake  structures  can 
sustain  only  a  certain  level  of  vorticity.  Beyond  some 
threshold  value,  the  wake  is  unable  to  support 
additional  vorticily  and  the  wake  structure  is  altered 
accordingly.  One  re-structure  method  allows  the 
vortex  to  separate  from  the  body.  Additional  vorticity 
generated  downstream  of  the  separation  point  is  then 
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convecled  through  a  newly  developed  structure,  in  this 
case  the  aftbody  wake. 

For  the  conditions  where  vortex  separation  docs 
not  occur  (low  static  angles  or  during  pitch  up  at  low 
mean  angles)  additional  vorticity  is  accommodated  by 
an  increase  in  vortex  size.  The  static  cross-sectional 
visualization  indicate  tliat  the  maximum  vortex  size  is 
bounded  by  the  diameter  of  the  body.  As  vorticity 
increases  with  increasing  angle  of  attack,  the  threshold 
limit  is  again  exceeded.  This  effect  is  demonstrated  by 
the  near  side  nose  vortex  separation  when  the  total 
wake  width  approached  the  body  diameter.  Convection 
of  additional  vorticity  was  accomplished  by  near  side 
vortex  separation,  and  the  production  of  an  additional 
vortex  system  which  permitted  both  vortex  structures  to 
grow. 

The  vortex  ring  discussed  earlier  is  also  indicative 
of  the  quantized  nature  of  the  wake  development.  As 
the  model  pitches  down,  the  vorticity  level  in  the 
attached  wake  increases.  At  some  point,  the  attached 
wake  can  no  longer  support  additional  vorticity.  The 
wake  casts  off  the  vortex  ring,  pinned  at  the  nose,  so 
that  the  downstream  wake  can  accommodate  additional 
vorticity. 

Tlie  concept  of  a  critical  vorticity  concentration 
level  above  which  the  vortex  ^stem  must  transition  is 
not  new.  Visser  and  Nelson*^  carefully  measured  the 
flow  field  properties  of  a  delta  wing  vortex  prior  to 
breakdown.  Their  results  led  them  to  postulate  that  the 
circulation  concentration  gives  rise  to  vortex 
breakdown.  Also,  Ng*^  considered  the  possibility  of 
critical  threshold  vorticity  concentration  in  his 
development  of  a  model  for  control  of  delta  wing 
leading  edge  vortices. 

Summary 

Additional  work  is  necessary  to  support  an  analytic 
model  for  threshold  vorticity  driving  wake  structure. 
Dynamic  force  balance  measurements  currently 
underway  at  the  University  of  Colorado  will  provide 
more  definitive  measurements  of  the  effect  of  model 
motion  history  on  the  wake  development.  Also,  hot¬ 
wire  investigations  are  scheduled  which  will  yield  nose 
vortex  strength  measurements.  These  quantitative 
measures  should  provide  further  insight  into  the 
unstea^  development  of  wake  structure  from  pitching 
slender  body  geometries. 
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Figure  2:  Typical  static  wake  development  with 
ncreasing  angle  of  attack. 


Figure  3:  Enlargement  of  shear  layer  instability. 


Figure  1 :  Schematic  of  model  dimensions 


^  '  Figure?:  Light  sheet  visualization  of 

Figure  6:  Near  and  far  side  wake  visualization  at  o  =  50“,  wake  cross-section  at  a  =  30'*  k  =  0 

k  =  0  (top),  0.4  (middle)  and  2.0  (bottom).  (top)  and  0.8  (bottom). 
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Figure  8:  Static  and  Dynamic  (top.k  =  0.4;  bottom,  k  =  2.0)  wake  regimes  as  a  function  of  angle 
of  attack. 
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e  10:  Nose  vortex  separation  location.  Closed  symbols  -  upstroke, 
open  symbols  -  downstroke. 


Figure  11:  Aftbody  structure  length.  Closed  symbols  -  upstroke 
open  symbols  -  downstroke. 


Figure  12;  Photograph  and  caricature  iiiustrating  vortex  ring. 
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ABSTRACT 

An  experimental  investigation  of  vortex 
development  and  interaction  with  downstream  vertical 
surfaces  is  presented.  A  1/48  scale  model  of  an  F-15C 
is  evaluated  using  both  colored  dye  and  fluorescent  dye 
laser-light-sheet  flow  visualization  in  conjunction  with 
hot-film  anemometry.  A  combination  of  side  view, 
plan  view,  and  cross-cut  views  equated  with  cross-cut 
velocity  and  turbulence  data  gives  a  detailed  accounting 
of  the  events  leading  up  to  vertical  tail  buffet. 
Specifically,  it  was  determined  that  the  primary  vortex 
progresses  from  one  originating  only  from  the  flow 
over  the  wing  to  that  emanating  from  the  forebody  and 
the  wing  flowfields.  The  vortex  tracks  downstream, 
encompassing  the  vertical  tails  with  a  highly  turbulent 
flow.  The  maximum  turbulence  intensity  in  the  vicinity 
of  the  vertical  tails  occurs  between  20®  and  25®  angle  of 
attack.  This  corresponds  to  regions  where  the  local 
velocity  magnitudes  are  near  freestream  values. 

INTRODUCTION 

It  is  only  recently  that  fighter  aircraft  have 
extended  their  operational  realm  to  prolonged  periods 
spent  at  high  angles  of  attack.  Coexistent  with  this 
expansion  of  the  viable  flight  regime  is  a  necessitated 
increase  in  aircraft  control  versatility  and  an  expanded 
interest  in  the  vertical  stabilizer  buffet  problem  which 
has  plagued  fighter  aircraft  in  the  past  Previous 
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research  has  studied  the  interaction  between  upstream 
generated  longitudinal  vortices  and  downstream  twin 
tails.^  On  the  F/A-18,  the  source  of  the  vortices  is  due 
primarily  to  bursting  of  vortices  shed  off  the  leading 
edge  extension.2  On  the  F-15,  the  source  is  principally 
due  to  shedding  off  the  wing  leading  edge.^  Forebody 
vortices  are  not  a  primary  contributor  to  the  tail  buffet 
since  the  maximum  buffeting  occurs  in  the  angle  of 
attack  range  from  18®  -  30®  in  both  cases.^ 

Mechanisms  designed  to  control  the  flow 
interacting  with  the  vertical  fins  have  been  explored 
recently.  Hebbar  and  Platzer^  found  that  addition  of 
LEX  fences  successfully  redistributed  turbulence  in  the 
vicinity  of  the  tails  on  a  YF-17  to  higher  frequencies. 
By  drawing  the  dominant  frequencies  away  from  the 
critical  low  fiequencies  of  the  vertical  tails  and 
spreading  them  over  a  wider  band,  the  buffet  effect  was 
reduced.  The  fences  were  also  successful  in 
diminishing  the  magnitude  of  the  unsteady  surfaces 
pressures  interacting  with  the  tails.  Bean  and  Wood^ 
used  tangential  leading  edge  blowing  of  the  leading 
edge  of  a  delta  wing  to  modify  the  fiowfield 
influencing  downstream  twin  tails.  The  blowing  shifted 
both  the  buffet  frequencies  to  higher  values  and 
extended  the  angle  of  attack  corresponding  to  peak 
buffet  intensity.  Increasing  the  blowing  level 
essentially  compressed  the  vortex  cores,  re-establishing 
flow  on  the  delta  wing  from  an  initially  stalled 
condition. 

A  basic  understanding  of  how  vortex 
development  is  interrelated  to  downstream  tail 
buffeting  is  a  necessary  precursor  to  development  of  a 
comprehensive  control  strategy.  Perhaps  the  most 
extensive  study  evaluating  the  flow  physics  involved 
with  vortex-tail  interaction  was  done  by  Washburn, 
Jenkins,  and  Ferman^.  In  their  investigation  of  the 
effect  of  tail  location  on  flow  parameters,  they 
determined  that  the  buffet  load  could  be  reduced  by 
moving  the  tails  toward  the  vortex  core  trajectory. 
They  also  m^ped  out  vortex  trajectories,  magnitude  of 


aerodynamic  force  coefficients,  and  tail  surface 
pressures  as  a  function  of  tail  location. 

.To  date,  there  has  not  been  a  comprehensive 
study  desired  to  track  and  evaluate  the  vortex  from  its 
inception  point  to  its  interaction  with  the  downstream 
vertical  tails  on  an  F-15.  By  thoroughly  understanding 
the  flow  physics  vdnch  lead  to  the  twin  tail  buffet 
problem  on  an  F-15,  the  possibility  of  developing  a 
successful  control  methodology  is  significantly 
enhanced.  This  study  investigates  coireladon  between 
multi-aspect  flow  visualization  and  hot-film 
anemometry  measurements  obtained  in  planes 
perpendicular  to  the  freestream  velocity.  Detailed 
velocity  and  turbulence  values  are  analyzed  in 
conjunction  with  both  vortex  characteristics  and 
trajectories  to  gain  a  thorough  grasp  of  the  specific 
physical  mechanisms  governing  the  twin  tail  buffet 
phenomena  on  an  F-15. 

EXPERIMENTAL  METHODS 

Facilities 

A  qualitative  and  quantitative  understanding 
of  the  flow  physics  governing  flow  over  the  F-15  was 
obtained  using  two  experimental  methods,  flow 
visualization  and  hot-film  anemometry.  Flow 
visualization  was  conducted  in  the  0.38  m  x  0.50  m  test 
section  of  the  Aeronautics  Laboratory’s  water  tunnel  at 
a  velocity  of  7 .6  cm/s,  or  a  chord  Reynolds  number  of 
11,000.  The  hot-film  anemometry  experiments  were 
conducted  in  the  Aeronautics  Laboratory's  open  return, 
low  speed  wind  tunnel  at  the  U.  S.  Air  Force  Academy. 
Die  0.91  m  x  0.91  m  test  section  is  designed  for  use 
with  flow  sensor  measurements,  and  has  freestream 
turbulence  intensity  levels  below  0.5%.  All 

measurements  were  taken  at  a  fieestream  velocity  of 
15.2  m/s,  or  a  chord  Reynolds  number  of  166,000. 

Model  Configuration 

The  models  used  for  the  experiments  were 
1/48  scale  models  of  the  F-15C.  Die  water  tunnel 
model  contained  four  dye  ports  on  the  forebody  and 
one  on  each  wing  glove  located  below  the  apex.  These 
were  used  for  both  colored  and  fluorescent  dye  flow 
visualization.  All  control  surfaces  were  fixed  in  the 
neutral  position,  and  the  aircraft  was  "clean",  with  no 
armament.  Die  model  was  sting  mounted  through  the 
starboard  engine  nozzle.  The  model  engine  inlets  were 
open  so  that  entering  mass  flow  could  exit  through  the 
port  nozzle.  'Die  wind  tunnel  model  was  similar  to  the 
water  tunnel  model  with  the  exception  of  the  engine 


inlets  and  the  sting  mount.  The  wind  tunnel  model  had 
a  solid  fuselage  core  for  added  rigidity  and  therefore 
the  inlets  were  closed.  The  wind  tunnel  model  was 
sting  mounted  along  the  centerline. 

Flow  Visualization 

Visualization  pictures  were  taken  at  5°,  10®, 
15°,  20°,  and  25°  angles  of  attack.  Both  colored  and 
fluorescent  dyes  were  used  for  flow  visualization  in  the 
water  tunnel,  with  dye  supplied  to  the  model  from  a 
containment  reservoir  using  compressed  air.  Die  rate 
of  dye  bleed  was  controlled  by  regulating  the  air 
pressure  through  a  needle  valve.  Dye  was  supplied  to 
two  regions  on  the  model:  1)  two  pair  of  dye  ports  on 
the  forebody,  with  one  pair  on  die  upper  surface  and 
one  pair  on  the  lower  surface,  and  2)  the  engine 
inlet/wing  glove  area,  where  the  dye  ports  were  located 
on  the  outer  side  of  the  inlet  and  just  below  the  apex  of 
the  wing  glove.  An  example  of  the  water  tunnel  model 
mounted  in  the  test  facility  is  given  in  Figure  1. 

An  SLR  35mm  camera  was  used  in 
combination  witfi  a  50mm  lens  to  record  the  rfata  In 
the  colored  dye  cases,  two  1000  Watt  flood  lamps  were 
used  to  illuminate  the  dye.  A  sheet  of  multimode  laser 
light  produced  by  a  12  watt  Argon-Ion  laser  was  used 
to  fluoresce  the  dye  in  the  fluorescent  experiments.  Die 
dye  used  was  a  mixture  of  100  mg  rhodamine  6G  per 
one  liter  of  methanol,  which  fluoresced  bright  yellow 
when  excited  by  the  laser  light  sheet.  The  laser  light 
sheet,  generated  by  an  Aerometrics  LVS  170  fiber  optic 
package,  made  it  possible  to  obtain  cross-cuts  of  the 
vortex  flowfield,  providing  an  excellent  means  for 
examining  the  vortex  structure  in  detail.  Die 
configurations  for  the  fiber  optics  and  the  laser  light 
sheet  are  shown  in  Figures  2  and  3.  Perpendicular 
cross-cuts  relative  to  the  fieestream  flowfield  were 
obtained  for  15°,  20°,  and  25°  angles  of  attack  at 
locations  of  x/1  =  0.62, 0.78,  and  0.93,  where  (x)  is  the 
distance  along  the  fuselage  centerline  and  (1)  is  the 
overall  model  length.  The  locations  coincided 
physically  to  the  engine  nacelle-fuselage  juncture,  just 
upstream  and  just  downstream  of  the  vertical  control 
surfaces.  Die  measurement  stations  are  shown  in 
Figure  4.  Die  laser  light  was  filtered  out  of  the 
photographic  data  using  a  short  pass  filter,  leaving  only 
the  fluorescent  light  for  the  photographic  exposure. 

Hot-Film  Anemometry 

All  hot-film  measurements  were  taken  at  a 
fieestream  velocity  of  15.2  m/s.  Data  was  taken  at 
points  corresponding  to  the  same  non-dimensional 
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fuselage  locations  as  the  flow  visualization.  The  21  by 
21  point  grids  were  set  up  perpendicular  to  the  tunnel 
floor,  coincident  with  the  flow  visualization 
photographs,  with  data  points  taken  at  5  mm  intervals. 
This  resulted  in  441  data  points  for  each  plane.  The 
sampling  frequency  for  the  data  was  128  samples  per 
second.  The  upper  surface  velocity  field  was  surveyed 
using  a  hot-film  probe,  set  at  a  30%  overheat  ratio,  in 
conjunction  with  a  TSI  IFA  100  Intelligent  Flow 
Analyzer.  All  data  acquisition  utilized  a  Keithley  500A 
measurement  and  control  system.  The  hot-film  probe 
was  moved  at  specific  increments  using  a  coinputer 
operated  traverse  coupled  with  the  data  acquisition 
system.  Turbulence  intensity  calculations  were  made 
for  the  tunnel  axis  u-component  of  velocity  only, 
represented  by: 


Turbulence  Intensity  = 


U 


RESULTS 

Colored  Dye  Flow  Visualization 

Colored  dye  flow  visualization  data  were 
acquired  at  angles  of  attack  of  5°,  10°,  15°,  20  ,  and 
25°,  at  zero  degrees  sideslip.  Both  side  view  and  plan 
view  photographs  were  obtained.  In  the  data  presented, 
blue  dye  was  used  to  visualize  the  vortices  forming  at 
the  apex  of  the  wing  glove  and  green  dye  was  used  to 
visualize  the  forebody  flowfield. 

Figure  5  shows  the  flowfield  about  the  F-15  at 
5°  angle  of  attack.  At  this  angle  of  attack  the  majority 
of  the  blue  wing  glove  dye  tracks  below  the  wing.  The 
dye  which  tracks  across  the  top  surface  of  the  wing 
shows  that  the  flow  remains  attached  to  the  top  surface 
of  the  wing  and  no  discernible  vortex  is  generated  at 
the  wing  glove.  The  flow  originating  from  the 
forebody  at  5°  angle  of  attack  separates  from  the 
fuselage  just  aft  of  the  canopy  and  then  tracks 
downstream  along  the  centerline  of  the  aircraft.  The 
forebody  flow  remains  close  to  the  top  surface  of  the 
aircraft  and  does  not  interact  with  the  vertical  tails. 

At  10°  angle  of  attack,  shown  in  Figure  6,  two 
vortical  structures  are  present  at  the  engine  inlet/wing 
glove  locations  on  both  sides  of  the  model  (Figure  6b). 
This  is  evidenced  by  the  two  distinct  regions  of  dye 
parallel  to  the  fuselage  centerline.  The  dominant  vortex 
is  formed  outboard  of  the  engine  inlet  at  the  apex  of  the 
wing  glove,  tracks  downstream  along  the  wing-fuselage 
juncture,  and  is  then  split  by  the  vertical  tail.  This 


vortex  remains  close  to  the  fuselage,  as  seen  from  the 
side  view,  and  breaks  down,  or  spreads  out, 
approximately  at  the  40%  root  chord  location.  The 
second  vortical  structure  forms  on  the  inboard  side  of 
the  engine  inlet,  tracks  downstream  along  the  forebody, 
and  then  coalesces  with  the  primary  wing-glove  vortex 
Approximately  40%  of  the  vertical  tail  span  is 
encompassed  by  the  vortex  wake  at  this  angle  of  attack. 
The  flow  begins  to  separate  from  the  wing  leading 
edge,  forming  a  vortex  which  flows  spanwise.  The 
forebody  flowfield  remains  laminar  and  attached  until 
the  wing  apex  and  then  transitions  to  turbulent  flow. 
As  the  forebody  flow  moves  downstream  the  majority 
of  the  fluid  iq>pears  to  enter  the  inboard  vortex,  which 
subsequently  breaks  down.  The  remainder  of  the  flow 
which  originated  at  the  forebody  enters  into  the  outer 
wing  glove  vortex,  as  seen  in  Figure  6d. 

As  the  angle  of  attack  is  increased  to  15°, 
shown  in  Figure  7,  the  vortices  which  form  about  the 
wing  glove  quickly  breakdown  into  a  large  scale 
turbulent  vortex  wake  structure.  The  flowfield  no 
longer  tracks  directly  downstream,  as  in  the  10°  case, 
but  rather  flows  spanwise  toward  the  wing  tips.  Low 
pressure  areas  atop  the  engine  inlets  and  main  wing 
contribute  to  the  breakdown  of  the  vortex  structure  by 
faycing  it  to  spread  spanwise.  The  large  vortical 
turbulent  wake  grows  in  height  slightly  in  comparison 
to  the  10°  angle  of  attack  case  and  impinges  on  the 
vertical  tails,  immersing  jqjproximately  75  percent  of 
the  vertical  tail  span  in  the  wake.  The  majority  of  the 
flowfield  produced  by  the  forebody  disassociates 
quickly  and  is  entrained  into  the  large  turbulent  vortex 
structure  resident  above  the  main  wing.  However,  as 
evident  in  Figure  7d,  a  segment  of  the  forebody  flow 
does  pass  between  the  votical  tails. 

More  fluid  from  the  forebody  and  wing  glove 
regions  is  entrained  into  the  low  pressure  areas  above 
the  wing  and  engine  inlet  at  20°  angle  of  attack,  shown 
in  Figure  8.  The  vortex  wake  becomes  larger, 
encompassing  most  of  the  wing  planform  as  the  vortex 
tracks  downstream.  At  this  angle  of  attack  very  little  of 
the  fluid  originating  from  the  forebody  passes  between 
the  vertical  tails.  The  forebody  flow  is  drawn  instead 
towards  the  fuselage  surface  just  aft  of  the  canopy  and 
flows  spanwise.  The  large  turbulent  vortical  structure 
resident  over  the  wing  begins  to  separate  from  the  aft 
fuselage  region,  turning  towards  the  favorable  pressure 
gradient  of  the  ftecstream  flowfield,  as  seen  in  Figiw 
8c.  This  results  in  the  entire  span  of  the  vertical  tails 
being  immersed  in  the  turbulent  vortex  flowfield 
(Figures  8a  and  8c). 

At  an  angle  of  attack  of  25°  the  vortex  is  very 
similar  to  that  of  the  20°  case,  except  the  vortex  wake 
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expansion  occurs  further  upstream.  The  turbulent 
vortical  structure  in  the  vicinity  of  the  vertical  tails 
continues  to  expand*  as  seen  in  Figures  9,  and  more 
fluid  from  the  turbulent  vortex  wake  is  now  present 
between  the  twin  vertical  tails. 

Fluorescent  Dye  Flow  VisuaHzation 

Fluorescent  dye  flow  visualization  data  were 
acquired  at  15®,  20®  and  25®  angle  of  attack.  Laser  light 
sheet  cross-cuts  of  the  flowfield  were  obtained  at  each 
angle  of  attack  for  the  three  fuselage  locations.  The 
light  sheet  was  perpendicular  to  the  tunnel  axis  and 
intercepted  the  top  of  the  fuselage  at  x/1  positions  of 
0.62  (position  1),  0.78  (position  2),  and  0.93  (position 
3).  Corresponding  data  is  presented  in  Figures  10 
through  12.  The  photographs  were  taken  from  the  aft 
end  of  the  model  looking  upstream  along  the  tunnel 
centerline.  Fluorescent  dye  was  introduced  into  the 
flowfreld  at  the  wing  glove/inlet  location  as  in  the  blue 
colored  dye  cases. 

The  15°  angle  of  attack  case  is  shown  in 
Figure  10.  At  position  1  (Figure  10a),  the  laser  light 
sheet  cross-cut  reveals  a  vortex  in  the  vicinity  of  the 
wing/fuselage  juncture  and  a  separated  turbulent 
flowfreld  in  the  region  above  the  main  wing,  forming 
off  the  wing  leading  edge.  At  position  2,  the  vortical 
wake  over  the  wing  has  widened  to  envelop  the  entire 
wing  span,  and  has  increased  in  height.  The  turbulent 
vortex  wake  continues  to  widen  in  the  spanwise 
direction  as  it  separates  from  the  wing  trailing  edge,  as 
seen  in  Figure  10c  where  the  laser  light  sheet  crosses 
the  tail  plane.  The  wake  remains  outboard  of  the 
vertical  tails  and  remains  below  the  tip  of  the  tail.  This 
frgure  coincides  with  the  colored  dye  visualization  in 
Figure?. 

Laser  light  sheet  cross-cuts  are  presented  at 
20°  angle  of  attack  in  Figure  11.  At  position  1  the 
flowfield  resembles  the  flowfreld  observed  at  15°  angle 
of  attack,  with  the  exception  that  a  more  clearly  defined 
vortical  flow  is  observed  above  the  inboard  section  of 
the  wing.  The  flow  separates  from  the  wing  leading 
edge  and  enters  the  large  vortical  wake  above  the  wing. 
At  position  2  (Figure  1  lb)  the  size  of  turbulent  vortex 
wake  has  increased  with  respect  to  the  15°  case,  and  at 
position  3  it  is  clear  that  the  region  of  turbulence 
envelopes  the  entire  span  of  the  vertical  tails.  Some  of 
the  vortex  wake  at  position  3  is  now  observed  inboard 
of  the  vertical  tails.  From  visual  real  time  observations 
during  flow  visualization,  it  was  evident  that  the 
complex  structure  observed  at  position  3  is  vortical  in 
nature,  and  remains  so  as  the  angle  of  attack  is 
increased  to  25°. 


The  turbulent  wake  continues  to  grow  in  size 
at  all  three  positions  for  the  25°  angle  of  attack  case,  as 
seen  in  Figure  12.  At  position  1  the  vortex  which  was 
observed  over  the  wing/fuselage  juncture  at  20°  angle 
of  attack  is  now  less  structured  and  continues  to  be  fed 
from  the  fluid  separating  from  the  wing  leading  edge. 
At  position  2,  near  the  wing  trailing  edge,  more  of  the 
vortical  wake  moves  inboard  toward  the  centerline  of 
the  model.  This  confrrms  previous  observations  from 
the  colored  dye  visualization  in  Figure  9.  In  the  tail 
plane,  position  3,  the  vortical  wake  has  substantially 
increas^  in  size  with  the  increased  angle  of  attack,  and 
has  moved  to  a  higher  position.  This  again  compares 
well  with  the  colored  dye  flow  visualization  case  at  25® 
angle  of  attack,  in  which  the  wake  in  the  vicinity  of  the 
empennage  separates  from  the  wing  trailing  edge  and 
bends  toward  Ae  direction  of  the  fieestream  flowfreld. 


Flowfield  Measurements 

Contour  plots  of  turbulence  intensities  above  5 
percent  and  non-dimensional  velocities  were  obtained 
at  angles  of  attack  of  15°,  20°  and  25°  in  the  same  cross 
sectional  planes  as  those  presented  in  the  fluorescent 
dye  photographs.  Consequently,  the  flow  visualization 
photographs  described  in  the  previous  section  can  be 
used  to  aid  in  the  understanding  of  the  quantitative  data 
obtained. 

It  is  important  to  point  out  that  in  this 
preliminary  investigation  of  the  vortex  flowfreld  about 
an  F-15,  only  the  x-component  of  velocity  in  the 
stability  axis  was  measured  with  the  hot-film  probe,  In 
Figures  13,  14  and  15,  the  contour  plots  of  turbulence 
intensity  and  non-dimensional  velocity  are  shown  for 
angles  of  attack  of  15°,  20°  and  25°  respectively. 

In  Figure  13,  it  is  readily  apparent  that  the 
region  of  turbulence  intensity  increases  in  both  width 
and  height  as  the  flowfield  propagates  downstream.  At 
position  1 ,  q)proximately  mid- wing,  a  maximum 
turbulence  of  25%  to  30%  occurs  near  the  surface  of 
fuselage  along  the  model  centerline.  These  high 
turbulence  intensity  pockets  occur  in  regions  of 
decreasing  velocity,  where  the  local  velocity  decreases 
to  60%  of  the  freestream  value.  The  non-dimensional 
velocity  contours  show  the  flow  accelerating  above  the 
wing,  which  is  expected  due  to  the  low  pressure  area 
which  exists  in  that  region.  At  position  2,  near  the 
wing  trailing  edge,  the  regions  of  maximum  turbulence 
intensity  occur  at  60%  -  75%  of  the  half-span,  most 
likely  due  to  the  flowfield  separating  from  the  wing 
trailing  edge.  Again,  a  decrease  in  local  velocity  occurs 
in  the  same  region  where  maximum  turbulence 
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intensity  is  observed.  In  Figures  13e  and  13f  a  line 
drawing  of  the  vertical  tail  is  added  to  the  contour  plot 
because  the  tail  intersects  this  cross-sectional  plane.  A 
turbulence  intensity  of  5%  is  measured  at 
approximately  the  60%  span  location  of  the  vertical  tail 
and  increases  to  25%  at  the  vertical  tail  base.  The  local 
velocities  in  these  regions  vary  between  100%  and  60% 
of  the  freestream  velocity.  Comparing  these  plots  with 
the  flow  visualization  in  Figures  5  and  10,  it  is 
apparent  that  the  40%  turbulence  intensity  measured 
outboard  near  the  70%  half-span  location  is  due  to  the 
vortical  wake  being  shed  off  of  the  wing  trailing  edge. 

Contour  plots  of  turbulence  intensity  and  non- 
dimensional  velocity  are  presented  in  Figure  14  for  the 
20°  angle  of  attack  case.  In  comparison  to  the  15° 
cases,  the  regions  of  high  turbulence  intensity  have 
grown  at  the  respective  positions,  and  again  increase  in 
size  as  the  measurement  location  fo'ogiesses 
downstream.  This  finding  is  well  substantiated  by 
reviewing  the  fluorescent  dye  photographs  at  20°  angle 
of  attack  in  Figure  11,  where  the  growth  of  die  vortical 
regions  with  downstream  distance  is  clearly  visible.  At 
position  1,  the  maximum  turbulent  intensity  has 
increased  to  45%  and  occurs  approximately  at  the  mid¬ 
span  of  the  wing.  This  coincides  with  reduced 
velocities  and  indicates  possible  flow  separation  on  the 
wing.  Above  this  region  of  turbulence,  the  fiowfield  is 
accelerated  due  to  entrainment  of  the  potential 
fiowfield  above  the  wing. 

Closer  to  the  model  centerline,  the 
characteristics  of  the  fiowfield  have  not  changed 
significandy  from  the  15°  angle  of  attack  case.  At 
position  2,  the  maximum  turbulence  intensity  is  55% 
and  occurs  at  40%  wing  half-span,  at  a  location  10%  of 
the  half-span  above  the  wing.  This  appears  to  be  the 
center  of  the  turbulent  vortex  wake  from  examination 
of  Figure  14b.  The  local  velocities  in  this  region  range 
from  30%  to  60%  of  the  freestream  value.  As  the 
turbulent  wake  flows  downstream,  the  region  of 
turbulence  intensity  grows  in  size  according  to  Hgure 
14e.  At  this  location,  the  minimum  turbulence  intensity 
level  plotted  (5%)  now  extends  to  the  tip  of  the  tail  but 
increases  to  25%  at  the  base  of  the  tail.  The  region  of 
maximum  intensity  is  much  closer  to  the  vertical  tail  at 
position  3  than  in  the  15°  angle  of  attack  case  and 
reaches  a  level  of  55%  at  its  peak.  These  observations 
were  confirmed  by  the  flow  visualization  results  in 
Figure  11.  At  this  location,  the  non-dimensional 
velocities  range  from  100%  where  the  turbulence 
intensity  is  15%  to  half  of  that  level  where  the 
turbulence  intensities  are  45%  to  55%. 

In  Figure  15,  turbulence  intensity  contours  and 
non-dimensional  velocity  contours  are  presented  for  an 


angle  of  attack  of  25°.  At  position  1,  the  low 
turbulence  intensity  boundary  of  5%  is  similar  to  that 
observed  in  the  20°  angle  of  attack  case,  but  the 
location  of  maximum  turbulence  intensity  is  located 
further  inboard  than  at  the  lower  angle  of  attack.  The 
maximum  turbulence  intensity  measured  is  45%,  which 
is  the  same  maximum  as  in  the  20°  case,  but  occurs  at  a 
half-span  location  of  approximately  20%  as  opposed  to 
40%  at  the  lower  attack  angle.  The  point  of  maximum 
turbulence  has  also  moved  from  just  above  the  surface 
to  15%  of  the  half-span  above  the  wing's  surface.  The 
local  velocity  corresponding  to  this  region  was 
approximately  95%  of  the  fieestream  velocity. 
Examining  Figure  15b,  there  is  a  region  of  accelerate 
flow  (110%  of  freestream)  that  seems  to  correspond  to 
the  vortex  feeding  sheet  shedding  off  the  wing  leading 
edge.  The  local  accelerations  and  deceleration  are 
difficult  to  explain,  but  could  be  caused  by  discrete 
vortices  in  the  shear  layer  propagating  downstream.^ 

As  the  large  rotating  turbulent  wake  reaches 
position  2,  the  region  of  measured  turbulence  intensity 
has  grown  significantly,  which  corresponds  well  with 
the  laser  light  sheet  photographs  presented  in  Figure  12. 
The  value  of  maximum  turbulence  intensity  of  45%  to 
55%  remains  similar  to  that  observed  in  the  20°  angle 
of  attack  case,  but  the  location  has  shifted  upward  and 
inward  as  it  did  at  position  1.  The  plot  of  non- 
dimensional  velocity  indicates  that  the  local  velocities 
in  this  region  of  maximum  turbulence  intensity  range 
from  40%  to  100%  of  the  freestream.  In  addition,  a 
large  region  of  reduced  velocity  occurs  just  above  the 
wing  trailing  edge.  At  position  3,  in  the  vertical  tail 
plane,  the  contours  of  measured  turbulence  intensity 
now  extend  well  above  the  tip  of  the  tail,  and  a  large 
region  of  45%+  turbulence  intensity  is  observed  along 
the  outboard  section  of  the  tail  and  near  the  tip.  The 
local  velocities  in  this  region  range  from  40%  to  60% 
of  the  freestream  velocity.  On  the  inboard  section  of 
the  tail  turbulence  intensities  of  35%  were  measured 
witii  corresponding  local  velocities  of  60%  to  70%  of 
the  freestream.  From  this  data  is  qipears  that  the 
vortex  wake  "spills  over"  the  tail.  This  is  confirmed  by 
the  flow  visualization  results  in  Figure  12c. 

DISCUSSION 

A  good  understanding  of  the  fiowfield 
development  about  an  F-15  is  obtained  through 
combined  analyses  of  the  colored  dye  flow 
visualization  photographs,  laser  light  sheet  fluorescent 
dye  cross-cut  photographs,  fiowfield  measurements  of 
non-dimensional  velocities,  and  turbulence  intensity 
calculations.  Plots  of  non-dimensional  velocity  used  in 
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conjunction  with  the  plots  of  turbulence  intensity  and 
flow  visualization  allow  for  a  more  complete  analysis 
of  the  flowfield.  For  example,  regions  of  high  velocity 
and  high  turbulence  intensity  are  going  to  impart  a 
greater  momentum  change  than  those  of  high 
turbulence  and  low  velocity  or  low  turbulence  and  high 
velocity.  It  is  this  momentum  change  in  the  fluid  whi^ 
leads  to  the  forces  causing  tail  buffeting. 

As  the  aircraft's  angle  of  attack  is  increased,  a 
complex  vortical  flowfield  develops  between  5“  and 
10®  angle  of  attack.  At  10®  angle  of  attack,  flow 
visualization  data  indicate  that  the  vortical  flows 
originate  from  the  engine  inlet  and  wing  glove  areas. 
There  are  two  distinct  vortices  observed  at  10®  angle  of 
attack;  a  dominate  vortex  structure  which  forms  about 
the  wing  glove,  and  a  secondary  vortex  which  seems  to 
originate  from  the  inboard  comer  of  the  engine  inlet 
As  these  vortices  track  downstream  they  breakdown 
and  combine  into  a  single  turbulent  vortex  wake  and 
impinge  on  the  votical  tails.  The  flow  originating  from 
the  forebody  becomes  turbulent  aft  of  the  canopy  and 
tracks  downstream  between  the  votical  tails  virith  some 
of  the  fluid  entering  the  turbulent  wake. 

At  IS®  angle  of  attack  the  vortex  flowfleld 
changes  considerably.  The  vortices  which  originate 
from  the  wing  glove/inlet  area  quickly  breakdown  and 
are  drawn  spanwise  over  the  wing.  In  addition,  flow 
separates  from  the  wing  leading  edge,  forming  a  large 
vortex  feeding  sheet.  The  fluid  moving  spanwise  from 
the  wing  glove  area  enters  this  feeding  sheet.  A 
classical  delta  wing  vortex  flowfleld  is  established 
which  is  relatively  large  due  to  the  low  sweep  angle  of 
the  wing.  The  flowfield  is  more  comparable  to  a 
broken  down  vortical  wake  structure  than  a  dghtiy 
wound  longitudinal  vortex.  This  vortical  structure 
separates  from  the  trailing  edge  of  the  wing  and  flows 
just  outside  of  the  vertical  tails.  Comparisons  between 
the  laser  light  sheet  photographs  and  velocity 
measurements  show  the  flowfield  is  highly  turbulent 
and  large  velocity  fluctuations  exist  in  ^e  vicinity  of 
the  vertical  tails.  These  conditions  coincide  with  the 
immersion  of  the  tail  in  the  vortical  flow.  This  explains 
at  least  one  source  for  the  vertical  tail  buffeting 

As  the  angle  of  attack  is  increased  to  20®  the 
flowfield  structure  remains  analogous  to  the  15®  angle 
of  attack  case,  but  changes  in  scale.  The  fluid 
originating  from  the  forebody  becomes  entrained  into 
the  separated  flowfield  above  the  wing.  The  vortical 
stracture  grows  in  height  due  to  influence  from  the 
potential  freestream  flow  above  the  model.  The  vortex 
wake  in  the  vicinity  of  the  vertical  tails  now  extends  to 
the  tip  of  the  tail,  and  the  turbulence  intensity  Hatg 
shows  that  the  regions  of  maximum  turbulence  have 


moved  closer  to  the  surface  of  the  tail.  This  presumably 
results  in  higher  dynamic  loading  on  the  vertical  tails  at 
this  angle  of  attack,  and  collaborates  well  with  findings 
from  other  researchers. 

At  25®  angle  of  attack  the  vortical  structure 
continues  to  expand,  and  moves  upward  in  relationship 
to  the  model  as  the  vortex  wake  propagates 
downstream.  The  fluid  originating  from  the  forebody 
enters  into  the  vortex  in  much  the  same  way  as  in  the 
20®  angle  of  attack  case.  Now  the  regions  of  maximum 
turbulence  in  the  tail  plane  are  much  closer  to  the  tail 
surface  and  measured  intensities  of  at  least  5%  extend 
well  above  the  tail  tip.  However,  the  local  velocities  in 
the  vicinity  of  die  tail  are  not  as  high  as  in  the  20®  case. 

Although  tail  loadings  cannot  be  obtained 
from  this  data,  some  conclusions  can  be  hypothesized. 
The  flow  measurement  data  indicates  that  at  20®  angle 
of  attack  the  local  velocities  are  higher  in  the  vicinity  of 
the  tail  than  in  the  25®  case,  but  the  levels  of  turbulence 
are  lower.  This  could  lead  one  to  conclude  that  the 
angle  of  attack  at  which  maximum  momentum 
fluctuations  are  occurring,  and  hence  buffeting,  would 
be  between  20®  and  25®  angle  of  attack.  Komerath,  et 
al,^  found  that  the  level  of  maximum  buffeting 
occurred  at  22°  angle  of  attack.  The  flow  visualization 
data  and  quantitative  data  presented  in  this  paper 
reinforce  those  findings. 


SUMMARY  AND  CONCLUSIONS 

Analysis  of  the  vortical  flowfield  about  an  F- 
15  and  its  interaction  with  the  twin  vertical  tails  was 
accomplished  through  qualitative  flow  visualization 
data  and  quantitative  hot-film  data.  The  combination  of 
the  wind  tunnel  and  water  tunnel  data  provided  an 
excellent  means  for  analyzing  the  complex  flowfield 
existing  about  the  F-15  at  moderate  to  high  angles  of 
attack.  The  major  findings  from  these  experimental 
investigations  follow: 

1)  At  10®  angle  of  attack  a  vortex  forms  at  the  apex  of 
the  wing  glove  outboard  of  the  engine  inlet.  This  vortex 
tracks  downstream,  parallel  to  the  fuselage  centerline 
and  impinges  on  the  vertical  tail  of  the  F-15, 
encompassing  approximately  40%  of  the  vertical  tail 
span.  The  flowfield  above  the  wing  separates  from  the 
wing  leading  edge,  forming  a  vortex  which  flowed 
spanwise  toward  the  wing  tip. 

2)  At  15®  angle  of  attack  the  vortex  originating  from 
the  wing  glove  area  and  a  portion  of  the  forebody  fluid 
is  now  drawn  into  a  large  vortical  structure  resident 
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above  the  wing.  This  vortex  results  from  separation  at 
the  wing  leading  edge,  and  causes  turbulence  intensities 
of  up  to  55%  near  the  wing  leading  edge.  As  this 
turbulent  wake  flows  downstream  it  passes  outboard  of 
the  vertical  tails,  resulting  in  turbulence  intensities  of 
20%  near  the  base  of  the  tail. 

3)  At  20®  and  25®  angle  of  attack  the  basic  structure  of 
the  vortex  flowfield  remains  similar  to  that  of  the  15® 
angle  of  attack  case.  The  significant  difference  in  these 
cases  are  that  the  vortical  flowfield  above  the  wing  has 
increased  in  size  in  both  the  spanwise  and  vertical 
directions  and  more  fluid  originating  at  the  forebody  is 
now  drawn  spanwise  into  the  vortex.  The  vortical  wake 
surrounds  the  entire  span  of  the  vertical  tail  at  20°  angle 
of  attack,  and  extends  above  the  tip  of  the  tail  at  25° 
angle  of  attack.  Turbulence  intensities  near  the  base  of 
the  tail  increase  from  25%  to  35%  as  the  attack  angle  is 
increased.  In  addition,  the  region  of  maximum 
turbulence  intensity  has  moved  much  closer  to  the 
vertical  tail  surface.  For  the  25®  angle  of  attack  case, 
the  maximum  turbulence  region  is  at  the  surface  of  the 
tail.  The  local  velocities  in  the  tail  vicinity  were 
slightly  higher  in  the  20®  case  than  in  the  25®  case. 
Based  upon  these  results  it  could  be  hypothesized  that 
the  maximum  momentum  changes,  and  hence  d3mainic 
loadings  on  the  tail,  occur  between  20®  and  25®  angle 
of  attack. 

The  combination  of  the  flow  visualization  data  and  the 
numerical  data  provided  a  means  of  studying  the 
vortical  flowfield  about  F-15  in  great  detail.  It  was 
determined  that  at  10®  angle  of  attack,  the  vortex  flow 
originates  at  the  wing  glove.  At  angles  of  attack  of  15® 
and  above,  the  primary  vortex  progresses  from  one  fed 
only  from  the  flow  over  the  wing  to  that  fed  from  the 
forebody  and  the  wing  fiowfields.  Ibe  vortex  tracks 
downstream,  encompassing  the  vertical  tails  with  a 
highly  turbulent  flow.  Through  the  continuation  of  this 
study,  to  include  a  more  detailed  survey  of  the  vortical 
wake,  the  ability  of  controlling  the  vortex  flowfield  to 
alleviate  buffeting  can  be  realized. 
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Figaic  1.  Water  Tunnel  Test  Configuration 

^  ■ 


Figure!.  Laser  Li^t  Sheet  Fiber  Optics 


Figure  3.  Laser  Light  Sheet  Configuration 


Position  1 


x/IsO.62 


Position  2 


Figure  4.  Flow  Visualization  and  Hot^FIlm  Measurement  Locations 
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c)  forebody  flow  visualization  -  side  view  d)  forebody  flow  visualization  -  top  view 


Colored  Dye  Flow  Visualization  at  a  =  20' 


a)  position  1,  x/1  =  0.62 


b)  position  2,  x/1  =  0.78 


c)  position  3,  x/1  =  0.93 


Figure  10:  Laser  Light  Sheet  Fluorescent  Dye  Flow  Visualization 
at  a  =  15“,  x/1  =  0.62, 0.78, 0.93. 
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OO 


a)  position  1,  x/1  =  0.62 


b)  position  2,  x/1  =  0.78 


c)  position  3,  x/1  =  0.93 


Figure  11:  Laser  Light  Sheet  Fluorescent  Dye  Flow  Visualization 
at  a  =  20°,  x/1  =  0.62, 0.78, 0.93. 
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O  T 


a)  position  1,  x/1  =  0.62 


b)  position  2,  x/1  =  0.78 


c)  position  3,  x/1  =  0.93 


Figure  12:  Laser  Light  Sheet  Fluorescent  Dye  Flow  Visualization 
at  a  =  25°,  x/1  =  0.62, 0.78, 0.93. 
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Igh  t 


y/b/2  Non -d  Imon  s  Ion  o  I  Span 
a)  turbuience  intensity,  position  1,  x/I  =  0.62 


y/b/2  Non  - d  lm©n  s  Ion  a  I  Span 
c)  turbulence  intensity,  position  2,  x/1  =  0.78 


y /b /2  Non-d Imons  lonal  Span 
e)  turbulence  intensity,  position  3,  x/1  =  0.93 


y/b/2  Non-d  I  mo  ne  lonal  Span 
b)  non-dimensional  velocity,  position  1,  x/1  =  0.62 


y /b /2  Non-d  I mens  lonal  Span 
d)  non-dimensional  velocity,  position  2,  x/1  =  0.78 


y/b/2  Non  - d  Ims n  6  lo n  a  I  Span 
f)  non-dimensional  velocity,  position  3,  x/1  =  0.93 


Igh  t 


y/b/2  Non -d  Imen  c  Ion  o  I  Spon 

a)  turbulence  intensity,  position  1,  x/1  =  0.62 


y-^b/2  Non -d  Imon  s  Ion  o  I  Spon 

b)  non-dimensional  velocity,  position  1,  x/1  =  0.62 


y/b/2  Non -d  Imon  e  Ion  e  1  Spon 

c)  turbulence  intensity,  position  2,  x/1  =  0.78 


y  /b  /2  No  n  “  d  Ifnonc  lono  I  Spon 

d)  non-dimensioiial  velocity,  position  2,  x/1  =  0.78 


V«r  t  Ice  I 


V«r  t  leo  1 1  To  1 1 
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y/b/'2  Non  ~d  Imon  e  Ion  o  I  Spon 

e)  turbulence  intensity,  position  3,  x/1  s  0.93 


y  No  n  —  d  l  fno  n  c  lono  I  Sp'o  n 

f)  non-dimensional  velocity,  position  3,  x/1  =  0.93 


Figure  14:  l^irbulence  Intensity  nnd  Non-dimensionnl  Velocity  Contours 
at  a  =  20%  x/I  s  0,62, 0.78, 0.93. 
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Abstract 

The  development  and  implementation  of 
a  feedback  control  system  for  static  airfoil 
performance  enhancement  is  discussed.  A  two- 
dimensional  NACA-0015  airfoil  is  evaluated 
under  active  control  supplied  by  a  spanwise  slot  at 
the  leading  edge.  A  feedback  control  loop 
couples  the  frequency  at  which  air  is  tangentially 
pulsed  through  the  slot  to  surface  pressure 
coefficients  measured  near  the  airfoil  leading 
edge.  By  optimizing  the  frequency  of  the  control, 
lift  enhancements  of  over  25%  (±4.1%)  and  static 
pressure  stall  delays  of  4°  in  angle  of  attack  are 
achieved.  Frequencies  determined  by  the 
feedback  control  system  correspond  to  Strouhal 
Number  of  0.2  near  the  angle  of  attack  region  of 
stall.  Post  stall  evaluations  show  that  only  very 
low  frequency  (5  Hz,  Stc  =  .04)  control  inputs  are 
successful  in  enhancing  lift. 

Nomenclature  List 

c  Airfoil  Chord 

Cl  L/l/2pU_^  (one  unit  span) 

Stc  Strouhal  Number  based  on  csina 

a  Angle  of  Attack  in  Degrees 

/  Control  Frequency 
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Introduction 

The  importance  of  both  understanding 
and  efficiently  using  active  control  techniques  to 
enhance  the  aerodynamic  performance  of  aircraft 
is  presently  well  accepted.  A  great  deal  of 
research  has  been  conducted  to  evaluate  the 
effectiveness  of  various  active  control  methods 
for  static  airfoil  flow.  The  reader  is  directed  to 
Gad  el  Hak  and  Bushnell'for  a  summary  of  these 
techniques. 

The  research  discussed  in  this  paper  is 
based  on  previous  successful  attempts  to  alter  the 
separation  characteristics  of  a  two-dimensional, 
static  airfoil.  Active  pulse-air  injection,  with  the 
frequencies  of  the  pulse  being  coupled  to  the 
natural  characteristics  of  the  flow,  is  effective  in 
either  delaying  separation  or  reattaching  a 
separating  boundary  layer  and  extending  the  stall 
angle  of  attack.^’^  The  benefits  of  this  type  of 
control  lie  in  the  low  amounts  of  mass  required  to 
effect  the  performance  enhancements  of  the 

airfoil.  By  taking  advantage  of  the  natural 
Kelvin-Helmholtz  instabilities  in  the  separating 
boundary  layer,  or  shear  layer,  researchers  have 
been  able  to  apply  techniques  previously 

developed  for  mixing  layers  with  beneficial 

results.**  By  forcing  the  large  scale  structures 
present  in  the  static  separating  boundary  layer  to 
pair,  the  shear  layer  growth  rate  increases, 
causing  it  to  interact  with  the  airfoil  surface  and 
effectively  reattach  to  the  airfoil. 

Although  a  thorough  physical 

understanding  of  the  flow  has  been  accomplished, 
a  clear  understanding  of  the  parametric  response 
of  the  flow  to  the  active  control  needs  to  be 
established.  In  addition,  a  model  must  be 
developed  to  extend  this  technique  to  a  practical 
application.  The  experimental  analysis  detailed  in 
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this  paper  attempts  to  accomplish  both  of  these 
goals.  A  feedback  control  law  which  can  be 
implemented  as  either  a  state  (modem)  or  classic 
system  has  been  designed  to  allow  both  quick 
analysis  of  various  parameter  combinations  as 
well  as  determination  of  the  optimum  control 
configuration  has  been  developed.  The  system 
uses  an  integrated  control  and  analysis  capability 
to  adjust  active  air  input  depending  on  the  status 
of  surface  pressure  values  on  the  upper  airfoil 
surface.  This  allows  optimization  of  the  surface 
pressure  magnitudes  to  effect  the  greatest 
aerodynamic  performance  on  the  airfoil. 

Experimental  Methods 

Experiments  were  conducted  in  the  Frank 
J.  Seiler  Research  Laboratory's  Low  Speed  Wind 
Tunnel  located  at  the  U.  S.  Air  Force  Academy. 
The  tunnel  has  a  0.91m  x  0.9  Im  test  section,  with 
a  freestream  turbulence  level  below  0.5%.  The 
majority  of  tests  were  conducted  at  a  freestream 
velocity  of  9.1  m/s,  corresponding  to  a  chord 
Reynolds  number  of  80,400.  A  NACA-0015 
model  with  a  60cm  span  and  a  15.2cm  chord  was 
evaluated  under  various  static  angles  of  attack. 
Circular  splitter  plates  30cm  in  diameter  were 
fixed  to  both  ends  of  the  airfoil  to  ensure  two- 
dimensionality.  Pressure  measurements  were 
obtained  along  the  upper  surface  of  the  airfoil 
using  close-coupled  Endevco  8502  miniature 
pressure  transducers.  The  transducers  were 
installed  inside  the  airfoil  to  minimize  any 
possible  lag-time  errors  in  the  surface  pressure 
readings.  For  each  case,  200  data  points  were 
taken  over  a  20°  angle  of  range  along  the  center 
span.  The  error  in  lift  calculations  for  all  data 
was  4.1%.  After  manipulating  the  control 
parameters  using  the  controller,  detailed  pressure 
measurements  were  made  using  a  Concurrent 
Data  Acquisition  System. 

Active  control  of  the  flow  was 
accomplished  by  injecting  air  through  a  spanwise 
slot  located  at  the  airfoil  leading  edge.  The 
system  was  designed  to  inject  air  through  the  slot 
in  a  tangential  manner,  either  using  a  constant 
stream  or  pulsing  the  air  at  frequencies  specified 
by  the  controller.  During  frequency  comparisons, 
the  blowing  momentum  coefficient. 


C|i=(Usiot^/U„^)(slot  width/c)(percent-on  time) 
was  kept  constant  at  0.004.  The  air  injection 
system  employed  two  solenoid  valves  supplying 
air  to  both  ends  of  the  slot.  A  control  circuit  was 
driven  by  the  controller  to  pulse  the  valves 
anywhere  from  0-45  Hz.  The  total  mass  injection 
was  kept  constant.  Details  of  the  air  injection 
system  can  be  found  reference  five^  and  are  also 
shown  in  Figure  1. 


Figure  1.  Detail  of  the  active  air  injection 
process 

The  feedback  control  system  employed  an 
Integrated  Systems,  Inc.  AC- 100  controller, 
allowing  for  eight  input  channels  and  16  analog 
outputs.  The  configuration  used  in  this 
experiment  is  shown  in  Figure  2.  The  AC- 100, 
manipulated  by  a  graphical  user  interface  utility 
called  MATRIXx,  set  the  type  of  active  control 
input  (constant  or  pulsed  air  blowing)  as  specified 
by  the  operator.  The  air  was  then  injected  into  the 
airfoil  flow,  and  the  responses  of  the  surface 
pressure  values  recorded  and  displayed. 
Adjustments  in  the  active  control  input  were  then 
adjusted  accordingly  and  the  process  repeated. 
The  controller  allowed  for  real-time  adjustments 
in  air  pulse  width  and  injection  amplitude, 
keeping  the  total  amount  of  mass  injected  per 
second  constant.  When  the  desired  levels  were 
reached,  data  acquisition  was  accomplished  using 
the  Concurrent  system,  and  lift  values  were 
calculated.  This  allowed  a  clear  determination  of 
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the  control  system's  efficiency  in  enhancing 
airfoil  aerodynamic  performance. 


Figure  2.  Configuration  of  the  feedback 
control  system  using  the  AC-100 

Results  and  Discussion 

The  response  of  the  flow  was  analyzed 
under  various  settings  of  pulsed  air  injection.  For 
the  first  part  of  the  study,  the  feedback  control 
system  was  configured  to  allow  real-time  user 
interface.  Frequency  values  corresponding  to  the 
previously  determined  vortex  passage  frequency 
and  its  subharmonics  were  evaluated.  Figure  3 
shows  the  lift  coefficient  values  for  three 
frequencies  which  were  tested,  compared  to  the 
unforced  case.  Quasi-steady  surface  pressure 
measurements  were  taken  up  to  an  angle  of  attack 
of  20°.  The  airfoil  was  pitched  slowly  at  a  pitch 
rate  of  0.4  deg/sec,  which  corresponded  to  a 
change  of  0.1  deg  per  data  point.  The  results 
closely  resemble  data  taken  previously  by 
statically  moving  the  airfoil,  and  are  therefore 
considered  to  be  quasi-steady.^ 

The  lift  values  for  all  four  curves  are 
close  prior  to  the  near-stall  region  after  10°  angle 
of  attack.  The  5Hz  injection  case  is  slightly 
higher  at  some  points  along  the  curve.  As  the 
airfoil  nears  stall  and  the  flow  starts  to  separate, 
effects  of  the  active  air  injection  become  greater. 
This  is  to  be  expected,  since  the  mechanism  by 
which  this  pulsed  air  control  works  is  by 
interaction  with  natural  instabilities  in  the  flow. 
These  instabilities  are  not  present  in  the  flow  until 
the  flow  begins  to  separate.  Around  8°  angle  of 
attack,  the  5  Hz  controlled  case  shows  a  marked 


increase  of  about  7%  over  the  natural  case.  This 
increase  exceeds  the  ±4.1%  error  for  these 
measurements.  This  lift  enhancement  continues 
until  approximately  15°,  where  the  20  Hz  control 
case  becomes  the  case  with  the  greatest  lift 
coefficients.  The  20  Hz  control  case  continues  to 
dominate  until  about  17°  angle  of  attack,  with 
increases  in  lift  coefficients  in  the  range  of  14- 
16%.  where  it  “exchanges  places”  with  the  5  Hz 
control  case. 


Figure  3.  Lift  Coefficient  vs  Angle  of  Attack, 
natural  and  actively  controlled  cases. 


Figure  4.  Value  for  optimum  control 
frequency  based  on  leading  five  pressure 
transducer  values  for  various  angles  of  attack. 
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Figure  5.  Lift  CoefBcient  vs  Angle  of  Attack, 
natural  and  optimally  controlled  case. 


Forcing  the  flow  at  5  Hz  imparts  the 
greatest  lift  increases  after  17°  angle  of  attack, 
with  increases  in  lift  exceeding  21%.  The  results 
clearly  show  that  the  most  effective  forcing 
frequency  changes  with  angle  of  attack.  This  is 
tobe  expected,  since  past  work  has  show  that  the 
characteristic  flow  frequencies  are  dependent  on 
both  angle  of  attack  and  freestream  velocity.  The 
benefit  of  a  feedback  control  system  lies  in  this 
dependency.  While  it  is  feasible  to  develop  a 
control  model  in  the  wind  tunnel  for  simple  airfoil 
or  wing  configurations,  it  becomes  more  difficult 
in  a  practical  application.  The  feedback  control 
system  can  continuously  adjust  to  changing  flow 
parameters  to  maintaining  top  aerodynamic 
performance. 

The  next  step  in  the  analysis  was  to 
employ  a  simple  algorithm  designed  to 
continuously  evaluate  the  leading  three  pressure 
transducers  on  the  airfoil,  determine  whether  they 
were  becoming  more  or  less  negative,  adjust  the 
pulsed  air  frequency  input,  and  re-evaluate  the 
surface  pressure  conditions.  The  control  input 
was  adjusted  so  that  the  mass  of  air  injected  over 
a  period  of  time  was  constant  while  the  frequency 
was  changed.  The  end  result  was  that  the  pulse 
width  decreased  proportionately  with  increasing 
frequency. 

The  results  of  this  examination  are  shown 
in  Figure  4.  The  optimal  control  frequency,  given 
in  Hz,  decreases  with  increasing  angle  of  attack 
until  18°,  where  it  reaches  a  constant  of  5  Hz. 


When  these  frequencies  are  converted  to  Strouhal 
numbers  using  chord  times  the  sine  of  the  angle  of 
attack  as  the  characteristic  height,  the  values  near 
stall  equate  to  Stc=0.2.  This  value  is  within  the 
range  determined  previously  to  instill  the  greatest 
forcing  response  to  a  downstream  facing  step.^  A 
static  airfoil  positioned  at  a  positive  angle  of 
attack  imparts  a  blockage  in  the  flow  similar  to  a 
downstream  facing  step. 

As  determined  by  Crow  and  Champagne* 
and  Kibens®,  when  forcing  is  tuned  to  the 
fundamental  frequency  or  the  corresponding 
subharmonics,  the  flow  "locks-on"  to  the  control 
frequency.  This  transmits  more  efficient  flow 
control  by  enhancing  vortex  pairing.  Employing 
the  relation  described  by  Katz'°,  Stc  =  /c(sina)AJ 
oo  ,  where  Stc  was  found  equal  to  0.18  -  0.19  for  a 
NACA-0012  airfoil  wake.  When  Lovato  and 
Troutt  compared  these  results  to  their 
characterization  of  the  static  separating  boundary 
layer,  they  found  that  the  lock-on  fundamental 
frequency  at  the  leading  edge  of  a  stall  airfoil 
corresponded  to  a  Stc  of  0.8.^  It  is  reasonable  to 
assume  that  the  structures  present  in  the 
separating  boundary  layer  undergo  pairings  prior 
to  reaching  the  trailing  edge  of  the  airfoil,  and 
thus  one  would  expect  the  lock-on  frequency  to 
be  less  at  the  trailing  edge  than  at  the  leading 
edge.  This  information  indicates  that  during  stall 
the  most  effective  control  frequency  is  the  second 
subharmonic.  Based  on  past  mixing  layer 
research,  this  would  force  every  four  vortices  to 
pair,  thus  reattaching  the  separating  boundary 
layer  to  the  airfoil.  In  the  post  stall  region,  a 
much  lower  frequency  is  needed  to  control  the 
flow,  possibly  a  much  lower  subharmonic. 

The  last  part  of  this  study  investigated  the 
effect  of  pulsing  amplitude  on  the  flow  response 
to  active  control.  While  the  previous  results 
clearly  showed  that  the  optimum  aerodynamic 
performance  enhancements  were  closely  tied  to 
the  natural  frequencies  present  in  the  flow,  the 
dependence  on  amount  of  mass  injected  also 
needs  to  be  determined. 

This  relationship  was  investigated  by 
maintaining  a  constant  injection  frequency  and 
varying  the  mass  injection  magnitude  by  adjusting 
the  “percent-on”  time  of  the  pulse.  Figure  6 
shows  the  results  for  a  pulse  frequency  of  5  Hz  at 
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three  different  injection  magnitudes.  The 
variance  in  magnitude  is  indicated  by  the 
percentage  time  the  air  was  being  injected  per 
second.  Even  an  injection  level  as  low  as  1%, 
corresponding  to  a  of  0.0004,  imparts  some  lift 
enhancement.  The  maximum  lift  increase  of  17% 
occurs  at  about  17°  angle  of  attack.  Performance 
enhancements  decrease  in  effectiveness  soon 
thereafter.  The  control  inputs  at  3%  = 

0.0012)  and  10%  (C^  =  0.004)  are  virtually 
identical,  resulting  in  lift  increases  of  over  25% 
near  the  16°  angle  of  attack  point.  There  is  also 
significant  stall  delays  at  these  levels. 


—  Natural  — 1%  —3%  —10%  ~| 


Figure  6.  Variance  of  Lift  Coefficient  with 
control  magnitude.  Forcing  frequency 
is  5  Hz. 

Conclusions 

The  experimental  analysis  here  has 
demonstrated  the  plausibility  of  a  feedback 
control  system  designed  to  optimize  aerodynamic 
performance  on  a  static,  two-dimensional  NACA- 
0015  airfoil.  The  control  system  has  been 
successfully  proven  to  impart  active,  pulsed-air 
control  on  the  upper  surface  flow  field  of  the 
airfoil  to  enhance  the  lift  characteristics  as  much 
as  25%  (±  4.1%)  in  some  cases.  Current  results 
show  that  the  control  efficiency  is  dependent  on 
forcing  frequency  and  angle  of  attack,  so  by 
altering  the  control  parameters  as  the  flow 
configuration  changes,  the  maximum 


aerodynamic  performance  enhancement  can  be 
obtained. 

Specific  conclusions  can  be  drawn  from 
the  experimental  analysis: 

•  The  most  effective  Strouhal  number  to 
enhance  aerodynamic  performance  is  near  0.2 
in  the  stall  to  post-stall  region. 

•  Controlling  the  flow  using  continuous, 
feedback  adjustments  to  frequency  input 
results  in  lift  enhancements  over  25% 
(±4.1%). 

•  Injecting  pulsed  air  at  blowing  momentum 
coefficients  as  low  as  C^  =  0.0004  can  impart 
significant  aerodynamic  enhancements  to  the 
flow,  even  delaying  the  angle  of  attack  of 
stall. 
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This  paper  examines  the  effects  of  trailing^dge  jet  entrainment  on  the  ^reamw«e 
Although  studies  have  examined  the  effects  of  leading-edge  suction  and  blowing  on  the  burst  location  of  ddto 
wine  vertices,  very  little  research  has  focused  on  the  effects  of  trailing-edge  jet  exhaust  on  ‘he  bui^  location. 
Sg  a  S  “eg  delta  wing  model  in  a  water  tunnel  with  dye  to  mark  the  vortex  core  it  was  po^jble  to  v^hxe 
how  the  locatL  of  the  vortex  breakdown  changes  with  trailing-edge  jet  velocity.  This  re^rch  h^  determined 
that  at  moderate  angles  of  attack  it  is  possible  to  delay  the  burst  locaUon  up  to  18%  of  the  'hjird  hy  ^ 

flow  velocity  from  the  exhaust  ports.  In  addition,  at  higher  angles  of 

asymmetric  separated  vortices  by  reattaching  the  flow  and  moving  the  burst  location  aft  on  the  wing.  ^ 


Nomenclature 

Vr  =  trailing-edge  jet  velocity  ratio,  Vy.JV « 

=  freestream  velocity 
Vjti  -  velocity 

Introduction 

ELTA  wing  vortex  dynamics  and  the  effects  on  a^craft  per¬ 
formance  have  become  an  important  research  topic  in  recent 
years.  The  most  significant  aerodynamic  aspect  of  della  wings  at 
moderate  and  high  angles  of  attack  is  the  formation  of  leading- 
edge  vortices.  These  streamwise  vortices  are  formed  as  the  flow 
separates  at  the  leading  edge  and  rolls  over  the  upper  surface.  At 
high  angles  of  attack,  approximately  half  of  the  lift  generated  on  a 
delta  wing  is  a  result  of  these  vortices  There  are  several  variables 
that  influence  delta  wing  vortex  dynarhics.  Some  of  these  variables 
include  angle  of  attack,  leading-edge  geometry,  wing  thickness, 
sweep  angle,  and  freestream  conditions.  The  most  significant  fea¬ 
ture  that  can  be  influenced  by  these  variables  is  the  bursting  pro¬ 
cess  of  the  vortices. 

Several  theories  governing  vortex  bursting  have  been  pro- 
posed.^"^  As  the  flow  separates,  it  curls  over  the  leading  edge  to 
form  a  well-defmed  vortex.  The  flow  in  core  of  the  vortex  acceler¬ 
ates  as  it  travels  downstream  and  can  reach  values  as  high  as  three 
times  the  freestream  velocity At  some  point  downstream,  the 
core  velocity  will  stagnate,  and  the  vortex  will  burst.  K  the  angle  of 
attack  is  increased  sufficiently,  the  vortex  burst  location  will  move 
upstream  over  the  wing  with  a  highly  turbulent,  wakelike  flow  re¬ 
placing  the  organized  flow  of  the  vortex  and  its  core. 

In  the  study  by  Sarpkaya,^  he  identifies  three  types  of  vortex 
breakdown.  The  most  common  form  of  vortex  breakdown  is  the 
spiral  breakdown.  After  some  distance,  the  flow  along  the  vortex 
core  rapidly  decelerates  until  it  reaches  a  stagnation  poini.  At  this 
bursting  location,  the  core  will  deviate  from  its  centerline  and  form 
a  helixlike  structure  with  a  diameter  much  larger  than  the  original 
core.  The  spiral  structure  will  usually  complete  three  rotations  be- 
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fore  it  breaks  down  into  the  wakelike  flow.'*  The  flow  in  the  center 
of  the  spiral  will  be  in  an  upstream  direction  tetween  the  stagna¬ 
tion  point  and  the  wake  region.  A  second  and  third  type  of  break¬ 
down  identified  by  Sarpkaya  are  the  bubble  and  the  double  helix. 
These  types  of  breakdown  are  not  as  prevalent  as  the  spiral  break- 

down.  .  ,  ,  t 

It  should  be  noted  that  at  low  angles  of  attack  the  burst  location 

may  be  well  behind  the  trading  edge  of  the  wing.  As  angle  of  at¬ 
tack  is  increased  and  the  burst  location  moves  forward  of  the  trail¬ 
ing  edge,  the  lift  generated  by  the  vortices  wUl  decrease.  In  addi¬ 
tion,  at  very  high  angles  of  attack  or  with  high  wing  sweep  angles, 
the  shedding  of  the  vortices  may  not  always  be  symmetrical  and/or 
may  be  separated.  As  a  result  of  this  asymmet^,  the  forces  on  ei¬ 
ther  side  of  the  wing  will  not  be  equal.  The  oscillation  of  the  vorti¬ 
ces  from  one  side  of  the  wing  to  the  other  can  induce  the  phenom¬ 
enon  of  wing  rock. 

To  control  delta  wing  aerodynamics,  several  studies  have 
looked  at  the  effects  of  suction  and  blowing  on  the  leading-edge 
vortices.  The  study  done  by  Visser  et  al.*  examined  the  effects  of 
blowing  near  the  leading  edge  and  back  to  30%  chord.  They  found 
that  the  blowing  was  most  effective  when  located  nearest  the  lead¬ 
ing  edge.  Another  study  was  conducted  by  Magness  et  al.  that  ex¬ 
amined  the  influence  of  suction  as  well  as  blowing.  The  suction 
was  applied  from  a  tube  located  at  a  distance  of  0.95  chord,  mea¬ 
sured  along  the  centerline,  and  placed  out  near  the  leading  edge 
under  the  burst  vortex  core.  The  blowing  was  introduced  to  the 
flow  from  a  port  on  the  surface  of  the  wing  typically  near  the  apex. 
It  should  be  noted  that  most  of  the  previous  studies  involved  small 
amounts  of  mass  flow  out  of  the  ports,  primarily  directed  at  bound¬ 
ary-layer  control.  A  study  by  Roos  and  Kegelman’  examined  the 
effects  of  a  suction  wand  dovmstream  of  the  trailing  edge.  They 
found  the  vortex  burst  location  could  be  changed,  but  for  the  spe¬ 
cific  cases  evaluated,  only  small  changes  in  the  integrated  forces 
on  the  delta  wing  were  measured. 

The  work  presented  in  this  paper  examines  a  different  aspect  of 
delta  wing  vortex  control.  Whereas  previous  studies  have  prima¬ 
rily  used  boundary-layer  blowing/suction,  this  study  concentrates 
on  the  effects  of  large  amounts  of  mass  injection  through  trailing 
edge  jets.  It  is  proposed  that  the  entrainment  from  jets,  or  engines, 
at  the  trailing-edge  of  the  wing  can  be  used  to  control  the  leading- 
edge  vortices.  Thus,  the  influence  of  the  trailing-edge  jets  on  vor¬ 
tex  bursting  and  asymmetry  was  examined  for  a  delta  wing  at  vari- 
ous  angles  of  attack. 

Experimental  Setup 

The  model  for  these  experiments  was  a  60-deg  swept  delta  wing 
fv;rr  n  Thf-  edf^es  of  the  wine  were  sharp  with  a  60-deg 
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Fig.  1  Delta  wing  model  with  60-deg  sweep,  60-deg  underside  lead¬ 
ing-edge  bevels,  and  high-aspect-ratio  rectangular  trailing-edge  Jets. 


Fig.  2  Plot  of  burst  location  vs  for  20-deg  angle  of  attack.  Each 
point  is  an  average  of  four  data  measurements  with  error  bars  for  the 
standard  deviation. 


underside  bevel.  The  jet  exhaust  ports  were  high  aspect  ratio  (8:1) 
rectangular  jets  along  the  trailing  edge  of  the  model.  This  configu¬ 
ration  was  chosen  due  to  geometrical  similarity  to  the  F-117 
stealth  fighter  and  because  the  60-deg  planform  has  been  used  in 
previous  studies.  An  Eidetics  0.46  x  0.46  m  water  tunnel  located 
in  the  United  States  Air  Force  Academy  Aeronautics  Laboratory 
was  used  for  these  studies.  The  water-tunnel  velocity  was  set  at 
7.62  ±  0.3  cm/s.  This  provided  a  Reynolds  number  per  unit  foot  of 
2.07  X  10^.  To  alleviate  wake  perturbations  the  model  was  sting 
mounted  from  the  underside  and  an  internal  cavity,  also  ported 
from ‘the  underside  of  the  wing,  was  used  to  supply  the  mass  flow 
for  the  jets.  To  insure  the  flow  would  exit  from  the  model  parallel 
to  the  freestream,  several  thin-walled  tubes  were  inserted  into  the 
exhaust  ports.  At  this  freestream  velocity,  trailing-edge  jet  to 
freestream  velocity  ratios  were  varied  from  0  to  8,  in  increments 
of  1.0.  The  chosen  velocity  ratios  and  trailing-edge  jet  exit  areas 
are  representative  of  actual  aircraft  parameters. 

Flow  visualization  of  the  delta  wing  vortices  was  accomplished 
by  staining  the  core  of  the  vortex  with  dye.  The  dye  was  injected 
into  the  flow  1.6  mm  from  the  apex  on  the  underside  of  the  leading 
edges.  The  dye  streaklines  were  illuminated  from  the  side  and  pho¬ 
tographed  with  a  35-mm  camera. 

A  transit  was  used  to  insure  that  the  model  was  set  up  in  the  tun¬ 
nel  at  0  deg  yaw,  pitch,  and  roll.  After  the  initial  conditions  were 
set,  the  model  could  be  pitched  by  using  a  crescent  traverse  with  an 
error  of  ±  0.05  deg  in  angle  of  attack.  The  angles  of  attack  used  for 


this  study  were  0, 10, 20,  and  30  deg.  For  most  combinations  of  Vr 
and  angle  of  attack,  the  experiment  was  repeated  so  that  two  sets  of 
data  could  be  compared  for  repeatability.  From  these  sets  of  data, 
measurements  were  made  for  the  burst  location,  core  position  with 
respect  to  the  centerline,  and  the  spreading  angle  of  .the  burst 
vortex. 


Results 

Vortex  Burst  Location 

With  the  model  at  0-deg  angle  of  attack,  the  dye  bled  over  the 
leading  edge  and  spread  into  the  upper  surface  boundary  layer  of 
the  wing.  There  were  no  leading-edge  vortices,  and,  in  addition,  no 
noticeable  flow  disturbances  due  to  the  dye  injection  were  seen. 
The  flow  appeared  to  be  relatively  symmetric  about  the  centerline. 
As  Vr  was  varied  from  0  to  8,  no  changes  were  observed. 

At  10-deg  angle  of  attack  and  Vr  equal  to  0.0,  prominent  lead¬ 
ing-edge  vortices  were  present.  The  burst  location  of  these  vortices 
was  aft  of  the  trailing  edge,  and  again  the  flow  appeared  to  be  rela¬ 
tively  symmetric.  The  bursting  process  followed  that  of  the  spiral 
mode,  as  discussed  earlier.  As  Vr  was  increased  to  8,0,  no  discern¬ 
ible  changes  in  burst  location  were  seen.  It  should  be  noted  at  this 
point  that  the  burst  location  aft  of  the  trailing-edge  is  probably 
dominated  by  the  blunt  finite  trailing-edge  thickness  of  the  model, 
as  it  significantly  affects  the  adverse  pressure  gradient  in  that 
region. 

When  the  model  was  placed  at  20-deg  angle  of  attack,  dramatic 
differences  were  seen  in  the  bursting  of  the  vortices.  The  burst  lo¬ 
cation  vs  Vr  is  plotted  in  Fig.  2.  Each  data  point  shows  the  aver¬ 
aged  measurements  for  the  left  and  right  burst  locations  in  both 
sets  of  data.  At  VJ.  =  0.0,  the  burst  location  had  an  average  value  of 
29%  chord.  As  Vr  was  incremented,  the  burst  location  moved  to 
47%  chord.  Details  of  this  process  can  be  seen  in  Figs.  3a  and  3b, 
which  show  the  leading-edge  vortex  system  at  Vr  =  0.0  and  8.0,  re¬ 
spectively.  In  addition  to  the  burst  location  movement,  these  pho¬ 
tographs  show  relative  symmetry  between  the  left  and  right  vorti¬ 
ces  and  clearly  depict  the  spiral  breakdown  mode. 

As  the  angle  of  attack  was  increased  to  30  deg  and  with  Vr  -  0.0, 
the  vortex  burst  location  moved  to  the  apex  of  the  wing.  In  addi¬ 
tion,  the  completely  burst  vortices  oscillated  from  one  side  of  the 
wing  to  the  other.  This  can  be  seen  in  Figs.  3c  and  3d,  which  were 
taken  at  different  times.  The  period  of  this  oscillation  was  about 
20  s.  Although  similar  oscillations  of  burst  vortices  at  the  apex 
have  been  observed  by  co-workers,  no  documentation  of  this  phe¬ 
nomenon  could  be  found  in  the  literature.  At  this  time  little  insight 
can  be  given  on  the  physical  mechanism  behind  this  process,  al¬ 
though  it  could  be  speculated  that  small  instabilities  in  the  test  fa¬ 
cility,  model  imperfections,  or  model  positioning  could  be  respon¬ 
sible.  As  Vr  was  increased  to  2.0  and  3.0,  it  appeared  that  the 
period  of  oscillation  was  increased  and  the  burst  vortices  began  to 
stabilize.  By  =  5.0,  small  vortex  cores  were  noticed  on  both 
sides  of  the  wing  (Fig.  3e).  Also,  these  vortices  were  stabilized 
without  any  oscillation.  The  bursting  location  was  extended  to 
17%  of  the  chord  and  remained  symmetric  when  Vr  was  increased 
to  8.0  (Fig.  3f). 

Vortex  Core  Position 

In  addition  to  measuring  the  change  of  burst  location  as  Vr  was 
increased,  two  other  variables  were  measured  to  determine  how 
the  vortex  dynamics  were  affected  by  the  trailing-edge  blowing. 
The  vortex  core  position  was  measured  by  the  angle  between  the 
centerline  of  the  wing  and  the  vortex  core.  At  both  10-  and  20-deg 
angle  of  attack,  where  a  vortex  core  could  be  seen,  the  angle  was 
17  ±1  deg.  There  was  no  positive  correlation  between  small 
changes  in  the  vortex  core  position  and  At  30-deg  angle  of  at¬ 
tack  and  with  Vr  greater  than  5.0,  the  vortex  core  angle  was  again 
measured  to  be  approximately  17  deg. 

Spreading  Angle 

The  spreading  angle  is  a  measurement  of  the  divergence  of  the 
burst  vortex  once  it  transitioned  to  a  turbulent  state.  With  spiral 
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Fig.  3  Photographs  of  delta  wing  vortices  where  a)  and  b)  show  burst 
delay  at  20  deg  with  increased  c)  and  d)  show  oscillation  of  burst 
vortices  at  30  deg,  and  e)  and  f)  show  stabilization  and  delay  of  burst 
as  %  increases:  a)  a  =  20  deg,  %  =  0.0;  b)  a  =  20  deg,  Vj.  =  8,0;  c)  a  =  30 
deg,  %  =  0.0,  time  1;  d)  a  =  30  deg,  =  0.0,  time  2;  e)  a  =  30  deg,  %  r 
5.0;  and  f)  a  =  30  deg,  %  =  8.0. 


vortex  breakdown,  the  angle  could  be  measured  using  the  burst  lo¬ 
cation  and  the  edges  of  the  spirals.  These  data  were  only  measured 
at  20  deg  angle  of  attack,  since  it  was  not  clearly  defined  at  10  deg 
where  the  burst  location  was  aft  of  the  trailing  edge,  or  at  30  deg 
where  the  burst  location  was  close  to  the  apex.  The  values  of  the 
spreading  angle  varied  from  24  to  32  deg,  but  again  there  was  no 
positive  correlation  between  these  values  and  Vp. 

Discussion 

The  results  of  trailing-edge  jet  entrainment  on  delta  wing  vortex 
kinematics  demonstrate  that  it  is  possible  to  control  vortex  burst¬ 
ing  and  asymmetry.  By  increasing  to  8.0,  the  burst  location  of 
the  leading-edge  vortices  was  delayed  by  up  to  18%  of  the  chord. 
By  delaying  the  burst  location,  a  more  cohesive  vortex  remains 
closer  to4he  wing  surface  and  thus  should  result  in  an  increase  in 
the  lift.  However,  the  total  amount  of  vorticity  generated  at  the 
leading  edge  does  not  change;  it  is  simply  redistributed  from  a  co¬ 
hesive  vortex  to  a  burst  vortex  that  is  spread  over  a  broader  extent. 
The  distribution  and  rate  of  dissipation  of  the  vorticity  in  the  burnt 
vortex  would  become  factors  in  determining  the  amount  of  lift 
generated.  In  an  aircraft  environment  at  higher  Reynolds  numbers, 
the  leading-edge  vortices  would  exist  in  a  fiilly  turbulent  fashion. 
Again,  the  bursting  process  will  redistribute  the  vorticity  over  a 
much  broader  region,  and  the  ability  to  change  the  burst  location 


could  affect  aircraft  performance  and  controllability.  Without 
detailed  measurements  of  the  vorticity  field,  it  is  difficult  to  specu¬ 
late  at  this  time  as  to  the  actual  effect  on  integrated  lift  or  lateral 

stability.  .  ,  i  ^ 

Perhaps  a  more  important  result  of  jet  entrainment  on  the  lead¬ 
ing-edge  vortices  is  the  ability  to  control  the  asymmetric  nature  of 
flow  at  high  angles  of  attack.  This  was  clearly  demonstrated  at  30- 
deg  angle  of  attack  where  the  flow  was  initially  completely  burst, 
separated,  and  asymmetric.  Not  only  were  the  burst  vortices  reat¬ 
tached  and  extended  over  the  wing,  they  were  stabilized  and  l>e- 
came  symmetric.  This  suggests  a  potential  flow  control  methodol¬ 
ogy  for  the  undesirable  condition  of  wing  rock,  .^though  the 
current  study  did  not  show  any  appreciable  change  in  the  vortex 
core  location  or  spreading  angle  due  to  jet  entrainment,  it  may  be 
possible  to  vary  these  characteristics  by  using  asymmemc  blowing 
or  by  changing  the  jet  location  and/or  using  differential  and  vec¬ 
tored  thrust. 

Modem  delta  wing  aircraft  fly  in  a  dynamic  environment  that 
includes  rapid  changes  in  angle  of  attack  and  coupled  rolling  mo¬ 
tions,  which  have  been  shown  to  affect  vortex  dynamics.  The  po¬ 
tential  to  control  the  resulting  vortex-dominated  flows  with  specif¬ 
ically  programed  thrust/engine  settings  alludes  to  performance 
enhancement  with  little  or  no  penalty.  ^ 

Conclusion 

This  study  examined  the  effects  of  trailing-edge  jet  entrainment 
on  the  flow  over  delta  wings.  Using  a  60-deg  delta  wing  model 
with  high-aspect-ratio  trailing-edge  jets  and  VjfV^  ratios  of  0.0- 
8.0,  significant  events  were  observed. 

It  has  been  shown  that  trailing-edge  mass  injection  can  be  used 
to  move  the  burst  location  of  the  leading-edge  vortices  back  by  up 
to  18%  of  the  chord.  Also,  in  addition  to  Relaying  the  burst  loca¬ 
tion,  trailing-edge  mass  injection  can  be  used  to  reduce  the  asym¬ 
metry,  commonly  found  in  the  leading-edge  vortices  over  delta 
wings  at  high  angles  of  attack.  By  reducing  the  asymmetry,  trail¬ 
ing-edge  jet  entrainment  may  provide  a  technique  to  prevent  unde¬ 
sirable  conditions  like  wing  rock.  The  delayed  burst  location  pro¬ 
vided  by  the  trailing-edge  blowing  could  possibly  result  in  higher 
values  of  lift.  However,  to  determine  the  actual  effects  on  lift,  ex¬ 
periments  employing  force  balances  and/or  surface  pressure  mea¬ 
surements  need  to  be  conducted. 
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INTRODUCTION 

Control  of  unsteady  separated  flow  fields  would  dramatically  enhance  aircraft  agility.  To  enable 
control,  real-time  prediction  of  these  flow  fields  over  a  broad  parameter  range  must  be  realized. 
However,  unsteady  separated  flows  exhibit  complex  spatial  and  temporal  dependencies  which  have 
impeded  accurate,  real-time  flow  field  prediction.  The  present  work  demonstrates  that  neural  network 
models  can  perform  real-time  predictions  of  three-dimensional  unsteady  separated  flow  fields  and 
aerodynamic  coefficients.  These  predictions  are  accurate  over  broad  parameter  ranges  for  novel 
conditions  and  diverse  pitch  motions. 


METHODS 

Unsteady  Surface  Pressure  Measurements 

Unsteady  surface  pressure  measurements  were  performed  in  the  Frank  J.  Seiler  low  speed  wind 
tunnel  using  a  rectangular  wing  bounded  at  the  root  by  a  circular  splitter  plate.  The  wing  had  a  NACA 
0015  cross-section,  0.5  ft  chord  length  and  semi-aspect  ratio  2.0.  Starting  at  0  degrees,  the  wing/splitter 
plate  configuration  was  pitched  at  constant  rate  about  the  quarter  chord  to  a  final  angle  of  60  degrees. 
Seven  nondimensional  pitch  rates  of  0.01,  0.02,  0.05,  0.075,  0.10,  0.15  and  0.20  were  investigated. 
Sinusoidal  pitclung  was  also  examined  for  nondimensional  frequency  0.25.  Test  section  velocity  was  held 
constant  at  9.14  meters  per  second,  corresponding  to  a  chord  Reynolds  number  of  69,000.  Fifteen 
Endevco  8507-2  pressure  transducers  were  close-coupled  to  15  pressure  ports  distributed  along  the  chord 
line  between  0.0  and  0.90  chord.  After  amplification  and  ami-alias  filtering,  these  15  signals  and  a  pitch 
angle  reference  signal  were  sampled  and  digitized  by  a  Masscomp  5500  data  acquisition  system.  A  more 
detailed  description  of  these  methods  has  been  previously  given  (Schreck  and  Helin,  1993a), 

For  all  nondimensional  pitch  rales,  prominent  temporal  variations  in  vortex  structure  and  surface 
pressure  topology  were  apparent.  In  addition,  the  surface  pressure  topologies  exhibited  distihet  spatial 
variations  over  tlie  span.  Dye  flow  visualization  confinned  that,  following  vortex  initiation  near  the  wing 
leading  edge,  vortex  growth  and  convection  were  accompanied  by  radical  three-dimensional  defonnalion 
of  the  unsteady  separated  flow  field  (Schreck  and  Helin,  1993a).  Tlicse  spatially  and  temporally  complex 
three-dimensional  unsteady  separated  flow  fields  fonned  the  experimental  basis  for  the  neural  network 
models  developed  herein. 

Neural  Network  ArchitecluF^aFd  Tiimlng 

The  digitized  unsteady  surface  pressure  records  and  pitch  angle  histories  were  used  to  define  a 
neural  network  model  of  the  tliree-dimensional  unsteady  separated  flow  field.  To  model  the  unsteady 
flow  fields,  a  time-series  paradigm  based  on  the  backpropagation  algorithm  was  utilized.  A  standard 
sigmoidal  activation  function,  l/(l+e"^),  was  used  and  the  input  to  the  network  was  tlie  pitch  angle  (a), 
the  angular  velocity  (da/dt)  and  the  initial  15  recorded  pressure  coefficients  at  tune  (tQ). 

Both  hidden  layers  were  comprised  of  32  units  and  the  output  layer  was  comprised  of  15  units.  Tlie 
targeted  outputs  were  the  15  pressure  coefficients  at  time  (t+At),  as  well  as  coefficients  of  lift,  drag, 
normal  force,  tangential  force  and  moment  about  the  quarter  chord.  Subsequently,  the  time  (l+At) 
network  predictions,  for  each  of  Uie  15  Cp  values,  were  fed  back  as  the  input  to  the  network  throughout 
the  pitch  history.  The  training  set  was  comprised  of  5  nondimensional  pitch  rates,  0.01,  0.02,  0.05,  01.0 
and  0.20,  recorded  at  37.5%  span  outboard  of  the  splitter  plate. 

During  training  the  5  data  sets  were  presented  randomly  with  the  stipulation  that  each  data  set  be 
presented  an  equal  number  of  limes.  The  initial  weights  were  set  randomly  and  training  was  performed 
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until  the  sum-squared  error  was  less  than  0.01  for  each  of  the  5  training  sets.  Thus,  post- training,  the 
network  should  predict  the  entire  surface  pressure  distribution  record  and  aerodynamic  coefficient  records 
given  only  the  pitch  angle  history  and  the  15  pressure  coefficients  at  time  (Iq).  This  type  of  network 
architecture  has  previously  been  described  in  detail  for  predictions  constrained  to  small  time  intervals, 
(At),  for  which  a  linear  equation  system  could  be  derived  from  the  weight  matrices  of  the  neural  network 
^(Faller  et  al.,  1993a;  Schreck  et  al.,  1993b). 


RESULTS 

To  evaluate  the  performance  of  the  neural  network,  the  predicted  surface  pressure  coefficients 
were  compared  directly  to  the  measured  data,  and  the  predicted  aerodynamic  coefficients  were  compared 
to  those  calculated  from  measured  surface  pressures.  Network  inputs  were  the  pitch  angle  (a),  the  angular 
velocity  (da/dt)  and  the  initial  15  recorded  pressure  coefficients  (Cp]-Cp|5)  at  time  (tQ).  Time  (t+At) 
network  predictions,  for  each  of  the  15  Cp  values,  were  fed  back  as  the  input  to  the  network  throughout 
the  pitch  history.  Model  performance  was  evaluated  graphically  by  co-plotting  the  measured  pressure 
profiles  and  the  predicted  data.  Tims,  post- training,  it  was  possible  to  determine  not  only  how  well  the 
models  predicted  the  training  data,  but  how  well  the  models  could  exUapolate  (generalize)  to  both  other 
nondimensional  pitch  rates  and  to  other  types  of  motion  histories. 

Unsteady  Surface  Pressure  Predictions _  _ 

The  graphical*  analyses  for  predicting  both  lire  training  data  as  well  as  for  data  sets  corresponding 
to  novel  pitch  rates  and  motions  are  shown  in  Figs.  1  through  4.  In  all  4  figures,  time-varying  surface 
pressure  at  port  1,  the  leading  edge,  is  at  the  bottom  of  the  figure.  Time-varying  surface  pressure  at  port 
15,  90%  chord  near  the  trailing  edge  of  the  airfoil,  is  at  the  top.  The  ordinate  is  the  pressure  coefficient 
(Cp)  and  the  abscissa  is  nondimensional  time.  The  measured  surface  pressure  data  are  shown  as  a  solid 
line  and  the  surface  pressures  predicted  via  the  neural  network  as  a  dashed  line.  The  magnitudes  of  the 
pressure  coefficients  are  accurate,  but  the  plots  have  been  offset  to  ease  viewing. 

The  analysis  for  nondimensional  pitch  rate  0.02,  1  of  the  5  training  sets,  is  shown  in  Fig.  1. 
Prior  to  suction  peak  occurrence,  the  network  accurately  predicts  surface  pressure  decreases  at  port 
locations  1  through  5.  However,  farther  aft,  at  port  locations  5  through  15,  the  network  failed  to  predict 
small  surface  pressure  fluctuations.  The  network  successfully  predicted  the  time  of  suction  peak 
occurrence  over  the  entire  wing  chord.  However,  suction  peak  magnitudes  were  slightly  underpredicied 
for  most  pressure  port  locations.  Network  predictions  of  constant  pressures  after  suction  peak  occurrence 
were  highly  accurate. 

The  analysis  for  predicting  nondimensional  pitch  rate  0.2,  also  1  of  the  5  training  sets,  is  shown 
in  Fig.  2.  Before  suction  peak  occurrence,  network  predictions  were  highly  accurate  for  all  pressure  port 
locations.  Again,  the  network  failed  to  predict  small  surface  pressure  fluctuations  both  before  and  after 
suction  peak  occurrence.  The  network  accurately  predicted  both  the  time  and  magnitude  of  suction  peak 
occurrence  for  all  pressure  ports.  And,  following  suction  peak  occurrence,  the  network  accurately 
predicted  surface  pressure  magnitudes. 

Figure  3  shows  a  similar  plot  for  nondimensional  pitch  rate  0.15.  This  record  was  not  used 
during  training.  The  neural  network  accurately  predicts  the  initial  surface  pressure  decline  for  the  first  2 
port  locations.  The  initial  stirface  pressure  decline  is  moderately  underpredicted  at  the  remaining  port 
locations.  For  the  leading  5  port  locations,  peak  magnitudes  are  slightly  underpredicted,  and  predicted 
occurrence  times  lag  measured  times  by  approximately  0.5  nondimensional  lime  units.  For  port  locations 
6  through  15,  suction  peak  time  of  occurrence  is  accurately  predicted  and  peak  magnitude  is  moderately 
underpredicted. 

Figure  4  shows  network  predictions  for  sinusoidal  pitch  motion  of  nondimensional  frequency 
0.25.  The  wing  was  pitched  sinusoidally  ±10  degrees  around  a  mean  pitch  angle  of  10  degrees. 
Nondimensional  time  0.0  corresponds  to  the  maximum  instantaneous  pitch  angle  of  20  degrees.  Since 
the  network  was  trained  only  on  constant  rate  pitch  motions,  iliis  record  represents  network  performance 
on  a  distinct  class  of  pitch  motions.  The  results  indicated  that  initial  surface  pressure  increases,  during  the 
pitch-down  phase  of  the  motion,  are  underpredicted  by  the  network.  During  this  phase  of  the  motion  the 
predicted  times  of  occurrence  lend  to  lag  near  die  leading  edge  and  lead  near  the  trailing  edge.  As  the 
wing  reaches  0  degrees,  at  nondimensional  time  7.0,  the  network  predictions  tend  to  collapse  to  the 
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measured  values.  During  the  sinusoidal  pitch  up,  nondimensional  time  7.0  to  12.0,  the  pressure  values  at 
ports  1  to  3  are  underpredicted  but  the  time  of  occurrence  is  accurately  predicted.  The  predicted  values 
fof  ports  4  through  15  are  both  underpredicted  and  tend  to  lead  the  measured  data. 

Unsteady  Aerodynamic  Coefficient  Predictioas 

The  graphical  analyses  for  one  of  the  constant  rate  training  sets  and  the  sinusoidal  pitch  motion 
are  shown  in  Figs.  5  and  6,  respectively.  In  both  figures,  the  abscissa  is  nondimensional  lime  and  the 
ordinate  corresponds  to  the  aerodynamic  coefficients.  The  measured  data  is  shown  as  a  solid  line  and  the 
aerodynamic  coefficients  predicted  by  the  neural  network  as  a  dashed  line.  Again,  the  magnitudes  of  the 
aerodynamic  coefficients  are  accurate,  but  the  plots  have  been  offset  to  ease  viewing.  Starting  at  the 
bottom  of  the  figure  and  proceeding  to  the  top,  the  aerodynamic  coefficients  shown  are  lift  (Cj),  drag 
(Cjj),  normal  force  (C^),  tangential  force  (q)  and  the  pitching  moment  about  quarter  chord  (C^). 

Aerodynamic  coefficient  predictions  for  a  nondimensional  pilch  rate  of  0.02,  1  of  the  5  training 
sets,  is  shown  in  Fig.  5.  Tlie  neural  network  accurately  predicts  these  force  and  moment  coefficients  both 
before  and  after  attainment  of  maximum  magnitude.  For  C\  ,  Q,  q  and  though,  maximum 
magnitude  is  slightly  underpredicted,  and  predicted  lime  of  occurrence  leads  measured  occurrence  time 
by  approximately  0.5  nondimensional  time  units.  Aerodynamic  coefficient  predictions  for  the  sinusoidal 
pitch  motion  are  shown  in  Fig.  6.  For  this  set  ot  conditions,  lift  and  normal  force  are  both  overpredicied 
and  lag  the  measured  data.  The  magnitude  of  and  are  closely  predicted,  but  tend  to  lead  the 

measured  data. 


CONCT.USTONS 

Three-dimensional  unsteady  separated  surface  pressure  fields  were  characterized  in  detail  over  a 
wing  undergoing  constant  rate  pitching  and  sinusoidal  pilch  oscillation.  For  each  data  record,  surface 
pressures  exhibited  prominent  temporal  and  spatial  fluctuations.  Further,  data  record  features  varied  in 
response  to  changes  in  nondimensional  pitch  rate  as  well  as  pilch  motion.  Using  a  subset  of  the  existing 
experimental  data  base,  a  neural  network  model  was  constructed  to  predict  measured  and  novel  unsteady 
surface  pressure  fields.  Network  inputs  were  tightly  constrained,  consisting  only  of  instantaneous  pitch 
angle,  pitch  rate  and  initial  surface  pressure  distribution.  In  spite  of  pronounced  disparities  in  surface 
pressures  and  highly  constrained  network  inputs,  neural  network  model  predictions  were  highly  accurate. 

Overall,  as  measured  by  an  average  deviation  over  lime,  the  results  indicated  that  the  network 
accurately  predicted  the  unsteady  surface  pressures  and  aerodynamic  coefficients  to  within  5%  of  the 
experimental  data.  This  was  verified  graphically  by  co-plotting  the  measured  and  predicted  quantities. 
Consistent  results  were  obtained  for  the  training  sets  as  well  as  for  extrapolation  (generalization)  to  novel 
constant  pitch  rates  and  sinusoidal  pitch  motions.  Tliese  results  clearly  indicate  that  the  neutal  network 
model  can  predict  the  three-dimensional  unsteady  surface  pressure  distributions  based  solely  on  the  pitch 
angle  information.  Tliese  results  also  clearly  demonstrate  that  highly  accurate  real-time  models  of 
unsteady  separated  flow  fields  can  be  developed  using  neural  networks. 
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Fig.  1.  Unsteady  surface  pressures  for  constant  rate  Fig.  2.  Unsteady  surface  pressures  for  constant  rate 

pitcliing,  nondimensional  pitch  rate  0.02.  pitching,  nondimensional  pitch  rate  0.20. 


Fig.  3.  Unsteady  surface  pressures  for  constant  rate  Fig.  4.  Unsteady  surface  pressures  for  sinusoidal 

pitching,  nondimensional  pitch  rate  O.IS.  pitching,  nondimensional  frequency  0.25. 
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Fig.  5.  Unsteady  aerodynamic  coefficients  for  constant  rate  pitching.  Fig.  6.  Unsteady  aerodynamic  coefGcipnts  for  sinusoidal  pitching, 

nondimensional  pitch  rate  0.02.  Coefficients  are  (from  bottom)  lift,  drag,  nondimensional  frequency  0.25,  Coefficients  are  (from  bottom)  lift,  drag, 

normal  force,  tangential  force  and  moment  about  the  quarter  chord.  normal  force,  tangential  force  and  moment  about  the  quarter  chord. 
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